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THE POSSIBILITY OF A PHASE CHANGE IN NUCLEAR MATTER! 


By L. E. H. TRAINoR AND W. R. DIxon 


ABSTRACT 
The possibility is discussed that Livesey’s hypothesis of a constant nuclear 
temperature at low excitations is indicative of a phase change in nuclear matter. 
INTRODUCTION 


In the study of nuclear physics where the basic theory is so little understood, 
it has sometimes been very fruitful to investigate analogies with situations 
in other fields of physics—situations which are presumably more familiar 
and better understood. For sufficiently large nuclei, thermodynamic concepts 
such as heat energy, temperature, and evaporation have been introduced, 
and much attention has been given to statistical models, such as the Fermi-gas 
model and the liquid drop model. These analogies have proved very helpful 
in the description of certain nuclear phenomena. In the present note we wish 
to draw attention to the possibility of another thermodynamic analogy for 
nuclei, namely, that of a phase change in the nuclear matter. This suggestion 
(Trainor and Dixon 1955) arose specifically out of a hypothesis by Livesey 
(1955) that the nuclear temperature is essentially constant for excitations 
below about 10 Mev. The idea of a phase change has been put forward inde- 
pendently by Weisskopf (1953, 1955) on other grounds. 


EVIDENCE FOR A CONSTANT NUCLEAR TEMPERATURE 


In the application of statistical thermodynamics to nuclei, the entropy S 
can be defined as the logarithm of the level density (thermodynamic probabil- 
ity) w, and a nuclear temperature T is introduced in the usual way to relate the 
entropy to the excitation energy F, 


(1) dE =Tds. 

An equation of state relating E and T then determines the entropy S(E) and 
the level density w(£). Simple models give 

(2) S~E 

where n lies between 5; (degenerate Fermi gas) and 1 (constant nuclear 
temperature), depending upon the nuclear model. 


'1Manuscript received September 28, 1955. 

Contribution from the Physics Department, Queen’s University, Kingston, Ontario, and the 
Radiology Section, Division of Applied Physics, National Research Council, Ottawa. Issued as 
N.R.C. No. 3856. 
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Lang and LeCouteur (1954) find that for nuclear excitations in excess of 
about 10 Mev. the experimental data are generally consistent with the model 
of a degenerate Fermi gas. Other experiments suggest, however, that at low 
excitation energies the entropy variation is stronger than the E} variation of 
the Fermi-gas model. For example, the energy distributions of particles 
emitted in (p,m) reactions (Gugelot 1951), in the inelastic scattering of 
neutrons (Graves and Rosen 1953), and in photonuclear reactions (Dawson 
1955, Dixon 1955) can be fitted using a constant nuclear temperature, which 
would imply a linear variation of entropy with excitation energy. (Cohen 
(1953) has suggested an even stronger dependence of S upon E.) Livesey 
(1955) has reviewed a wide range of experimental data, and suggests that 
they may be brought together if one assumes that the nuclear temperature is 
effectively constant below 10 Mev. excitation at a value of the order of that 
given by the Fermi-gas model at 10 Mev., viz. T = 10 A} Mev. 

The evidence for a constant nuclear temperature is by no means conclusive, 
since reliable experiments are difficult to perform and the interpretation of the 
results is far from straightforward. For example, in deriving nuclear tem- 
peratures from the energy distributions of emitted particles, direct effects 
must be excluded. If they are erroneously included, there are too many high- 
energy particles, and the derived temperatures are too high for low energies 
of excitation in the residual nucleus. The subtraction of direct effects is par- 
ticularly important in considering the energy distributions of charged particles, 
since the Coulomb barrier discriminates severely against the evaporation 
process. 

Aside from questions concerned with the applicability of the compound- 
nucleus description (for example, the changing mean free path of nucleons in 
nuclear matter (Feshbach, Porter, and Weisskopf 1954), and the influence of 
direct effects insofar as these cannot be separated out), there are difficulties 
in the determination of nuclear temperatures within the framework of that 
description. Blatt and Weisskopf (1952) write the energy distribution J(e) of 
neutrons evaporated from a compound nucleus as 


(3) I(e) de = const. e o,(e, E) w(emax—e); 


where émax is the maximum energy with which a neutron may be emitted, 
W(€max—€) is the level density in the residual nucleus at an excitation energy 
E = émsx—e, and oa,(e, Z) is the cross section for formation of a compound 
nucleus when a neutron of energy e is incident upon the residual nucleus in 
the state of excitation E. Unfortunately this cross section cannot be measured 
directly, since experimentally a bombarded nucleus is initially in its ground 
state. It is usually assumed that for neutrons a, is independent of both e and E. 
A marked energy dependence of this quantity however could significantly 
affect all conclusions concerning level densities and nuclear temperatures. 
The possibility should also be noted that the experimentally derived tempera- 
ture might depend on the particular reaction being studied, through a selection 
of only certain of the levels of a given energy in the residual nucleus. 
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THE POSSIBILITY OF A NUCLEAR CHANGE OF PHASE 


Although there is no conclusive evidence to support the hypothesis of a 
constant nuclear temperature at excitations below about 10 Mev., it is never- 
theless interesting to speculate upon its possible thermodynamical significance. 
We have suggested (Trainor and Dixon 1955) that Livesey’s hypothesis may 
indicate a first order phase change in nuclear matter, in which energy, entropy, 
and volume increase without change in temperature and pressure. 

If the notion of a phase change in nuclei is a sensible one, two quantities 
are of immediate interest, viz. the latent heat and the temperature at which 
the phase change occurs. To a first approximation one can consider that the 
composition and density of nuclear matter are independent of atomic mass 
number A, from which it follows that both the latent heat per unit volume and 
the transition temperature should be independent of A. The excitation energy 
at which the phase change sets in (which must be low according to the experi- 
mental results), and the excitation energy at which the phase change ceases, 


TEMPERATURE T 
TEMPERATURE T 





EXCITATION ENERGY E ATOMIC MASS NUMBER A 
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Fic. 1. Illustrating possible behavior associated with a nuclear change of phase. (a) Nuclear 
temperature T versus excitation energy E for given mass number A. The dashed curve is for 
the larger value of A. (b) Nuclear temperature T versus mass number A for given excitation 
energy E. The dashed curve is for the larger value of E. 


should both increase linearly with A. This general scheme of things is illustrated 
in Fig. 1(a), where the temperature T is plotted against excitation energy for 
two hypothetical values of A. 

Since the transition temperature has been supposed essentially independent 
of A, and since it can be measured at any energy in the transition region, it 
can be determined experimentally in a variety of ways. It is not so easy to 
measure the latent heat. In principle, the latent heat could be measured for a 
particular nucleus by determining the graph of temperature versus excitation 
energy, but this is difficult to do in practice. However if the latent heat per 
unit volume is independent of A, its value can also be inferred from a graph 
of temperature versus A for a given excitation energy, as illustrated in Fig. 
1(6). For a particular nucleus there is an excitation energy which is just 
sufficient to ‘‘melt’’ that nucleus. With the same excitation energy, nuclei 
heavier than this critical one will be only partially ‘‘melted”’, and the measured 
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temperature will be the transition temperature. The same excitation energy 
in nuclei lighter than the critical one, however, will be sufficient to complete 
the phase change and to heat the nucleus to some point beyond the transition 
temperature. The latent heat per unit volume is then just the ratio of the 
excitation energy used in the experiment to the volume of the critical nucleus. 

Present experimental results however do not permit an accurate determina- 
tion of the latent heat per unit volume. In Fig. 2 the experimental tempera- 
tures found in the photoproton work of Dawson (1955) and in the photoneutron 
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Fic. 2. Experimental values of the nuclear temperature. The solid curve has been drawn 
through the photonuclear results. 


work of Dixon (1955) are plotted as a function of A. The analysis of the 
photoproton results has involved taking into account the direct emission of 
protons, according to the theory of Courant (1951), and the analysis of the 
photoneutron results has involved a small correction for (y, 2m) reactions. 
The temperatures determined from the (, ”’) work of Graves and Rosen 
(1953) and from the (p, ~) work of Gugelot (1951) are also plotted in Fig. 2. 
The data of Graves and Rosen are uncorrected for multiple neutron reactions 
which tend to give too many low-energy neutrons and too low a value for the 
nuclear temperature. In the case of the (p,m) and (n,m’) experiments the 
temperatures shown have been determined for excitation energies of the 
residual nucleus of about 10-12 Mev., while for the photonuclear experiments 
the average excitation energies are about 5-7 Mev. The results of Gugelot 
(1954) on the inelastic scattering of protons have not been included in Fig. 2 
because they are not considered to give reliable values for the temperatures 
owing to the presence of direct effects (Eisberg 1954). 

It is not clear from the experimental results shown in Fig. 2 that the break 
in the curve for the photonuclear results occurs at a smaller value of A than 
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for the particle reactions, as is to be expected from the relative excitation 
energies. It is also not clear whether the transition temperature is actually 
independent of A, as in the simplest picture, or whether it perhaps varies as 
A-+, as in Livesey’s hypothesis. Assuming that the transition temperature is 
independent of A, however, a rough determination of the latent heat per 
nucleon can be made, and one obtains a figure of about 0.1 Mev. per nucleon 


from the photonuclear results. 

In conclusion, it should be emphasized that the above discussion has been 
quite speculative. There is not yet sufficient experimental evidence to deter- 
mine if the idea of a phase change in nuclear matter is useful for the description 


of certain nuclear phenomena. 
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GROWTH TWINS IN GERMANIUM! 


By G. F. BoLiinG, W. A. TILLER,? AND J. W. RUTTER 


ABSTRACT 


The occurrence of growth twins in germanium was investigated as a function 
of the thermal, constitutional, and crystallographic conditions of solidification; 
it was found to be very markedly dependent upon these growth conditions. In 
every case the nucleation of a twin could be traced toa condition of supercooling, 
either absolute or constitutional, during solidification. A new type of substructure 
was also observed to form under certain growth conditions. 


INTRODUCTION 


During recent years, a great deal of interest has been focused upon the 
production of single crystals of the semiconductors germanium and silicon. 
The demands of industry necessitated the immediate production of crystals, 
and consequently very little time was devoted to the understanding of why 
these semiconductors solidified as they did in a manner different from the 
solidification of metal crystals from the melt. 

The most apparent difference between the solidification of semiconductors 
and metals is the frequency of occurrence of growth twins in semiconductors. 
Ellis and Fageant (1954) investigated the orientation relationships in cast 
germanium, and have shown that all major regions in a progressively solidified 
ingot of germanium were related through successive orders of octahedral 
twinning. It was suggested that the entire ingot grew from one nucleus formed 
in initial freezing, and that all subsequent orientations originated through 
nucleation on the surface of the existing solid. Ellis and Treuting (1951) 
showed theoretically that twinned regions can occur with relative ease in the 
diamond cubic structure. This fact is consistent with the existence of twins ina 
cast structure, but it does not explain how or why they occur. 

The investigation of Billig (1954), who studied growth twins in germanium 
and silicon, gives further insight into the phenomenon. He concludes, in 
accordance with earlier investigators, that the observed twinning initiates 
at the solid—liquid interface and is not a solid state transformation. His study 
shows that a crystal attempts to grow with one of its {111} planes aligned as 
closely as possible with its maximum temperature gradient. If necessary, it 
will twin to do this. His observations indicate that twinning is not a haphazard 
event, and that it can be promoted by an unevenness in growth rate. 

Experience at the University of Toronto derived from growing germanium 
crystals from the melt supports the following. If a specimen is grown under 
linear heat flow conditions such that the solid—liquid interface is approxi- 
mately a plane, then a single crystal will usually be obtained. However, when 
the heat flow is not linear, the isotherms in the liquid adjacent to the interface 
become curved. The interface will attempt to adopt the curvature of the 
freezing-point isotherm. Twinning is observed each time before this can be 
accomplished. 


1Manuscript received October 31, 1955. 

Contribution from the Department of Metallurgical Engineering, University of Toronto, 
Toronto, Ontario. 

2Now at Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania. 
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The present work is an attempt to explain the occurrence of twinning in 
terms of thermal, crystallographic, and constitutional conditions. The observa- 
tions are qualitative in nature but lead to an explanation for twinning in 
germanium. 


EXPERIMENTAL PROCEDURES 


Crystals of zone-refined germanium of approximately 10 ohm-cm. resistivity 
were grown from the melt in a horizontal carbon boat using the Chalmer’s 
technique (1953). Approximately linear heat flow was produced by soldering 
a water-cooled copper block to seed crystals of desired crystallographic 
orientation. The crystals were grown under various growth rates and different 
temperature gradients. The solid—liquid interface of each crystal was exposed 
and preserved by separating the solid from the liquid using the accelerated 
decanting technique of Elbaum and Chalmers (1955). The curvature of the 
interface could then be observed, and the orientation of the twin crystals 
determined by X-raying the interface. 

Measurements of the temperature configuration at the solid—liquid interface 
were made by placing a protected thermocouple into the melt. For a desired 
growth rate and imposed temperature conditions, this thermocouple was 
frozen into the crystal and a time-temperature plot of this trial run was 
recorded. The crystal was remelted, and then resolidified under identical 
conditions. Decanting was effected just before the solid—liquid interface 
reached the thermocouple. In this way, both the growth rate and the tempera- 
ture gradient in the liquid at the interface were known from the trial run, and 
the interface was preserved. 

EXPERIMENTAL OBSERVATIONS 


Thermal Effects 

Crystals were grown using the same seed at rates of growth R = 0.16, 
0.5, 1.5 cm./min. It was observed that as the rate of growth increased, the 
interface was further outside the furnace and, therefore, the amount of lateral 
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Fic. 1. Temperature distributions in the germanium near the interface for various growth 
rates. Corresponding interface shapes for these cases are shown. 
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heat loss from the liquid increased. The distance-temperature plots for these 
crystals were as shown in Fig. 1A. It can be seen that, as the rate of growth 
increases, the temperature gradient in the liquid, G, decreases from a positive 
value, through zero, to a negative value; that is, at high growth rates a small 
zone of absolute supercooling exists ahead of the interface. 

In the system that was used, it was impossible to vary R and G independently. 
For example, when R was increased by increasing the rate of furnace movement, 
G decreased. However, both factors could be measured separately, and similar 
experimental results were obtained on different samples provided that the ratio 
G/R was the same, even though R and G individually were different for each. 
As a consequence, it is sufficient to use the parameter G/R to describe the 
general growth conditions. 

A plane interface shape could be maintained at high values of R if G were 
also high, keeping G/R high. With a plane interface, no twins occurred. 
As the lateral heat loss increased, the interface altered shape from a plane to a 
concave parabolic form, as shown schematically in Fig. 1B. Beyond any 
slight curvature, as the ratio G/R decreased, the following events were 
observed: 

(i) Growth twins appeared, a few lamellar twins appearing for small 
curvatures and twin crystals for large curvatures. The density of twinning 
increased inversely as the ratio G/R. 

(ii) All the growth twins extended to at least one outer surface. 

(iii) Under the condition of G approximately zero at the interface, a sub- 
structure appeared on the bottom of the crystal projecting into the liquid. 
This structure is centered by a lamellar twin and is shown in Fig. 9. The 
length of projection ahead of the general interface increased as the ratio G/R 
decreased. This structure is not a simple dendrite because its axis is not the 
dendrite direction. 

(iv) Before shape changes of the interface became too pronounced, the 
twin crystals were examined by X-rays. It was observed that the interface of 
each twin was (within a few degrees) of the {h, h — 1,1} family. These are 
families of planes perpendicular to one of the {111} planes. 


Crystal Orientation Effects 

A single crystal seed was used to grow specimens turned in orientation 
+10° respectively from the seed axis, as shown in Fig. 2. A value of G/R 
was used such that the solid—liquid interface was close to a plane and an 
absence of twins would be expected in the resulting crystal. One specimen 
became aligned within 2° of a [111] direction. It exhibited a regularly spaced 
structure on both sides of the crystal. The structure habit was like twinning 
and appeared lamellar. Some ‘‘fragmentation” appeared in the center of the 
crystal where the lamellae from each side met. 

The other crystal was rotated awav from the particular [111] direction 
mentioned above. The structure appeared again but was more favored on one 
side, as is shown diagrammatically. 
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Fic. 2. Two samples with their growth axes related to <111> directions. A regular struc- 
ture appears as shown. 


Constitutional Effects 

Seed crystals of 35 ohm-cm. resistivity germanium, having a [321] axis, 
were used to grow crystals containing from 0.014 atomic per cent to 0.75 atomic 
per cent additions of gallium. These crystals were grown under values of G/R 
such that no twins would appear in the pure germanium. 

Internally bounded growth twins like those in Fig. 3 began to appear at 
0.043 atomic per cent gallium, and when the gallium content was increased to 
0.13 atomic per cent twin lamellae appeared on two different {111} planes as 
shown in Fig. 4. As the gallium content was increased, the twin density 
increased greatly. 

By scratching the surface and also by casting the shadow of a graticule on 
the surface, it was observed that the twin lamellae projected a very short 
distance ahead of the rest of the interface. 

It was observed that the germanium liquid does not completely wet high- 
resistivity solid germanium. This is evident from the observation that a drop 
of liquid clings to the solid during decanting, as shown in Fig. 5. The gallium 
content increases in a clockwise motion from the upper right in the figure. 
It can be seen that, as the gallium content increases, the contact angle de- 
creases toward zero so that the wetting of the solid by the liquid must increase. 

On the basis of considerations by Walton et al. (1955), it would be expected 
that a transition from a smooth interface to a cellular interface should occur 
as the gallium content is increased. However, the interface exhibits a pox-like 
structure which tends to cluster at twin boundaries, as shown in Fig. 6, long 
after one would expect a cellular interface. Near 0.2 atomic ..«. cent gallium, 
the crystal interface has broken down to exhibit an array ~’ well-developed 
pyramidal cells, as shown in Figs. 7 and 8. 

Crystals were also grown with antimony additions. At 0.1 atomic per cent 
antimony, a huge dendrite which projected 14 cm. into the liquid was preserved ; 
it is shown in Fig. 10. X-ray investigation showed the dendrite direction 
to be the [112] direction, and the upper surface to be a (111) plane. Multiple 
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beam interferometry showed that the upper surface was flat to within 200 A. 
The lower surface of the dendrite was also a (111) plane, but it was very 
irregular. 


DISCUSSION 


The solidification of a semiconducting element like germanium differs 
essentially from the solidification of a metal in that for germanium a metallic 
liquid is in contact with a covalent solid rather than a metallic solid. For 
germanium, there is thus less bonding between solid and liquid than for a 
metal, so that the solid—liquid interface has a higher interfacial energy. 
Because of the directional nature of the bonds, interfaces of different orienta- 
tions will differ greatly in energy and the crystal will find it easier to attach 
atoms to the solid in certain directions rather than in others. 


Thermal and Crystal Orientation Effects 


From the observations it appears that germanium finds it easiest to grow 
in the (h, h—1, 1) directions and that it tries at all times to have an interface 
of the {h, h—/, 1} type. This is consistent with the observations of Billig (1954, 
1955). He finds that germanium prefers to grow by the extension of that set 
of {111} planes most closely perpendicular to the interface. In his study of 
dendritic growth in germanium, Billig shows that dendritic growth in crystals 
with the diamond structure proceeds mainly by the extension of one or more 
sets of the {111} family of planes, and is especially favored along one of the 
(112) directions in these planes. 

The present observations could therefore be accounted for by consideration 
of Fig. 11(a). Consider a germanium crystal growing with AA’ as the interface 
with one {111} plane perpendicular to both AA’ and the page. AA’ will be of 
the family {k, h—1, 1} which is preferred. Then if BB’ B’ represents the freezing 
point isotherm distorted by lateral heat losses, the region of liquid between 
AA’ and BB’’B’ is supercooled. The greatest supercooling is at the outer 
edges A’B’. If this supercooling becomes large enough, a twin will nucleate at 
the outer surface along one {111} plane not perpendicular to the page. The 
interface of the twin can be of the {4, h—/, 1} family of AA’ but at an angle to 
it. The twin crystal can then proceed inwards to absorb the supercooling and 
allow the interface to adjust to the freezing point isotherm. 

To explain the observations for other crystal orientations, Fig. 11(5) must 
be considered. The preferred {h, h—/,1} plane is represented by CC’. It is 
oblique to the direction of general heat flow. DD’’D’ represents the freezing 
point isotherm. If CC’ can exist, there will be supercooling and it will be 
greatest at the edge CD. A twin lamella can nucleate here to favorably alter 
the interface shape, but, as it grows into the crystal, its usefulness disappears 
and another twin must nucleate and so forth. At C’D’, to a smaller extent, 
this cyclic arrangement will also appear. If the thermally determined solid— 
liquid interface is a {111} plane, it will be as far away as possible from the 
preferred type. The appearance of a twin structure would be symmetrically 
favorable at the outer edges and a balanced cycling would appear. 














PLATE | 





Fic. 3.) An internally bounded structure parallel to the trace of a {111} plane appears on 
the decanted interface when some gallium is added = (300). 

Fic. 4. The structure of Fig. 3 appears on two different {111} traces at higher gallium 
concentration. There is a preference for the pock-marked background in some regions (100X). 








PLATE II 





Fic. 5. Germanium interfaces (13). The gallium content: increases for the samples, 
viewed in a clockwise sense from the upper right corner. The “contact angle” decreases indicat- 
ing increased wetting with gallium addition. 

Fic. 6. A twin boundary at the interface shows a preference for the pock marks (1000X). 














PLATE III 





Fic. 7. With large amounts of gallium added, the interface exhibits pyramidal cells (150 X ). 
Fic. 8. A single pyramid obliquely illuminated (500). 





PLATE IV 





Fic. 9. A germanium sample viewed on its lower surface which soliditied in contact with 
the boat. Projecting structures centered by a twin lamella were revealed at the interface by 
decanting (7X). 

Fic. 10. A very flat dendrite grown with antimony added to the germanium (7X). 
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Crystals of unfavorable orientations may then be grown only if a critical 
supercooling for the nucleation of a twin is not reached. This could be achieved 


by controlled temperature gradients in the liquid. 


Solute Effect 

On the basis of considerations concerning the distribution of solute at the 
growing interface by Tiller et al. (1953), it would be expected that even under 
the addition of 0.01 atomic per cent of gallium, a zone of constitutional 
supercooling must exist adjacent to the solid—liquid interface. In metallic 
systems, the interface would alter shape from a smooth interface to a cellular 
interface, as studied by Rutter and Chalmers (1953). However, the interface 
does not develop an array of hexagonal projections, and this may indicate 
either that the growth mechanism is unable to form these projections, or that 
the increase of interfacial energy which must accompany the formation of this 
substructure is too great to allow its formation. 

Thus, constitutional supercooling does exist, resulting in the occurrence of 
twins. As the degree of constitutional supercooling increases, the instability 
is great enough to allow even the nucleation of a second less favorably oriented 
twin system. Constitutional supercooling would most probably be equal over 
the whole interface. The critical supercooling necessary for twin formation 
could then occur first away from the outer edges, giving rise to internally 
bounded twins. This occurrence of internally bounded twins lends support 
to the necessity of supercooling as a prerequisite to twin formation. 

As can be seen from the addition of gallium, the liquid wets the solid much 
better at higher alloy concentrations than at low. This increased bonding 
between liquid and solid indicates that the interfacial energy of the germanium 
has decreased; that is, it is becoming more like a metal. When the alloy 
concentration is high enough, the interface does break down to an array pf 
projections which are similar to the cellular projections previously mentioned. 
The occurrence of well-developed projections which exhibit a pyramidal 
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form indicates either that the surface energy has become low enough to allow 
the transition, or that the growth mechanism has changed. It is difficult on the _ 
basis of few experiments to separate these two effects. Similar structures were 
observed indirectly by Ellis (1955) in the growth of germanium-silicon alloys. 


CONCLUSIONS 


(i) The growth mechanism is such that the germanium exhibits a strong 
tendency to grow in the (h, h—1, 1) directions. 

(ii) The occurrence of twins requires the existence of supercooling, either 
constitutional or absolute. The nucleation of twins requires a certain minimum 
of supercooling for their occurrence. The frequency of twinning increases as 
the supercooling increases. 

(iii) The interfacial energy of solid germanium against liquid germanium 
was lowered by gallium additions. 

(iv) At a certain degree of constitutional supercooling, a smooth interface 
becomes unstable with respect to an interface which exhibits a regular array 
of pyramidal projections. 
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THE PHOTODISINTEGRATION OF OXYGEN NUCLEI 
INTO FOUR ALPHA-PARTICLES! 


By W. K. Dawson? Anp D. L. LIVESEy?® 


ABSTRACT 


A total of 323 four-prong stars found in photographic plates exposed to brems- 
strahlung of maximum energy 70 Mev. have been identified as due to the dis- 
integration of O"* nuclei into four alpha-particles. The excitation function of the 
reaction exhibits pronounced peaks at 22.5 Mev. and 32.5 Mev. as well as 
maxima at intermediate energies, and the integrated cross sections for various 
energy ranges have been estimated. The mechanism of the process has been 
investigated and most of the disintegrations proceed via known excited levels in 
C® and various states of Be’. 


1. INTRODUCTION 


When a photographic plate is exposed to high energy X-rays many types of 
multiple nuclear disintegration occur, leading to the production of identifiable 
stars when the light nuclei of the gelatin are involved. Several authors have 
published accounts of the photodisintegration of C nuclei into three alpha- 
particles and of O'* into four alpha-particles. The following properties of these 
reactions have been established with some degree of certainty: 

(a) the excitation functions exhibit a series of narrow resonances grouped 
into larger peaks which alone are detected when the energy resolution is not 
of a high order (see Goward and Wilkins 1952), 1953), Livesey and Smith 
(1952); the same kind of behavior is shown by the (vy, m) reactions in these 
nuclei according to Katz et al. (1954) and Spicer et al. (1954); 

(b) the mode of disintegration of these even-even light nuclei is predomi- 
nately via excited levels in the intermediate nuclei C” and Be® according to 
Goward and Wilkins (1955), Livesey and Smith (1953), and Hsiao and Telegdi 
(1953) ; 

(c) the preferred states of C’® and Be® are determined to some extent by 
isotopic-spin selection rules, in that above about 25 Mev. electric-dipole 
absorption leads to T = 1 states which then break up via T = 1 states of the 
intermediate nuclei in agreement with the predictions of Gell-Mann and 
Telegdi (1953). 


The oxygen reaction 
O' + hy = 4He* — 14.5 Mev. 


has been studied by the authors quoted above, also by Dawson and Bigham 
(1953) and by Millar and Cameron (1953). The method of identifying the four- 
pronged stars produced by this reaction is well-known, and the variation of 
yield with photon energy for 700 stars has been published by Goward and 
Wilkins (19520). Analysis of the modes of disintegration of the nuclei involved 
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has been carried out by Livesey and Smith (1953), by Goward and Wilkins 
(1952a), and by Millar and Cameron (1953). However, the behavior of the 
reaction above 30 Mev. has not been studied in detail, and no detailed analysis 
of a large number of stars has yet appeared. This paper describes the results 
of analyses carried out on 323 oxygen stars produced in Ilford C-2 plates ex- 
posed to bremsstrahlung of maximum energy 70 Mev. from the electron 
synchrotron at Queen’s University. In a later section the deductions are com- 
pared with those of other authors. 


2. EXPERIMENTAL PROCEDURE 


Ilford C-2 nuclear emulsions, 200 yw thick, coated on 1 in. by 3 in. glass 
backings were irradiated in the unfiltered 70 Mev. X-ray beam from the 
electron synchrotron at Queen’s University. All irradiations were monitored 
by Baldwin-Farmer BD-2 ionization chambers placed near the beam edge. 
These chambers were calibrated against a Victoreen thimble placed behind 
7.3 cm. of lucite at the position the plates later occupied. The plates were 
placed at a distance of 2 meters from the synchrotron target with the emulsion 
face perpendicular to the beam direction. Plate dosages varied from about 
10 to 50 roentgens. 

A standard temperature-change method was used for processing the plates. 
They were first immersed in a 2% sodium carbonate solution at room tempera- 
ture, then placed in a refrigerator kept at 3° C. and left there for several hours. 
Next they were transferred to a cold developing bath also at 3° C. consisting 
of 1 gm. hydroquinone, 20 gm. sodium sulphite, and 1 gm. potassium bromide 
in 2 liters of distilled water. (This developer has been found to give good 
discrimination between different types of particles.) The cold soak continued 
for at least 12 hr. and was followed by a development period of about 30 min. 
in the same type developer at 20° C. Development was arrested by immersing 
the plates in a 2% acetic acid solution for 30 min. Fixation started with a 
30 min. soak in fixer at 3° C. followed by a transfer of the plates to fixer at 
20° C. where they were left for about five hours. Kodak formula F-5 which 
contains a hardening solution was used for the fixer. A mechanical rocker was 
used for all processing carried out at 20° C. Following a five-hour wash in 
running water the plates were immersed in a 30% glycerin solution for about 
45 min. before being placed out to dry. 

Measurement of the shrinkage factor was carried out as follows. Immediately 
after an irradiation the thickness of the glass backing plus emulsion was 
measured at 10 different places on the plate using a bench gauge which read 
consistently to 0.0001 in. After processing but before the plates were put out 
to dry, the emulsion was removed from each corner of the plate. Thus it was 
possible to measure the glass thickness at these points with the bench gauge 
and to determine the original emulsion thickness at the intermediate points. 
The processed emulsion thickness was measured at the same intermediate 
points using the microscope fine-focus adjustment and the same set of ob- 
jective and eyepieces that would later be used for star measurement. Under 
reasonably favorable conditions this method allows the shrinkage factor to be 
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determined to within about 5%. The values of the shrinkage factor for the 
various plates used lay between 0.9 and 1.4. This method has the advantage 
that it calibrates the fine-focus adjustment of the microscope against the bench 
gauge. Day-to-day checks of the processed emulsion thickness were made in 
order to allow for variations occasioned by changing humidity and temperature. 

Measurement of the stars was carried out using a combination of objective 
and eyepieces which gave a field of view about 110 yu in diameter. The stars 
were identified using the momentum balance principle described by Goward 
and Wilkins (1950). 

The correction for stars not measured because any prong is not entirely 
contained in the emulsion was determined experimentally as follows. A plot 
of the number of stars versus star thickness (the distance between two planes 
parallel to the emulsion surface which just enclose the star) was made for each 
of the energy groups indicated in the yield curve. In group A, for example, 70% 
of the stars have thicknesses less than 10 yu while all the stars in this energy 
group have thicknesses less than 20 yu. Hence, if all the stars in group A were 
situated 2.5 uw (one-half the average thickness) from the emulsion surface about 
70% would be entirely contained in the emulsion. At a distance of 7.5 u from 
the surface approximately all group A stars would be contained in the emulsion. 
A curve of fraction of stars contained in the emulsion versus distance of the 
vertex from the emulsion surface was made for each energy group. The area 
under each curve out to one-half the emulsion thickness is divided by the 
area that would have been obtained had there been no loss, giving the fraction 
of stars entirely contained in the emulsion. The inverse of this quantity is 
then the required correction. Values of the correction factor varied from 1.02 
for group A to 1.12 for group E, in the yield curve (Fig. 10). 


3. RESULTS AND DISCUSSION 


The first step in the analysis of the oxygen stars consisted of the calculation 
of the balance of momentum in each of the three dimensions as already 
mentioned, and a star was accepted only if the net error in momentum balance 
(4) was less than 50 photon units. In these units the energy of a photon in Mev. 
is numerically equal to its momentum. The energy of the photon responsible 
for each star was found under the assumption that none of the energy was lost 
in excitation of the product nuclei. The number of stars per 0.5 Mev. interval 
was then plotted against photon energy, as is shown in Fig. la, and this yield 
curve exhibits several peaks similar to those reported by Livesey and Smith 
(1953) and by Goward and Wilkins (19525). The estimated error on each energy 
determination is of the order of +0.5 Mev., and this severely limits the reso- 
lution of the different energy groups. In order to improve the resolution slightly 
and to select a group of stars suitable for detailed analysis, a second yield curve 
is plotted in Fig. 1b with data from stars in which the value of A did not exceed 
35 photon units. This histogram is based on figures from 263 stars and it shows 
well-defined peaks close to 22.5 Mev. and 32.5 Mev. with possible groups at 
29 Mev., 35 Mev., and between 26 and 28 Mev. For the purpose of further 
analysis the stars were divided into six groups, marked on the diagram, but, 
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All stars with A ¢ 50 


No. of stars 


stors with A < 35 


millibarns 





Approximate cross - section in 


16 20 25 30 35 40 45 
Photon energy in MEV 


Fic. 1. (a) Yield curve for 323 stars with A< 50 photon units. (6) Yield curve for 263 
stars with A < 35 photon units, used in detailed analysis, with stars involving Be’ in the ground 
state shown shaded. (c) Estimated cross section curve for the reaction after corrections for 
escape of particles and the bremsstrahlung spectrum. 


owing to the limited resolution of the method, these groups are not sharply 


defined. 

The actual excitation function for the reaction is derived from the yield 
curve by applying corrections for the escape of tracks from the emulsion 
boundaries and for the shape of the spectrum of bremsstrahlung emitted from 
the synchrotron target. The Schiff spectrum was used in this work (see Trainor 
and Brown 1955), and the dosage in roentgens as measured by the Victoreen 
chamber was used to find the integrated cross sections for the various groups. 
In Table I the energy limits of the groups and the corresponding cross sections 
are shown. The total integrated cross section up to 35.5 Mev. is approxi- 
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TABLE I 
ENERGY GROUPS AND INTEGRATED CROSS SECTIONS 


Energy range Integrated cross section 
Group in Mev. in Mev-millibarns 
A Up to 25.0 1.0 
B 25.0-26.5 0.3 
¢ 26.5-28.0 0.2 
D 28.0-30.5 0.5 
E 30.5-33.5 0.6 
F 33.5-35.5 0.2 


mately 2.8 Mev-mb. Our result of 1.5 Mev-mb. is in agreement with the 
estimate of 1.4 Mev-mb. for the energy range up to 28 Mev. given by the data 
of Goward and Wilkins (19520). 

The variation of the cross section with photon energy is shown in Fig. Ic 
and, although there are considerable errors both in the energy resolution and 
in the absolute value of the cross section, the function shows marked resonance 
effects with the maximum cross section occurring near 32.5 Mev. At low 
energies the dotted line allows for losses of low-energy stars with tracks too 
short for measurement. The position of group A agrees fairly closely with the 
peak in the “giant resonance’’ curve shown by the (y, ”) and (vy, p) reactions 
in oxygen, but, as will be shown later, this group is probably complex in 
structure. Group B agrees in position with a peak in the excitation function 
published by Goward and Wilkins (19526), but the peaks corresponding to 
groups D and E have not previously been detected. It is interesting to note 
that the cross section tends to increase over the energy range between 25 and 
30 Mev., where the cross sections of the (y, 7) and (y, p) reactions are believed 
to decrease rapidly. This effect may be due to some competition between the 
various reactions, as has been discussed by Goward and Wilkins (1953) in 
connection with the photodisintegration of C™ nuclei. 

A notable feature of many stars, especially those belonging to group A, is 
the occurrence of two alpha-particle tracks of similar length very close to- 
gether. This happens much more often than a random orientation of the tracks 
would allow, and it is clear that each pair of tracks arises from the intermediate 


nucleus Be’ in its ground state. After allowing for a small number of stars of 
this kind as arising fron chance juxtaposition of tracks, 49 stars were allotted 
to the ‘‘ground-state”’ category, and these are shown heavily shaded in Fig. 16 
and also in Fig. 4. This yield curve exhibits a single well-marked peak at 


22.3+0.2 Mev., but a small proportion of the stars at all the energies investi- 
gated belong to this category. 

In any star showing a pair of tracks derived from Be’ in the ground state it 
is possible to calculate two values of the excitation energy of C” as an inter- 
mediate nucleus. Of these only one in each case can be significant, so the histo- 
gram of all calculated C” excitation energies, shown in Fig. 2, should consist 
of a continuum of non-significant results with superimposed peaks where any 
level of C!? makes an appreciable contribution to the process. The only clear 
evidence of a peak is at 10.8+0.2 Mev., but examination of the individual 
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_ Fic.2. Histogram of calculated excitation energies of C” as intermediate nuclei in 49 stars 
involving Be® in the ground state. One-half of the values only may be significant. 


stars shows that a single level cannot explain all the stars observed, and the 
known level at 9.6 Mev. in C!? (see Ajzenberg and Lauritsen 1955) must be 
invoked to account for about 30% of them. The present evidence for levels at 
9.6 Mev. and 10.8 Mev. agrees with the results of Goward and Wilkins (1952a) 
and those of Livesey and Smith (1953). 

The stars not involving the ground state of Be* were analyzed to determine 
if any levels of C!® were operative, but here each star yielded four values of the 
excitation energy, of which only one can be significant. In general it was found 


c™ excitation energies ip stars not involving 
the ground state of Be 


group A 





10 18 20 26 
c'® excitation in MEV 


Fic. 3. Histograms of calculated excitation energies of C' in stars not involving Be® in the 
ground state: (a) for group A stars; (b) for groups B+C; (c) for groups D+E+F. One-fourth 


f the values may be significant. 
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that the groups B and C behaved similarly, likewise the groups D, E, and F, 
so in Fig. 3 we show three histograms of C™ excitation values, for groups A, 
B+C,‘and D+E-+F respectively. In the case of group A there is a possible 
peak at 12.3+0.3 Mev., but this does not rise significantly above the dotted 
curve estimated for a random orientation of tracks. In groups B+C, however, 
there is clear evidence of a peak between 15 and 16 Mev., and it is probable 
that here the well-known T = 1 levels of C” at 15.1 and 16.1 Mev. are in- 
volved. This is in agreement with the theory of Gell-Mann and Telegdi (1953) 
and with previous experimental results by Hsiao and Telegdi (1953). The 
D+E-+F histogram also shows evidence of the level at 16 Mev., and a second 
peak between 18 and 19 Mev. This latter peak is probably due to the levels at 
18.3 and 18.9 Mev., which are known to lead to the emission of three alpha- 
particles (Goward and Wilkins 1953). 

The Be® excitation energies for different pairs of tracks in all the stars in- 
cluded in Fig. 16 have been worked out, and the appropriate histograms are 
shown in Fig. 4. Here the significant values number only one-sixth of the total 


Be’ excitation energies in stars with A € 35 


groups B+C 


group A 


No. of stars 


groups D+E +F 





8 oO 5 10 15 
Be excitation in MEV 


Fic. 4. Histograms of calculated Be* excitation values in stars belonging to different energy 
groups: (a) for group A stars; (6) for groups B+C;: (c) for groups D+E+F. One-sixth or 
one-third of the values may be significant according to the mechanism postulated. 


if the reaction proceeds by successive emission of single alpha-particles, and 
one-third of the total if two Be* nuclei are formed from the oxygen nucleus, 
but no clear evidence of the latter reaction has been obtained and the foregoing 
evidence on C® levels seems to rule out the possibility except in a small pro- 
portion of cases. The diagram for group A stars shows a well-marked peak for 
the ground state stars, and it is of interest to note that the average ground 
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state energy relative to two free alpha-particles calculated from these data is 
0.09+0.02 Mev., in agreement with more precise determinations quoted by 
Van Patter and Whaling (1954). There is no evidence of the 3 Mev. level of 
Be® occurring in this group, but it is almost impossible to detect this level in 
the presence of a large number of non-significant values distributed about a 
mean close to 2.9 Mev. The behavior of groups B+C is markedly different 
from that of group A, and it seems to be quite clear that the 3 Mev. level does 
not occur in these stars. The histogram shows a peak near 5 Mev. which is of 
particular interest in the light of recent evidence for Be® levels in this region. 
Goward and Wilkins have reported the presence of a 4.1+0.1 Mev. Be® level 
in C(y, He*)Be®* stars, and Gibson and Prowse (1955) have published evi- 
dence of Be’ levels at 4.1 and 5.3 Mev. in the Li’(d, n) Be®* reaction. However, 
the present data do not give definite information on the point and they do not 
confirm the suggestion of Livesey and Smith that the supposed level at 4.1 Mev. 
is responsible for group B stars. Again, in groups D+E+F there is no clear 
evidence of Be® levels except a peak near 3 Mev. which seems to be significant 
when compared with the trend of the expected random curve in its vicinity. 

These results are in good agreement with previous deductions concerning 
this reaction. Above 28 Mev. the dominant mode of disintegration is via levels 
in C” at 16.1, 18.3, and 18.9 Mev. leading to the broad 3 Mev. level in Be’. 
It is believed that the level at 16.1 Mev. in C” is of the type (2+, T = 1) 
(Inglis 1953), and this would readily yield that Be® level with the emission of 
an s-wave alpha-particle, but one cannot rule out the possibility that some of 
the levels quoted are of the type 1-. The stars in group A show evidence of 
at least three different modes of disintegration and no clear indication of the 
spin relations has been obtained. The remaining groups B and C show some 
trace of the T = 1 level in C at 15.1 Mev., and this is believed to be of the 
type 1+, in which case transitions to the 2+ level in Be® at 3 Mev. require the 
emission of a d-wave alpha-particle. It is not surprising, therefore, that this 
mode does not occur very often, but, on the other hand, no other even-parity 
levels in Be® (which alone can give two alpha-particles) are favored. Until the 
difficult question of the existence of Be® levels in the 4 to 5 Mev. region is 
settled, one cannot draw definite conclusions concerning the stars in groups 
B and C. 

4. CONCLUSIONS 


A great deal of this work has been in the nature of a confirmation of previous 
results with rather more accurate data. We summarize here the salient points 
emerging from the analysis: 

(i) the excitation function of the reaction shows definite peaks at 22.5 and 
32.5 Mev. and possible peaks at 26.0, 27.5, 29.5, and 35.2 Mev.; 

(ii) a small proportion of stars at all energies proceed via the ground state 
of Be’, and the yield of this type of star shows a sharp peak at 22.3 Mev.; 

(iii) the stars in group A (below 25 Mev.) may disintegrate via levels in 
C” at 9.6 and 10.8 Mev., leading to the ground state of Be*, or by other 
mechanisms not ascertained; 
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(iv) the stars in groups B and C (25 to 28 Mev.) disintegrate to some extent 
via the T = 1 levels in C” at 15.1 and 16.1 Mev.; 

(v) the high energy stars (28 to 35 Mev.) proceed chiefly via C” levels at 
16.1, 18.3, and 18.9 Mev. and the broad 3 Mev. level in Be’. 
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THE FAR ULTRAVIOLET (C—X) BANDS OF N,+! 


By P. G. WILKINSON 


ABSTRACT 


Emission bands of the N2+ C—X system in the vacuum ultraviolet have been 
obtained under high resolution. The rotational analysis confirms the assignment 
as a *?2{ —?*2f transition and the rotational constants obtained are in sub- 
stantial agreement with previous investigations. The unusual intensity distribu- 
tion in the v’ = 3 progression and perturbations in the v’ = 0, 1, and 2 levels of 
the . — are in agreement with the theory of inverse predissociation in the 
v’ = 8 level. 


A. INTRODUCTION 


Three band systems of N2+ have been reported, the C?2+ — X*Z* system 
in the vacuum ultraviolet, the B*?2+ — X?2Z+ system in the near ultraviolet, 
and the A?II, — X?*Z+ system, first discovered in auroral emissions in the near 
infrared. The rotational structure of the ultraviolet bands has been studied 
by Childs (1932), Coster and Brons (1932), Crawford and Tsai (1935), Parker 
(1933), and Douglas (1952). Douglas showed that the B?Z> state has a dis- 
sociation energy at 70358 cm~', which value gives strong support to the 
9.756 ev. dissociation energy of nitrogen. The auroral bands were obtained 
under low resolution by Meinel (1951), and soon after duplicated in the 
laboratory (Dalby and Douglas 1951; Douglas 1953). Watson and Koontz 
(1934) made a vibrational analysis of the C—X system which was corrected 
by Setlow (1948) and Miescher and Baer (1952). Takamine et al. (1939) and 
Tanaka (1953) have made intensity studies of these bands excited under 
different conditions. The present paper reports a high resolution study of a 
number of the C—X bands, and several intensity anomalies are discussed. 


B. EXPERIMENTAL 


The C—X bands were excited in an iron hollow cathode discharge tube 
operating at 0.2 to 0.6 amp. in a mixture of flowing tank helium and nitrogen 
(1%) at a pressure of 10 mm. The photographs were taken on a 21-ft. normal 
incidence vacuum spectrograph (Wilkinson and Mulliken 1955) employing a 
Bausch and Lomb concave diffraction grating (30,000 lines per inch, 150 mm. 
by 45 mm.). The plate dispersion was 1.33 A per mm. in the first order, and 
under favorable conditions a resolving power of 90,000 was obtained in 
the 1500A-2000A region using Ilford Q2 photographic plates. The exposure 
times were about one hour. The plates were reduced in the usual way employing 
Fel and Fel] lines as wavelength standards (Green 1939). 


C. RESULTS 
An enlargement of the 4-10 and 3-9 bands is shown in Fig. 1. The bands 


possess only P and R branches, which are resolved almost to the head, and the 
expected intensity alternation (even N strong, odd N weak) is exhibited. The 


1Manuscript received October 26, 1955. 
Contribution from the Laboratory of Molecular Structure and Spectra, Department of 
Physics, The University of Chicago, Chicago, 37, Illinois, U.S.A. 
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classification of the transition as ?2*+ — *Z+ (Setlow 1948) is confirmed but no 
spin fine structure doubling is observed under the resolving power employed. 
The bands possess exactly the same structure as a '>—'Z transition except 
that the lines are numbered JN instead of J (Hund’s Case 6). The wave numbers 
of the 3-6, 3-8, 3-9, 3-10, 4-10, and 4-11 bands are given in Table I. The data 


TABLE I 
WAVE NUMBERS OF THE C—X BANDS 


3-6 3-8 3-9 3-10 
N R(N) P(N) R(N) P(N) R(N) P(N) R(N) P(N) 


4 52256.32* 52224.53* 50384.19* 50357.64* 
5  58079.02* 58047.00 54162.88* 54125.06*  56.32* 18.64" 84.19% 52.05 
é 5 23 

















5 79 .02* 40.60 62. 88* 3.83 56.32* 16.62 84.19* 46.03 
7 79 .02* 33.52 62.88* 18.59* 56.32* 12.10* 84.19* 39.59 
8 79 .02* 24.36 62. 88* 09. 86 56.32* 05.12* 84.19* 33.21 
9 76.06 17.83 58.92 03 .65* 53.92 52197.76* 84.19* 27 .94* 
10 69.41 06 . 86 56.60 54095.40* 52.06 89 .97* 81.46 18.72 
1l 65.33 57999.26 53.81 84.98 49.61 82 .09* 79.28 10.62 
12 61.29 86.41 50.38 75.45 46.71 71.89 76.89 02.19 
13 53.20 75.75 46.40 65.39 43.19 62.47 74.16 50293.39 
14 50.67 63.92 41.82 55.14 39.10 52.57 70.73 83.95 
15 44.64 48.27 36.66 43.95 34.79 42.13 66.95 74.40 
16 37.77 38.94 31.21 32.49 29.82 31.28 63.07 63.94 
17 30.26 24.34 25 .06* 20.36 24.53* 19.95 57 .64* 53.30 
18 22.21 11.59 18.59* 07.78 18.64* 08.38 52.97 42.09 
19 13.66 57896.77 11.38  53994.69 12.10* 52095.73 47.40 30.58 
20 04.10 81.56 03 . 65* 81.34 05.12* 83.09 41.31 18.75 
21 57994. 21 65.45  54095.40* 66.83  52197.76* 69.64 35.10 06.25 
22 83.72 49.21 86.74 52.30 89.97* 55.92 27.94* 50193.39 
23 72.19 32.28 77.61 37.14 82.02* 41.64 20.69 80.10 
24 60.83 14.63 67.64 21.41 73.25 26.94 13.28 66.35 
25 57796.52 57948.27 57.74 05.13 64.11 11.79 04.64 52.19 
26 77.42 35.55 46.54 53888.28 54.26 51996.10  50295.99 37.48 
27 57.93 22.16 34.13 69.70 43.98 79.91 86.28 22.75 
28 38.05 07.97 23.39 53.28 33.38 63.28 77 .33 07.13 
29 17.16 57892.88 11.11 35.30 22.21 46.35 
30 57696 .15 78.07 53998.22 16.30 10.56 28.76 50075. 16 
31 74.98 84.72 53796.90 52098.34 10.81 
32 51.91 70.97 77.29 85.91 51892.36 
33 58.65 57.05 72.97 73.52 
34 41.36 36.11 59.44 54.20 
35 15.07 44.92 33.94 
36 10.39 53693.03 30.92 13.95 
4-10 4-11 
N R(N) P(N) R(N) P(N) 
1 50518.27* 50509 .34* 
2 52341 .52* 20.41* 05 .43* 
3 36 .95* 22.15* 01.12* 
4 52359 .09* 31.96 24.49* 50496.48* 
5 60.57* 26 .57* 24.49* 91.33 
6 60.57* 20.63 24 .49* 85.75 
7 60. 57* 16.14 24.49* 79.76 
8 60. 57* 09 .52* 24.49* 73.44 
9 59 .09* 00.16 24.49* 66.66 
10 57.28 52292.45* 22.15* 59.36 
11 52.29 86 .35* 20.41* 51.78 
12 49.17 77 .90* 18.27* 43.56 
13 45.55 66.32 15.66 35. 26 
14 41.52* 56.32 12.58 26.36 
15 36 .95* 46.71 09 .34* 16.65 
16 33.52 34.79 05 .43* 07.14 
17 26. 57* 24.53 01.12* 50396.88 
18 22.54 12.10 50496.48* 86.18 
19 14.28 52197.76 
20 09 .52* 89.97 
21 02.35 73.25 
22 52292 .45* 62.47 
23 86 .35* 
24 77 .90* 31.28 


*Lines uncertain because of overlapping. 


represent an improvement over previous work in that P and R branches are 
resolved, making possible independent determinations of the rotational 
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constants in the upper and lower states. A summary of the rotational constants 
found, along with the band heads and origins, is given in Table II. The con- 


TABLE II 


ROTATIONAL CONSTANTS AND BAND POSITIONS OBTAINED FROM ANALYSIS OF THE C—X BANDS 
(BH—band head; »p—band origin) 


v’ vy” B’ D’.108 2 D’’.108 vo 
(em) (cm)? = (cm”) = (cm™) (cm™) (cm=) 

3 6 1.499 4.5 1.805: 6.16 58070. 30 58079 .0-2 
3 8 1.4995 4.6 1.7654 6.25 5§4152.8 54162.8, 
3 9 1.499» 4.6 1.7456 6.20 52246. 25 52256. 32 
3 10 1.499, 4.6 1.724, 6.2 50373 .0¢ 50384. 1, 
4 10 1.478 4.5 1.724, 6.2 52350.1,  52360.57 
4 11 50512.9 50524.8 
4 9 54222.4* 54232 .0 
5 10 54294. 8* 54305 .1 
5 11 52459 .5* 52470 .6 


*Band origin computed from band head position and rotational constants. 


stants B and D were computed from the combination relations, A,F’(N) = 
R(N)—P(N) and A:F’(N) = R(N—-1)—P(N+1) (Herzberg 1950). The 
rotational constants obtained for the X?Z+ ground state are consistent with 
the data of Crawford and Tsai (1935) and Douglas (1952). Since the usual 
B, versus v plot is nonlinear in this case, it is necessary to introduce a (v+ 4)? 
term in order to fit the data properly. The resulting equation, B, = 1.9253— 
0.0174;3(v-+ 4)-0.00016,(v+4)?, fits the present data within +0.0005 cm™ 
and previous work within +0.002 cm—!. The value of B, (1.926 cm-") is smaller 
than previously given (Douglas 1952), and the corresponding value of r, 
(1.118 A) is somewhat higher. The combined data of Crawford and Tsai 
(1935), Douglas (1952), and the present work are given in Fig. 2. The only 
rotational constants for the upper state previously available are those of 
Setlow (1948): B;’ = 1.493 and B,’ = 1.480 cm, which compare with 1.499 
and 1.478 obtained in the present work. 
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Fic. 2. B, graph for the on electronic state of N3*+. The value 0.0160(0-+-4)-0.000164 
X (vw +4)? has been added to order to allow a large-scale linear plot. 
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D. DISCUSSION OF RESULTS 


Vibrational perturbations are observed in the v’ = 0, 1, and 2 levels of the 
C?=+ state as may be seen in Fig. 3. Tanaka's (1953) band head data were used 





Fic. 3. AG curve for the C?2{ state of N2*. The curve was constructed from the data of 
Tanaka (1953). None of the points should be in error by more than a few cm™!. 


to construct the AG,,; versus v plot in the figure. Although the use of band 
heads undoubtedly introduces an error, this is only about 2.5 cm~',? whereas 
the AG,4 curve in Fig. 3 deviates 33 cm~! at v’ = 0, 16 cm™ at v’ = 1, and 
7 cm~' atv’ = 2 from a linear relationship. Douglas (1952) has considered that 
the unusual shape of the AG,,; versus v curve in the B?* state is caused by 
interaction with the first three vibrational levels of C*2+ and predicts that 
Bo, By, and Bz will be anomalously large in C?Z+. Setlow (1948) reported 
By = 1.62 cm~, which is unusually large, but this value must be considered 
uncertain since it was obtained by analysis of several only partially resolved 
bands (0-5, 0-6). 

There are a number of peculiarities in the observed intensity distribution 
of the C—X bands which depend on the excitation conditions. Transitions 
from v’ = 3 and 4 are by far the strongest especially in helium-nitrogen 
mixtures. This was explained by Watson and Koontz (1934) as resonance 
excitation by Het ?S since the ionization potential of helium (24.580 ev.) 
nearly coincides with the energy of the v’ = 3 level of the C? Lf state (24.330 
ev. above v = 0 of No X'2F). The intensity distribution in the various bands 
as given by Tanaka (1953) is shown in Table III (top) and as found under 
the present conditions (bottom). In this work the excitation of the v’ = 3 and 
4 progressions was by far the strongest and the 3-9 band was the strongest 
in the system; in addition, these progressions appeared only in the pressure 
range 10-20 mm. Above 20 mm. the bands were very weak or absent. If the 
mechanism involves resonance excitation by Het (24.580 ev.) then the 


2This error was estimated using reasonable rotational constants obtained from Setlow (1948), 
Douglas (1952), and this research. 
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TABLE III 


INTENSITIES IN THE C—X BANDS. DATA GIVEN BY TANAKA? ARE AT THE TOP AND DATA FOUND 
IN THE PRESENT WORK ARE AT THE BOTTOM? 














ae 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
0 0 0 
1 0 0 0 
2 0 0 9 
3 3 3 2 5 4 3 6 6 6 7 3 
4 2 3 2 5 5 0 3 5 6 2 
5 0 0 4 5 1 
6 0 3 3 
© 0 1 0 
8 0 
3 4 6 10 6 
4 1 6 5 
5 0 0 





*Y, Tanaka (1953). The 0-5, 0-6, 1-6, 1-7, 1-8, 2-7, 2-8, and 2-9 bands were obtained in 
disruptive discharge of nitrogen and nitric oxide. All other bands were obtained in nitrogen 


and helium mixtures. 
’Intensities obtained in hollow cathode discharge through He - 1% Nz. 


v’ = 4 level (24.575 ev.) should be preferentially excited rather than the 
v’ = 3 level (24.330 ev.) (Tanaka 1953). Douglas has suggested (1952) that 
the excitation mechanism may be a case of inverse predissociation. This is 
supported by the observations that the dissociation energy of the B?Z+ state 
is 70358 cm~! (above v = 0 of N.+ X?=*),* which is 269 cm~' below the 
energy of the v = 3 level of C?2*, and that under certain discharge conditions 
(see above) band progressions with v’ = 3 and 4 do not appear. The existence 
of vibrational perturbations in the v’ = 0, 1, and 2 levels as well as the appear- 
ance of strong NI and NII lines in the spectra also lends support to this 
theory (Tanaka 1953). It is reasonable to suppose that both processes (reso- 
nance excitation by Het+ and inverse predissociation) may be operating in 
tubes containing helium. Predissociation in the v’ = 3 level is consistent with 
either the 7.373 or the 9.756 ev. dissociation energy of nitrogen although a 
potential hill of 0.029 ev. must be involved in either process. 

Although the perturbations in the C?2+ state are undoubtedly real, and it 
seems possible there could be interaction with B?2* (Douglas 1952), it appears 
that the perturbations and predissociation in C?2* could also be understood 
from an interaction with a slightly attractive state with its asymptote at 
approximately 70600 cm~! (8.75 ev.). If the 9.756 ev. dissociation energy of 
nitrogen is assumed (corresponding to 8.721 ev. for the dissociation energy 
of N.*), then the products of the predissociation of N,*+ at 8.75 ev. must be 
N*+(@P)+ N(4S). Consideration of selection rules gives the result that possible 
predissociating states are 72+, 42+, ®2*, 711, ‘11, and ®Il,. 

If the Franck—Condon principle is employed to predict relative band 
intensities (Setlow 1948) one can see immediately that 3 — 9 is the most 
probable transition and that 3 — 8 and 3 — 10 should be weaker and about 
equal in intensity in agreement with experiment (Table III). One would also 


’\ small discrepancy exists between the experimental dissociation energy of Ne* as derived 
from the B —X bands (Douglas 1952) and that derived from the dissociation energy of nitrogen 
(Douglas and Herzberg 1951) and the ionization potentials of the nitrogen atom and the 
nitrogen molecule. The latter gives 70396 cm~! and must be considered the more reliable. 
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predict the 3 — 7 and 3 — 6 transitions with decreasing intensities. In actual 
fact the 3-7 band was not observed in the present work and the 3 > 6 
band appears nearly as strong as the 3 — 8.4 


I wish to thank Professor R. S. Mulliken of the University of Chicago and 
Dr. Y. Tanaka of the Air Force Cambridge Research Center, for many helpful 
discussions. This research was assisted in part by the Geophysics Research 
Directorate of the Air Force Cambridge Research Center, Air Research and 
Development Command, under Contract No. AF19(604)-1019 with the 


University of Chicago. 
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HIGH-ENERGY FORBIDDEN 8-RAY TRANSITIONS FROM Cs™ 
(2.3 yr.), Co® (5.3 yr.), Sc** (84 d.), and Hg?® (47.9 d.)! 


By J. L. WoLFson? 


ABSTRACT 


Beta-ray transitions of energies 1478+6 kev. from Co® (5.3 yr.) and 1475+6 
kev. from Sc** (84 d.) were observed with intensities of (1.0+0.2) X10~ and 
(3.6+0.7) X10~* f-rays per disintegration respectively. Forbidden spectrum 
shapes were obtained for both transitions. No evidence was obtained for the 
8-ray transition of energy 1451 kev. from Cs™ (2.3 yr.), nor for the transition 
of energy 473 kev. ion Hg™ (47.9 d.), and upper limits of 5X10-* and 3X10~ 
8-rays per disintegration respectively are placed upon the intensities of these 
transitions. The log fot values are S14. 5, f4. 0+0.1, 13.0+0.1, and >11.3 for 
the transitions referred to from Cs™, Co, Sct, and Hg? respectively. 


1, INTRODUCTION 


The magnitude of the matrix elements operative in beta decay is expressed 
in terms of the ft value (or log ft), the quantity ft being proportional to the 
inverse square of the nuclear matrix element. The ft values are known for a 
large number of 8-ray transitions (Goeppert-Mayer 1955), but the vast 
majority are for allowed and first-forbidden transitions. In particular, only a 
relatively small number are known for forbidden transitions of the second 
order, partly because of the scarcity of simple spectra of high orders of forbid- 
denness, and partly because, among the more complex modes of disintegration, 
the transitions of higher orders of forbiddenness must compete with the much 
more probable allowed or first-forbidden transitions, and their experimental 
observation becomes difficult. Nevertheless, since these forbidden transitions 
have higher energies than those transitions of lower orders of forbiddenness 
with which they compete, it is possible to observe at least the higher energy 
portion of their spectra, and thence deduce their intensities. In the present 
work this has been attempted with the §-ray transitions proceeding from 
the decay of Cs!*4 (2.3 yr.), Co® (5.3 yr.), Sc*® (84 d.), and Hg? (47.9 d.). 
The transitions sought from the decay of the first three nuclides are expected 
to be second forbidden, whereas that from the decay of Hg? is expected to 
be first forbidden. 

In all of the transitions studied it was of importance to correct for or to 
avoid the detection of conversion or secondary electrons from the relatively 
high energy and high intensity gamma-ray transitions in the daughter nuclei. 
The second of these alternatives was chosen by restricting the region studied 
to that portion of the 8-ray spectrum of energy greater than that of the most 
energetic y-ray emitted from the source. While this procedure rendered the 
detection of the weak 8-ray transitions rather difficult, it did possess the virtue 
of yielding unambiguous results. 


1Manuscript received November 4, 1955. 

Contribution from Physics Division, Atomic Energy of Canada Limited, Chalk River 
Laboratory, Chalk River, Ontario. Issued as A.E.C.L. No 

*Present address: Jewish General Hospital, Montreal, Quebec. 
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A report of the results obtained in these studies was presented at a meeting 
of the Canadian Association of Physicists at London, Ontario, in May, 
1953. 


2. APPARATUS AND METHOD 


The method employed in each of these studies, as indicated in Section 1, 
involved a detailed examination of the spectrum sought at energies in excess 
of the most energetic y-ray emitted by the daughter nucleus. A double-lens 
magnetic spectrometer was used, adjusted for a line width of about 0.9% in 
momentum at half-maximum intensity. The transmission of the spectrometer 
was about 0.4% of 4x steradians. The low level of the background counting 
rate achieved by the spectrometer, as discussed elsewhere (Elliott e¢ a/. 1954), 
constitutes a most valuable feature of the instrument. The background was 
further reduced in these investigations by employing pulse height selection, 
using a single-channel pulse height analyzer. In the studies with Cs!*4, Co®, 
and Sc‘**, a Nal(TI) crystal, sufficiently thick to stop all electrons whose 
energies lay in the range of interest, was used as detector. For examination 
of the spectrum from decay of Hg’, the Nal crystal was replaced by an 
anthracene crystal 1 mm. thick. 

The sources were all prepared by neutron irradiation in the NRX reactor. 
Since only those regions of the spectra near their end-points were examined, 
the thickness of the sources was not as important a consideration as it usually 
is in studies of spectrum shapes. Nevertheless, considerable effort was made 
to obtain the highest possible specific activities, as well as uniform deposition 
of the source material. In general the sources were prepared by evaporation 
of salts from nitric acid solution on backings of 2.3 mgm./cm.? aluminum 
foil. The sources were all approximately 2 mm. in diameter. 

The intensity of a transition whose spectrum was detected was found by 
extrapolation of the spectrum to zero energy, followed by comparison of the 
area underneath the spectrum with that of the K internal conversion line of 
the most energetic y-ray transition. The intensity of the 8-ray transition was 
then deduced from the known A conversion coefficient and the quantum 
intensity. This method has the advantage of minimizing the effects of source 
thickness, since all electrons observed are of high energy. 

3. RESULTS 
Ce** (2:9 9.) 

The source was about 9 mgm./cm.” thick.and about 0.75 me. in activity. 
The results of an examination of the spectrum above 1300 kev. are shown in 
Fig. 1. The arrow in the right-hand figure at about 6320 gauss-cm. indicates 
the end-point expected for the 1451 kev. beta spectrum. There is no evidence 
for this spectrum. If the disintegration scheme of Keister, Lee, and Schmidt 
(1955) is correct, an upper limit of 5X10-° 6-rays per disintegration can be 
placed on the intensity of this transition. The corresponding log fot value is 
> 14.5. These figures are based on Keister, Lee, and Schmidt’s (1955) values 
for the intensity of the 1367 kev. y-ray transition and its K conversion co- 
efficient. 
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Fic. 1. The po spectrum of Cs™ (2.3 yr.) above 1300 kev. The K, L, and M internal 
conversion lines of the 1367 kev. y-ray transition in Ba are shown on the left. The counting 
rates at energies greater than 1367 kev. are shown on the right. The arrow indicates the end- 


point of the expected 1451 kev. 8-ray spectrum. There is no evidence for the existence of this 
spectrum. The vertical bars through the experimental points indicate the standard deviations. 


Co® (5.3 yr.) 

The source was about 9 mgm./cm.? thick and approximately 1 mc. in 
activity. The results of an examination of the beta spectrum above 1270 kev. 
are shown in Fig. 2. A forbidden Fermi plot, using the shape factor C27 
(Greuling 1942) corresponding to a spin change of three without parity change 
(Section 4), is shown in the inset. The extrapolated end-point of the plot is 
1478+6 kev. in satisfactory agreement with the expected end-point of 1481+3 
kev. Further, the straight-line shape of the plot is consistent with a spin 
change of three without parity change for the transition.* The intensity of the 
transition is (1.0+0.2) X 10-4 B-rays per disintegration. This value is based on 
a total conversion coefficient of 1.32 10-4 for the 1332 kev. y-ray transition, 
as required for 2 radiation. Using the curves of Moszkowski (1951) and 
Davidson (1951), the corresponding log fot and log fot values are 14.0+0.1 
and 13.6+0.1 respectively. 

Sc*® (84 d.) 

The source was about 6 mgm./cm.? thick and about 2.7 me. in activity. 

The results of the examination of the beta spectrum above 1050 kev. are 


shown in Fig. 3. An allowed Fermi plot of the data is shown in the inset. 
The plot is somewhat concave toward the energy axis, and has an extrapolated 





*While the data found here are not sufficiently extensive to rule out the possibility of a 
Al = 2 transition, the evidence is overwhelmingly in favor of the AJ = 3 transition (Section 4), 
and the corresponding shape factor has been used in finding the intensity of the transition. 
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Fic. 2. The 1481 kev. 8-ray spectrum of Co® (5.3 yr.). The K internal conversion line of 
the 1332 kev. y-ray transition in Ni® is shown on the left. The Z and M internal conversion 
lines are not resolved from the K line. The inset shows a forbidden Fermi plot of the data 
obtained at energies greater than 1350 kev., using the shape factor Czr applicable to a AJ = 3 
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Fic. 3. The 1477 kev. -ray spectrum of Sc** (84 d.). The K internal conversion line of the 
1119 kev. y-ray transition in Ti** is shown on the left. The L and M internal conversion lines 
are not resolved from the K line. The inset shows an allowed Fermi plot of the data obtained 
at energies greater than 1150 kev. The vertical bars indicate the standard deviations. 
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end-point of 1475+6 kev. in agreement with the expected end-point of 1477+8 
kev. (provided that the 2+ state in Ti‘ is at 890 kev.). This finding will be 
discussed in Section 4. 

Since it is expected that the spectrum is characterized by a spin change of 
two without parity change (Section 4), the shape of the spectrum is not 
expected to be unique. Inasmuch as only a very small portion of the spectrum 
has been observed, no attempt has been made to derive the shape factor Cop 
which will fit the experimental points. The shape factor S} ([p2+0.69?], as 
for Cl8* (Wu 1955)) has been found to give a satisfactory fit,* but this is not of 
great significance since any factor involving the square of the momentum will 
have much the same effect for the restricted high-energy portion of the spec- 
trum observed. When the S} factor is applied for the purpose of extrapolating 
the spectrum to zero energy, an intensity of (3.6+0.7) X10-5 B-rays per 
disintegration is found for the 1477 kev. transition. This value is based on a 
total conversion coefficient of 9.86 10-5 for the 1119 kev. y-ray transition as 
required for £2 radiation. The corresponding log fot value is 13.0+0.1. 


Hg? (47.9 d.) 

In the attempt to detect a 473 kev. 8-ray transition from Hg?, Hg enriched 
to 95% in the isotope Hg**? was used in the irradiation in order to obtain the 
highest possible specific activity. The low energy of the 6-rays, as compared 
with those from Co® and Sc**, makes the source thickness of greater impor- 
tance in this case. The source was about 6 mgm./cm.? thick and approximately 
0.5 mc. in activity. 

The results of the study are given in Fig. 4. Two runs taken 113 days apart 
are shown. It is immediately evident that the spectrum above 310 kev. decays 
at a rate much slower than that corresponding to the A internal conversion 
line of the 278.3 kev. y-ray transition. Eight very weak K internal conversion 
lines (not shown in the figure) have been observed in the spectrum. These 
have been identified with those resulting from the decay of 270 d. Ag!!®. 
Corresponding to two of these K lines, two ZL conversion lines have been 
observed, the A—L difference in each case being in excellent agreement 
with the A —L difference in Cd. Two @-ray transitions of energies 520 kev. 
and 2.9 Mev. occur in the decay of 270 d. Ag!'®, of which the former is intense 
and the latter weak. The 8-ray spectrum shown in Fig. 4 is consistent with the 
assumption that it results from the decay of 270 d. Ag'!®. The major part of 
the spectrum ends at about 520 kev. (shown by the arrow in Fig. 4), while the 
portion above 520 kev. can be accounted fer by the weak 2.9 Mev. component. 
The curious hump in the region of 520 kev. is probably caused by the conver 
sion lines in that vicinity, 

The decay of the spectrum is also consistent with a half-life of 270 days. 


*S>— p*+aq? for the five known twice-forbidden non-unique spectra, where \ = 0.6 for 
the Cl spectrum, —2 for the Tc spectrum, and > 10 for the Cs'®, Cs!87, and I spectra 
(Wu 1955). A good experimental value for \ can clearly only be found if the low-energy portion 
of the spectrum is available for study. The arbitrary choice of ’ = 0.6 for the 1475 kev. 
spectrum is justified only by the necessity of making some choice, and Sc" is similar to C6 
in that both are low Z odd-odd nuclei. If 4 > 0.6 the intensity of the 1475 kev. transition will 
have been underestimated, by about 20% if A = 1, and almost a factor of 2if X = 2. 











WOLFSON: FORBIDDEN TRANSITIONS 261 







K 
278.3 kev { 






Hg” @ SPECTRUM ABOVE 31/0 kev 


x- TAKEN 113 DAYS AFTER POINTS MARKED BY o 


COUNTS PER MINUTE 
nn 


473 kev, 520 kev. 


10000 


Bp (GAUSS-CM) 


Fic. 4. The §-ray spectrum of Hg** (47.9 d.) above 310 kev. Two runs taken 113 days 
apart are shown. The K internal conversion line of the 278.3 kev. y-ray transition in TI 
is shown on the left. The 8-ray spectrum is shown on the right. The expected end-point of the 
473 kev. B-ray spectrum of Hg** is marked by an arrow. The spectrum observed is due to 
270 d. Ag"® which occurs as a contaminant in the source. The arrow at 520 kev. indicates the 
end-point of the 520 kev. 8-ray component of Ag"® (270 d.). The vertical bars shown indicate 
the standard deviations. Where no bars are shown the standard deviations are about the same 
as, or less than, the magnitudes of the symbols marking the experimental points. 


About 10% of the spectrum can however be assigned to a component of half- 
life longer than 270 days without being inconsistent with the experimental 
data, but there is no evidence for a component with a half-life as short as 
47.9 days. 

From the above it appears certain that the spectrum observed is due largely 
if not entirely to 270 d. Ag!!° contained as an impurity in the sample material. 

From the results shown in Fig. 4 an upper limit of 3 X 10~‘ 8-rays per disinte- 
gration can be placed on the occurrence of the 473 kev. 8-ray transition. This 
figure is based upon a K conversion coefficient of 0.13+0.01, and K/(L+M) = 
2.75, as measured by the author. The log fot value corresponding to this limit 
is 11.3. If the transition is of the AJ = 2 (yes) type (Section 4), then, using the 
curves of Davidson (1951), a value of > 10.2 is obtained for log ft. 

Asummary of the results of these investigations is given in Table I. 


TABLE I 
Energy of §-ray Intensity of 
Nucleus transition (kev.) 8-ray transition log fot 
Cs'™ (2.3 yr.) 1451 <5x10- >14.5 
Co® (5.3 yr.) 1478 (1.0+0.2) X10~ 14.0+0.1* 
Sc** (84 dp 1475 (3.6+0.7) X10-* 13.0+0.1 
Hg*® (47.9 d.) 473 <3x<10~ >11.3t 


~ *log fet = 13.6+0.1. 
tlog fit > 10.2, if spectrum is of AJ = 2 (yes) type. 
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4. DISCUSSION 


Cs'*4 (2.3 yr.) 

The transition sought in the decay of 2.3 yr. Cs!*4 is that proceeding from the 
ground state of Cs'*4 to the 605 kev. state in Ba!*4 with an energy of 1451 kev. 
A 2+ state in Ba'* at 605 kev. has been proposed by both Cork et al. (1953) 
and Keister et al. (1955). Further, Klema (1955) has just reported that experi- 
ments conducted at Oak Ridge are in agreement with the disintegration 
scheme proposed by Keister et al. However the possibility remains that the 
2+ state in Ba'*‘ is of energy 797 kev. rather than 605 kev., in which case the 
transition sought, of energy 1259 kev., would not have been detected, since 
examination of the spectrum was restricted to energies greater than that of 
the 1367 kev. y-ray. 

The spin of the ground state of Cs!*4 has been determined by Bellamy and 
Smith (1953) to be four, and from shell structure considerations its parity is 
expected to be even. The §-ray transition to the 2+ state in Ba!*4 is therefore 
expected to be characterized by a spin change of two without parity change; 
i.e. it is a transition of the second order of forbiddenness. 

From what has been said, the upper limit of 5 X 10-5 6-rays per disintegration 
for the 1451 kev. transition is significant only if a 2+ state in Ba!*4 exists at 
605 kev. However, data shown in Fig. 1 are otherwise useful in demonstrating 
the extremely low level of scattering within the spectrometer. The counting 
rate at energies just greater than that of the y-rays shows no contribution 
whatsoever from conversion electrons scattered into the detector. This obser- 
vation supplements previous findings (Elliott et al. 1954) of the excellent 
performance of the spectrometer with respect to background levels in general. 


Co® (5.3 yr.) 

The transition sought in the decay of Co® (5.3 yr.) proceeds from the ground 
state of Co® to the 1332 kev. level in Ni®*. The energy of the intense transition 
to the 2504 kev. level in Ni® has recently been measured by Bolla et al. 
(1954), who obtain a value of 309+3 kev., whence the transition to the state 
in Ni® at 1332 kev. has an energy of 148143 kev. From shell structure, the 
parity of the ground state of Co® is expected to be even. An assignment of 
spin five to this state is consistent with most experimental evidence (Gold- 
haber and Hill 1952). An assignment of 2+ to the 1332 kev. level in Ni® 
is in accord with all experimental evidence (Goldhaber and Hill 1952). The 
1481 kev. 8-ray transition is therefore expected to be characterized by a spin 
change of three without change of parity. The spectrum shape will accordingly 
be unique (Wu 1955). 

Keister and Schmidt (1954) have recently reported observing this spectrum 
with an intensity of 0.15% of the total number of disintegrations. The shape 
of the spectrum obtained by them does not however correspond to a spin 
change of three, but rather to a spin change of two. To account for this they 
propose a disintegration scheme for Co® in which the ground state of Co® 
is assigned a spin of four and the 10.7 min. state of Co® a spin of one, both 
with even parity. This latter proposal is not in accord with the failure of 
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Deutsch and Scharff-Goldhaber (1951) to observe a 6-ray transition from the 
10.7 min. level in Co® to the ground state of Ni®. 

The results obtained in this investigation differ from those reported by 
Keister and Schmidt (1954). The intensity of the transition found by these 
workers, 0.15% of the disintegrations, is more than an order of magnitude 
greater than that found in this study. Further, over the restricted range of 
the spectrum here studied, the shape of the spectrum has been found to be 
consistent with a AJ = 3 (no) transition, in agreement with the requirements 
of all other experimental results (Deutsch and Scharff-Goldhaber 1951; 
Goldhaber and Hill 1952). 

The spectrum of the secondary electrons does not affect the results obtained 
in the present work since measurements were taken only at energies greater 
than those of all y-rays emitted by the source. This however does not hold for 
the investigations of Keister and Schmidt (1954), who extended their observa- 
tions to much lower energies, and it is possible that their spectrum shape and 
deduced intensity may be accounted for by insufficient correction for the 
intensity of the secondary electron component in their measured spectrum. 


Sc*® (84 d.) 

An assignment of 4+ is probable for the ground state of Sc** and an assign- 
ment of 2+ for the first excited state in Ti*® (Goldhaber and Hill 1952). 
The high-energy 8-ray spectrum is therefore expected to be characterized by a 
spin change of two without parity change, and accordingly will be second 
forbidden (Wu 1955). 

Keister and Schmidt (1954) and Whalen et al. (1953) have reported the 
detection of a B-ray component of energy about 1.25 Mev. with an intensity 
of 0.1% according to the former, and <0.1% according to the latter. 

The results obtained here differ from those obtained by Keister and Schmidt 
(1954) as to both intensity and energy. The intensity of the high-energy 
B-ray component as found by Keister and Schmidt is about 30 times greater 
than that reported here. The energy of the transition given by them, viz. 1.25 
Mev., requires that the 2+ state in Ti*® be at 1.12 Mev., whereas the energy 
of the transition found here, 1475 kev., requires the 2+ state to be at 890 kev. 
In this connection, Heydenburg and Temmer (1955) have recently reported a 
value of 890 kev. for this state in Ti‘*, as found in Coulomb excitation experi- 
ments. The remarks made above regarding Keister and Schmidt’s results 
for the high-energy spectrum detected by them in the decay of Co® (5.3 yr.) 
apply also to their results in the present case. 


Hg? (47.9 d.) 

The transition sought in the decay of 47.9 d. Hg?®* proceeds from the ground 
state of Hg*®* to the ground state of Tl?°*. The main transition is from the 
ground state of Hg?, via a 195+2 kev. 8-ray transition, to a 278.3+0.5 kev. 
state in Tl*°*. Both of these measurements have been made by the author. 
The energy of the high-energy 8-ray transition is therefore 473.3 kev. From 
shell structure considerations the odd neutron in Hg? is in either a 3/2 
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or an fs,/2 state. The parity of the state is therefore expected to be odd. The 
parity of the ground state of T1?°*, from shell structure, is even, and the spin 
of the state has been measured to be 1/2 (Jackson 1932). The 6-ray transition 
to the ground state of T1*°* is therefore expected to involve a change of parity, 
and possibly a spin change of two. In this case the spectrum has a unique 
shape (Wu 1955). Wapstra et al. (1954) have recently reported an upper limit 
of 1.5% for the occurrence of this transition. 

The lower limit of 10.2 obtained here for log f:t of the transition, while 
somewhat high, is nevertheless not inconsistent with a AJ = 2 (yes) transition. 
The presence of the 270 d. Ag!!° contaminant unfortunately rendered the 
experiment less sensitive than it would otherwise have been. 
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AN EXAMINATION OF VERWEY’S MODEL FOR THE 
LATTICE STRUCTURE OF THE FREE SURFACE 
OF ALKALI HALIDE CRYSTALS! 


By G. C. Benson, H. P. SCHREIBER,? 
AND D. PATTERSON? 


ABSTRACT 


The Verwey model for the free surface structure of alkali halide crystals at 
0° K. has been used to estimate the effect of surface distortion on the surface 
energy of these crystals. It appears that the contributions from distortion are 
considerably larger than previously reported. In the case of some lithium salts 
the corrected surface energy becomes negative; the mechanical stability of the 
solid is considered in the light of this. Several aspects of the Verwey model have 
been examined. These include the effect of various repulsive potentials and forms 
of distortion. 


1. INTRODUCTION 


The Born—Mayer lattice model has been widely adapted for the compu- 
tation of surface energies, particularly for the alkali halides. The problem is 
complicated by the presence of an unsymmetrical electric field at the crystal 
surface, which gives rise to some form of surface deformation. The effect of 
these deformations on the surface energy has been discussed by various 
authors (Lennard-Jones and Dent 1928; Dent 1929; Moliére, Rathje, and 
Stranski 1949; Glauberman 1949). Perhaps the most realistic model to date 
has been proposed by Verwey (1946) who, following Madelung’s (1919) sug- 
gestion, has examined the displacements of the positive and negative surface 
ions of alkali halides perpendicular to the surface plane. Polarization of the 
surface ions has also been taken into account, and a lowering of the surface 
energy due to surface relaxation predicted. As more accurate experimental 
values of the surface energies of certain alkali halides may be obtained in the 
foreseeable future, it seems timely to re-examine the theoretical work. 

This paper is divided into the following sections. First, an expression is 
developed for the change in energy of an infinite crystal as it is separated 
along a plane into two semi-infinite hemi-crystals, the surfaces of the hemi- 
crystals being allowed to deform in the Verwey manner. The development of 
this expression seemed justified by a first examination of the Verwey model, 
which yielded a surprisingly large lowering of the surface energy of NaCl and 
NaBr. In fact it is shown later that in the case of LiF and Lil, the final surface 
energy is negative. This unexpected result may be due to the inapplicability 
of the model, but there did not seem to be an a priori reason to reject a negative 
surface energy as impossible, apart from the obvious consideration of mechan- 
ical stability of the crystal. The expression, though not necessary to our main 
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argument, is useful in discussing the plausibility of the negative surface energy; 
thus, it seems worth while to present the mathematical analysis attending its 
derivation. The normal surface energy as discussed by Verwey can be obtained 
by considering the change in energy when the two hemi-crystals are at an 
infinite distance. In the second section, the ionic displacements calculated by 
Verwey are checked for NaCl and NaBr, and the corresponding distortion 
energies and surface energies calculated and discussed. In the light of the 
results of this section, a third section is devoted to an examination of the 
effects on the surface energy of assuming various possible forms of the re- 
pulsive potential between ions, and of assuming various displacement and 
polarization characteristics. Finally, the distortion and surface energies of 
LiF and Lil are calculated, and the results discussed with the aid of the above- 
mentioned model. 


2. THE CHANGE IN ENERGY WHEN AN INFINITE CRYSTAL 
IS SEPARATED INTO TWO HEMI-CRYSTALS 


The model (shown in Fig. 1) assumes that the Verwey distortion is restricted 
to the first layer S of each hemi-crystal. The structure of the inner layers J is 
regular. The lattice nearest neighbor distance a is taken as a convenient unit 


© 2 Os r= 2 Oe 


Fic. 1. Model for the calculation of the change in energy when an infinite crystal is sepa- 
rated into two hemi-crystals. 


of length. At any separation na, measured between the positions of the un- 
‘distorted first layers of the hemi-crystals, the displacements of tne positive 
and negative ions of the left-hand surface layer are* z,a and z_a; these have 
positive value if the ions have moved to the right of their undistorted positions. 
The corresponding displacements of the ions of the right-hand surface are 
then —z,a and —z_a. For simplicity, only the polarization of the negative 
ions of the surface layers is considered. These ions have moment y in the 
direction Z for ions in surface 1, and —wp for ions in surface 2. 


*Our notation 24a corresponds to Verwey's displacement 2;. 
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For n = 1, the two hemi-crystals coalesce to form a perfect infinite crystal. 
Here all displacements and dipole moments must vanish. At any separation 7, 
U(n) is the total energy of the system relative to the free ions and stated per 
ion pair in the combined surface planes of parts 1 and 2. This quantity is of 
course infinite, but for any two values of the separation a finite change in it 
can be computed. Adopting the regular infinite crystal as a reference state, 
the change in energy is denoted by AU(n), where 


(1) AU(n) = U(n)—U(1). 


The energy U(n) is computed by summing interactions between pairs of ions 
and pertinent internal energy contributions from the ions. In general, the 
interaction energy for two ions i and j with charge e; and e; and dipole moments 
ui and y; is taken to be 


(2) ui = Ce fiy —eu(Liy wyriy Hes (Ci wots 
Dotan B(Fiswe)E es a) 
BU; ( 143 44 +R(ri;), 


where r;; is the position vector of the jth ion relative to the ith, and R(r;;) is 
the repulsive energy. The dipoles are considered to be permanent and subse- 
quently a correction for the quasi-elastic energy is introduced as an internal 
energy contribution. 

It is convenient to further divide U(n) into two parts: Uy2.(n) and U1,(n). 
The former is composed of terms representing the interactions of ions in part 1 
with those in part 2; the latter represents interaction terms between pairs of 
ions in either of the two parts together with all pertinent internal energy 
contributions. Hence we can write 


(3) AU(n) = AU ii(n) +AU 12(m) « 


Because of the symmetry when the surface is a {100} plane, Ui:(n) can be 
computed by considering terms arising from part 1 alone. Stating this per ion 
pair in the surface of the hemi-crystal 


(4) AUn() = Un(n)—Un(1) 

n MN Si Si 2 

be U4 s+ b us+3 do ty s+4 D0 u_s+5- ’ 
j J J J 2a 


where the first two sums represent respectively the interactions of a positive 
and negative surface ion in S; with all ions in J;, the next two sums are 
analogously defined for interactions between surface ions, and the last term 
is the quasi-elastic energy of a negative surface ion with polarizability a. 

The electrostatic lattice sums in equation (4) can be treated by the method 
of van der Hoff and Benson (1953), leading to the equation* 


*It may be noted that 7 does not appear explicitly in equation (5), but implicitly through 
the displacements and dipole moment. 
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(5) AUiu(n) = Tints 5 (z_—2,)° 
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where T = 9.0336217 is a constant representing a lattice sum originally com- 
puted by Topping (see van der Hoff and Benson 1953). Terms in (z_—2,) of 
higher than square order have been neglected. 

In evaluating the repulsive contributions ARi:(n), it is assumed for the 
moment that the Born—Mayer form is applicable, so that R(r;;) may be equated 
to b:; exp(—ri;/p) = bei; exp[(ri+r;)/p] exp(—rij/p). Here r; and r; are the 
two ionic radii and the c;; are the Pauling factors for ++, +—, and —— 
interactions. These have the respective values 2.00, 1.375, and 0.75 for lithium 
halides, and 1.25, 1, and 0.75 for other alkali halides. The value of p is 
0.345 X 10-8 cm. The value of 6 is determined for each alkali halide by mak- 
ing the spacing at the minimum of the total potential energy curve for the 
uniform dilatation of the crystal agree with the experimental value of the 
nearest neighbor distance. Following Verwey, only repulsions between nearest 
neighbor ions and next-nearest neighbor negative ions are considered. 

Thus, 


(6) ARiu(n) = by-[exp[—(a/p) (1+24)]+exp[— (a/p)(1+2_)] 
+4 exp{—(a/p)V[1+ (z_—2)]} —6 exp(—a/p)] 
+4b_[exp{ — (a/p)V/[1+ (1+2_)’]} —exp[— (a/p)v2]]. 


Some modifications of the repulsive contribution will be discussed later. 
In a similar way, it may be shown that 


eas 
(7) AUi2(n) = - [sx D1 V (i +12"). {1 +exp[— av (11° +/2°)]} 


xf expl— #(1-+9—a5)V (i109) 
+exp[—x(1+n—2.)V(h' +h*)]—2 exp[—24v (hi +1’) 
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—2[exp[— w(n—24—2_)V (li +12”)]—exp[— rv(ui+iyy} | 


+4x a of Xr Lo lexp[— a(n—24.—2_)V (ly +/2')] 
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—(—1)"exp[— #(9—22_)V (L'+h’)]] 


_ » exp[—x(1+n-2)V(h +H] 
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+2x° a a Lo (- 1) "J (h?-+12") exp[— #(n—22_)V (Li? +12’) ]+ 4Ra(n), 
li+ leeven 
where 
ARi2(n) = b4_[exp[—(a/p)(n —24 —2_)] —exp(—a/p)] 
+2b__[exp{ —(a/p) V[1+(n—2z_)*]} —exp[—(a/p)V2]] . 


Following Born and Heisenberg (1924), the variables z,, z_, and yz in the 
foregoing equations may be treated as independent. Then for any fixed 7, the 
conditions for equilibrium are 


9AU(n) _ 9AU(n) _ 9AU(n) _ 


O24 Oz_ Ou 


A solution of these conditions yields the displacements and dipole moment. 
As a check on internal consistency it can be shown that at 7 = 1, the equations 
are simultaneously satisfied by z; = z_ = u = 0, and also that at this point 
dAU(n)/dn = 0. 

At infinite separation U,2(n) vanishes and 
(8) AU}:(@) = — U,,(1) ° 


Equation (8) is the surface energy per ion pair of the undistorted hemi- 
crystal. The term AU;;(@) then is the ‘distortion energy’’, representing the 
correction to the surface energy arising from the displacement and polar- 
ization of the ions in the surface layer. 

For the purpose of numerical calculations the electrostatic sums were simpli- 
fied in a manner essentially equivalent to that adopted by Verwey. Thus, 
only the first term of each sum was considered and exp[—2+/(/,?+/,")] was 
neglected when compared with unity. These two simplifications have a com- 
pensating effect. The simplified expression for AU1;(n) is 


AUuls) © él (-s35) F (—s)" 
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A similar simplified form of AUi2(n) can be derived. At infinite separation, 
equilibrium conditions demand that 
dAUi;(@ ) dAU (©) dAU (©) 
ee BE renee SE : eee ernee enn ae 0, 
O24. Oz Ou 


giving rise to the equations 


(a) 0 = dAUui(@) _ 7 Th. —13 


O24 a 2/2 
48x expl-#(1+5,)v2] | 
e <* a OARi1(@ ) 
— Ali 5) as 


(10) (6) O= gakule) = 2) T (2-—24) 


+87 exp[— r(l+s)v2 | 


+5 | G2) T+82°4/2. exp[— r(1+s-)v2 | 


OARi(”) 
+ 0z_ , 
_ 9AUn(@) _ _e o_o = 
(c) O= .. = (i 273 T(2_—2;) 


+8 exp[— (1 +e)vai feel oot] i 


The first two conditions have been obtained by Verwey by separately equating 
to zero the net force on a positive and negative surface ion, while the last 
equation corresponds to nu = aF, where E is the total electric field at a negative 
surface ion. Equation (10) (c) may be solved for u as a function of 2, and z_, 
and the result employed in (10) (a) and (6). The numerical solution of these 
equations is tedious as Verwey has noted; in this work values of 2, and z_ 
were estimated and refined by the Newton-Raphson (Whittaker and Robinson 
1926) technique employing the derivatives of equations (10) (a) and (0). 
3. A CHECK OF VERWEY’S DISPLACEMENTS FOR NaCl AND NaBr. 
THE SURFACE ENERGIES OF THESE SALTS 
In order to carry out this computation in a manner analogous to Verwey, 
a further adjustment of the repulsive interaction terms had to be made. For 
surface displacements in which the interpenetration of ions was such that 
r<r,+r_, the exponential power law was replaced by a ‘‘harder”’ inverse 12th 
power law (i.e. R(r) = A+B/r'*). The transition was effected by matching 
the force and potential statements of the two expressions at r = r,+r_. The 
adjustment constant A was added in order to make the two terms identical 
at the chosen separation. The additional repulsive law constants were ob- 
tained from the equations 
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a tae by exp[—(ra+r_)/o], 
aoa —= exp[— (¢-+r-)/a). 


With the aid of equations described in the preceding section, the displace- 
ments of the surface ions of NaCl and NaBr were evaluated for 7 = ~. The 
distortion energies (AU1:(~)) and the surface energies were also obtained. 
The surface energy of a {100} surface at any separation 7» of the two hemi- 
crystals is defined by 
(12) o100(n) = AU(n)/2a* 


where the factor 1/2a? converts the result to a unit area basis. o;99( ©) is then 
the normal surface energy of a {100} face. These results together with the 
values obtained by Verwey are shown in Table I. Values of the pertinent 
constants are also shown. 

TABLE I 
DISPLACEMENTS OF THE SURFACE IONS, DISTORTION AND SURFACE ENERGIES OF NaCl anp NaBr 








AUi(2) a100( 2) 


a ¥: a a, By 3 
cm. X 107-8 cm.*X 107" 25 z. debye erg-cm.~? 
NaCl This work 0.95 1.81 2.81 3.1 —.0586 .0173 0.9348 -—77.1 52 
Verwey 0.98 1.81 2.81 3.1 —.0534 .018 0.90 — —_ 
NaBr This work 0.95 1.95 2.97 4.2 —.0701 .0205 1.2082 —83.9 32 
Verwey 0.98 1.95 2.97 4.2 —.0700 .0212 1.20 —32 90 
(—84) (38) 


The agreement between the two sets of displacements is good,* but the 
energy evaluation differs in the case of NaBr by about 50 erg/cm.? We are 
most indebted to Dr. Verwey, who has kindly re-evaluated his result for NaBr 
and has reported in private communication a distortion energy value of 
—84 erg/cm.? (bracketed terms in Table 1) in agreement with our result. The 
surface energies are not in complete agreement because of differences in the 
formulation of the electrostatic term. Verwey simplified this to be proportional 
to exp(—72/2) while in this work two additional terms were retained 
(i.e. exp(—2V/2) —2 exp(—24+/2)+exp(—32r+/2)). The surprisingly large 
distortion energies lower the values of the surface energies beyond the scope 
of the discussion at the end of Verwey’s paper. In the case of NaCl the value 
is much smaller than that indicated by experimental data (Benson and Benson 
1955) on the heat of solution of small particles. 

It must also be pointed out that the approximate corrections used by Shuttle- 
worth (1949) for the distortion energy of the alkali halides and based on Ver- 
wey’s NaBr result must be regarded as erroneous. 

_ “Exact agreement could not be expected because of the different values of the ionic radius 


of Nat used. For consistency, the two Pauling radii were used in this work rather than the 
Goldschmidt radius for Na* used by Verwey. 

+A direct comparison with these results is not valid since the experimental work was done 
at 298° K. while the theoretical model pertains to a crystal at absolute zero. However a 
rough estimate of the temperature correction which can be obtained from the data of 
Patte rson, Morrison, and Thompson (1955) indicates that this effect is relatively small. 
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4. DEPENDENCE OF SURFACE ENERGY ON MODEL CHARACTERISTICS 


Verwey has pointed out that the magnitude of the distortion energy is 
sensitive to the choice made for the function representing repulsive inter- 
actions between the ions. Furthermore, all other models used to describe the 
surface configuration lead to absolute values of the distortion energy sub- 
stantially smaller than we obtained. The dependence of the distortion and 
surface energies on repulsive interaction and displacement characteristics has 
therefore been investigated. 

(a) Displacements, distortion and surface energies for NaCl and NaBr 
were computed for three forms of the repulsive interaction. The data utilizing 
the inverse 12th power relationship have already been presented. Verwey’s use 
of the inverse 12th power law was justified by the work of Verwey and de Boer 
(1940) where, in formulating the molecular energy of the diatomic alkali 
halides, the inverse power law was found best suited to represent repulsive 
interactions when the ions were separated by a distance smaller than the sum 
of the ionic radii. In a more recent paper, Rittner (1951) has shown preference 
for an altered exponential power law utilizing a value of p = 0.31 X10-* cm. 
A second set of data has therefore been computed using a combination of two 
exponential repulsive power laws which were matched in the same way as 
described in Section 3. Finally the single exponential form (with p = 0.345 
X10-* cm.) was used to represent all repulsive interactions. The results are 
given in Table II. 


TABLE II 


DISPLACEMENTS, DISTORTION AND SURFACE ENERGIES FOR NaCl AND NaBr 
USING THREE FORMS OF REPULSIVE INTERACTION 


py AU) 7100( ©) 








Repulsive terms 2+ z_. 
erg-cm.~? 
NaCl Exponential + inverse 12th —.0586 +.0173 — 77.1 52 
Two exponential —.0938 —.0059 — 89.7 39 
Single exponential —.1105  —.0164 — 96.1 33 
NaBr Exponential + inverse 12th —.0701 +.0205 — 83.9 32 
Two exponential —.1034 —.0002 — 96.1 20 


Single exponential —.1201  —.0103 — 103.2 13 


Although the energy values are rather sensitive to the choice of repulsive 
law function, none of the three methods of evaluation leads to an appreciably 
lower distortion energy. The present results then are more likely characteristic 
of the surface configuration than of any particular combination of repulsive 
interaction expressions. 

(6) The dependence of the surface energy on surface configuration was 
investigated for the case of NaCl only. For simplicity, repulsive interactions 
were uniformly represented by the single exponential form. This seems justified 
because there is as yet no clear statement of preference for any one method of 
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formulation; moreover, the effect discussed in the preceding subsection 
appears to be of secondary magnitude. The following surface configurations 
were considered :* 


(i) Positive and negative surface ions in a plane, which is free to move in 
Z-direction. No polarization. (2, = z_, u = 0.) 
(ii) Surface ions move separately. No polarization. (2; # z_, u = 0.) 
(iii) No displacements. Negative ions polarized. (2; = z. = 0, w_ finite.) 
(iv) No displacement. Both types of surface ions polarized. (2, = z_ = 0, 
uy and pw finite.) A value of ay, = 0.1810-*4 was used. 


(v) Same as (i) but negative surface ions polarized. (2; = 2_, u— finite.) 


(vi) Same as (i) but both types of surface ions polarized. (Lennard-Jones 
and Dent model, 2, = z_, un, and yw finite.) 


(vii) Verwey model (2, ¥ z_, u— finite). 
The results are shown in Table III. The displacement in model (vi) is about 


TABLE III 
DEPENDENCE OF THE SURFACE ENERGY OF NaCl ON SURFACE CONFIGURATION 


53 4Un(©) 7100( ©) 


Model £5 Z ra by 
debye erg-cm.~? 
(i) —.0052 —.0052 0 0 0 129 
(ii) —.0193 —.0022 0 0 — 4 125 
(iii) 0 0 0.385 0 —22 107 
(iv) 0 0 0.414 —0.051 —24 105 
(v) —.030 — .030 0.440 0 —26 103 
(vi) —.0355  —.0355 0.484 —0.059 —27 102 
(vii) —.1105 —.0164 1.131 0 —96 33 


half that given by Lennard-Jones and Dent (1928). This confirms criticism of 
their calculations by Moliére, Rathje, and Stranski (1949). It is evident from 
the table that it is the combination of ion polarization and separate displace- 
ment which is instrumental in creating the large distortion energy. Thus the 
Verwey model appears to lead to the most stable surface structure of NaCl. 
This must, however, be regarded only in the light of the restricting assump- 
tions inherent in the model. Inclusion of the van der Waals attractive forces 
for example may greatly alter these results. Further, no final conclusions can 
be reached until the problems of a more realistic formulation of the polariza- 
bility of ions in crystal lattices, and the connection between the repulsive term 
and the state of polarization of the ions, have been resolved. Mayer (Gomer 
and Smith 1953) has recently drawn attention to the importance of these 
questions. 


*Equations (9) and (10) of Section 2 must be modified according to the needs of each model. 
In cases where y, is finite, additional appropriate terms must of course be added to equation (9), 
and equations (10) augmented by a fourth equilibrium condition. 








274 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 
5. DISTORTION AND SURFACE ENERGIES OF LiF AND Lil 


The importance of the distortion energy contribution to the surface energy 
appears to increase as the size difference between constituent ions increases. 
It was therefore interesting to calculate the distortion and surface energies 
of LiF and Lil. Verwey’s values of the displacements and dipole moments were 
employed, and the inverse 12th power repulsive relationship applied under 
the appropriate conditions. The distortion energies of the halides were: 
—203 erg/cm.? for the fluoride, —237 erg/cm.? for the iodide, leading to re- 
spective surface energies of —93 and —142 erg/cm.? 

The negative surface energies could be taken to indicate the inapplicability 
of the classical Born—Mayer model to surface energy computations in these 
alkali halides. This may be a result of the inherent oversimplifications men- 
tioned at the end of Section 4. The present results must in any case be con- 
sidered with much caution. Not only have the important van der Waals inter- 
action terms been omitted, but, in order to conform with Verwey’s calcu- 
lations, positive—positive ion repulsive interactions and third nearest neighbor 
repulsive terms have been neglected, although a crude calculation shows these 
terms to be still appreciable in the cases in point. It seems fairly certain, 
however, that inclusion of the additional repulsive terms would not alter the 
sign of the surface energy. 

It has been pointed out in the introduction that no a priori reason exists for 
rejecting negative surface energies, if mechanical stability of the crystal can 
be maintained. The surface energy of LiF has therefore been computed for 
two small separations. In the course of this computation it was found that 
Verwey’s displacements at 7 = @ did not satisfy the equilibrium conditions 
as well as might be expected. A new set of values was therefore computed and 
tabulated in Table 1V together with values of the surface energy as a function 
of n. It will be noted that at » = @ the change in distortion energy from that 
based on Verwey’s displacements is trivial. 


TABLE IV 
SURFACE CHARACTERISTICS AND SURFACE ENERGY OF LiF AS A FUNCTION 
OF CRYSTAL SEPARATION 


n 1.0 i 1.5 Co 











24 0 .0323 — .0686 — .0697 
Z 0 .0326 — .0039 +.0014 
p_, debye 0 .0037 .5198 .5616 
. AU i1(n), erg-cm.~? 0 7.0 — 206 —207 
2a? 
1 - - 
dq? 4U:12(n), erg-cm.~? 0 0.2 109 110 
a100(n), erg-cm.~? 0 te — 97 — 97 


Initially there is an increase in the energy which goes through a maximum 
for a separation of n in the vicinity of 1.1 and then decreases rapidly to its 
final negative value in the region 7 = 1.5. It is possible, therefore, that the 
initial increase of energy serves to preserve the mechanical stability of the 
infinite crystal at 0° K. 
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In summary it appears that the Verwey model leads to surface energies of 
the alkali halides which are smaller than can be justified by the few experi- 
mental data available or by extrapolation from the liquid state. Undoubtedly 
the picture of separate displacement of positive and negative surface ions is a 
more realistic approach than used by previous workers. However, quantitative 
results for the distortion and associated energy change calculated by the 
classical Born—Mayer method with simple polarization corrections must be 
considered with some caution. It is suggested that the difficulty is a funda- 
mental one and requires a more detailed treatment of the polarization and 
repulsive terms. 
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THE RADIATION FORCE ON A SPHERICAL OBSTACLE IN A 
CYLINDRICAL SOUND FIELD! 


By T. F. W. EMBLETON 


ABSTRACT 


The corresponding problems of the radiation forces acting on a spherical 
obstacle in a plane or spherical progressive sound field have been examined 
previously by King (Proc. Roy. Soc. (London), A, 147: 212. 1934) and Embleton 
(J. Acoust. Soc. Amer. 26: 40. 1954) respectively. In these cases the sound field and 
scattering obstacle both have symmetry of revolution about the line joining the 
center of the obstacle to the origin of the sound field, but for a cylindrical sound 
field there exists only a much lower degree of symmetry. A general expression 
has been obtained for the radiation force in terms of the complex amplitudes 
of spherical harmonics required to synthesize the incident sound field—for the 
cases of greatest symmetry this reduces to the simpler expression previously 
obtained. The first 20 non-zero amplitudes have been evaluated for a cylindrical 
sound field and it is shown that the force is one of attraction near te the source, 
becomes zero at a certain distance, and is a force of repulsion at a greater distance. 
Qualitatively, this force is the same as for spherical waves but for any size of 
obstacle and frequency of the sound field the point of zero force is always nearer 
to the source in a cylindrical wave. 


I. INTRODUCTION 
The importance of dealing theoretically with the radiation force in a sound 
field where there is no simple relation between the mean kinetic and mean 
potential energy densities has been dealt with in a previous paper (Embleton 
1954) where the field considered was spherical. This present paper deals with a 
cylindrical sound field where, in addition, there is a lower degree of symmetry 
of the field and consequently terms remain here which have vanished from 
previous analyses. The method used is similar to that employed by King 
(1934) who considered a spherical obstacle in a plane field—other work using 
essentially King’s method but of no direct connection with the present paper 
is that of Awatani (1955) who derived the radiation force of plane waves on a 
circular disk. 
II. SYMBOLS 
6p pressure variation from undisturbed state. 
po density of propagation medium in undisturbed state. 
t time coordinate. 
¢@ velocity potential. 
@,, p1, M_ radius, mean density, and mass of spherical obstacle. 
c phase velocity of sound. 
p, 0,2 spherical coordinates; equilibrium position of center of sphere as 
origin, O; the line OL from O perpendicular to the line source as 
the zero direction of 6; the plane containing the line source and O 
as the zero plane of Q, see Fig. 1 
r,l distance of any element 6/ of line source from the center of the 
sphere O and the point L respectively. 
ro mean distance of the center of the sphere from the line source 
(the length of the line OL). 
‘Manuscript received November 14, 1955. 
Contribution from the Division of Applied Physics, National Research Council, Ottawa, 
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III. THEORY 
(a) Boundary Conditions 
The pressure variation in the propagation medium (most commonly, air) 


may be written entirely in terms of the velocity potential using the form of 
Bernoulli’s equation applicable to non-steady flow 


(1) ap = pods" 6-5 polV4)* 

This equation is correct to the order of (V@)?/c? and so ¢ may be expressed 
as any solution of the simple wave equation, (V?+k?)¢ = 0, where the wave 
number k defines the frequency of the sound field. Equation (1) may be 
integrated over the surface of the sphere to obtain the net force acting on the 
sphere. The last two terms represent a unidirectional force when integrated 
over the surface; the first term will give rise to dynamical equations of motion 
of the sphere if it is free to move. If u is the linear velocity of the sphere 
determined in this way with respect to the origin O of the velocity potential 
field, then in place of ¢ we must have ¢—u.V¢. The last term of this expres- 
sion gives rise to a contribution to 6p of the same order of magnitude as the 
last two terms of Equation (1), which now becomes 


. ] . 1 
(2) 5p = pod—pou.V o+5 ° e-5 po(V >)’. 


From the degree of symmetry of the cylindrical sound field and the spherical 
obstacle the only resultant forces acting on the sphere will be along the line 
OL, Fig. 1. It is assumed that the sphere is rigid, that its surface is impervious 


P( p,8,2) 





Fic. 1. 


to the sound field, but that it is free to move. The sphere will execute small 
oscillations along the line OL under a law of force given by 


(3) Mu = —2a; f J pod sin 8. cos 6. d0. dQ. 
0 Jo 
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Integrating Equation (3) with respect to time, the oscillatory motion of the 
sphere is obtained as 


° 1 
(4) Mu = —201%90 f J on -du..ao 


where » = cos 0. 

Since the fluid medium must remain in contact with the sphere at all points 
on its surface, the velocity potential must satisfy the boundary condition 
that 


(5) -(2¢) = u cos 6. 


(b) Solution of the Wave Equation 

From the nature of the boundary conditions it is most convenient to adopt 
a solution of the wave equation in terms of spherical coordinates, in which 
the wave equation is separable, based on the center of the sphere as origin. 
The velocity potential may then be written as 


(6) ie ze {AnVa(kp)+An'fa(kp)} » (kp)* - Sa(u, 2) 


where S, is a surface harmonic of order n, a function of both @ (or uw) and Q. 
¥n(z), ®,(z), and f,(z) are spherical wave functions equal to (—d/z dz)" operat- 
ing on (sin z)/z, (cos z)/z, and e~‘*/z respectively. From the nature of these 
functions the term 


(7) — 2 Anva(be) . (kp)".Sn 


represents the incident velocity potential and the coefficients A, may be made 
to satisfy the conditions of a line source whose nearest point to the origin is 
at a distance ro. Equation (6) can be made to satisfy Equation (5), the bound- 
ary condition over the surface of the sphere, by suitable choice of the coeffi- 
cients A,’ of the scattered waves in terms of the coefficients A, of the incident 
waves. This has been described in sufficient detail in previous papers (King 
1934; Embleton 1954); here it is only necessary to state that A,’ was written 
as 


i G, (a) 
F,, (a) —1G, (a) 


where a = ka,. The form of F,(a) and G,(a) was then determined in order to 
satisfy the boundary condition. H,(a) and ¢, which are required later may 
now be defined from F, (a) and G,(a) by 


(9) F,(a)+i G,(a) = H,(a) .e™. 


The surface harmonic S,(u, 2) may be expanded in a series of terms varying 
as cos m2 and sin m2 (Lamb 1932) and since the velocity potential at any 
point P is expressed in terms of S, over the whole sphere, m may be considered 
an integer. The surface harmonic takes the form 


(8) A,’ = A, 
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(10) Sn(u, 2) = Pr(u) + > (Bn cos mM2+C,, sin m2) P™(u), 
m=1 


where P,,(u) is a Legendre polynomial and Py (u) are the associated functions. 
The sound field originating from a uniform line source lying in the zero plane 
of the angle 2 is symmetrical with respect to this plane and hence the non- 
symmetrical terms in 2 must vanish. Hence, the expansion (10) reduces 
to 


P,(u)+ x B, cos m2. Pr(u) 


m=1 


or A,S, = = Ax cos m2. Pu). 


m=0 


(11) Sa 


Thus the incident sound field of Equation (7) becomes 
(12) $i = 2 dy An. ¥n(kp) . (kp)" . cos mQ . Pr'(u). 
(c) Analysis of the Radiation Force 
In general the coefficients Ay of the incident sound field are complex and 
may be written 
(13) An = |Az| exp[i(wtt B,')]. 
Using Equations (8), (9), (11), and (13), Equation (6) may be put in the 
form 


(14) ¢= ona > R* cos m2 . P7(u)+sin a - T rcos m2. P,"(u), 








n=0 m=0 n=0 m=0 
where 
. |An| cos(Br'+én) _ |An| sin (Bn'+€n) 
(15) Ri; = H, (a) . a a and y gy = H, (a). oe 


When Equation (14) is substituted into Equation (2) an expression is 
obtained for the instantaneous pressure 6p at any point on the surface of the 
sphere. This may then be averaged with respect to time to obtain the mean 
pressure at any point and integrated with respect to @ and Q to obtain the net, 
mean force on the whole sphere, i.e. the radiation force. 

The first-order pressure variation P,; (from the first term of the right-hand 
side of Equation (2)) is 


(16) P, = puck} cos wae a Tx cos m2. P (uw) 


n=0 m=0 


—sin on = Rt cos m2. . Pau) 


n=0 m=0 


The time average of this expression is zero, so P; contributes nothing to the 
mean pressure at any point on the sphere. 

Let P, denote the contribution of the term 4(po/c?)¢? when integrated 
over the surface of the sphere and P, its time average. From considerations 
of symmetry P, will act along the line joining O and L. We have 
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(17) Po= ap ff ¢u.du. da 
i 
(18) = —a,'pok” J - {0s at 5 a Tx cos mQ . Px'(u) 
n=0 m=0 


sical Yo RE cos mo. Pru) Yn. dp . dQ. 


n=0 m=0 


Taking the time average of this we obtain 


(19) P, = et J f l} - . cos m2. Pruyt 
ae a Tx cos m2. Pru} | . dp. dQ. 


n=0 m=0 
This expression may be evaluated by making use of the integrals (1) and (2) 
listed in the Appendix, and hence 


B— _ona%o Set) _ (pop? 4 p97 
(20) ss — 27a po Dy (2n+1)(2n+3) (RaRnsit TaTn+1) 
© ! 
— na’po >, > a Ra RnvitTn Tn): 


4=0 Ani (n—m)! (2n+1)(2n+3) 


Similarly let P, denote the contribution of the term — 4} 0(V¢)? when 
integrated over the surface of the sphere and P, its time average. As before 
we have 


® 1 
(21) Py = arto f J. (V>)'u. du. dQ 


cy nf SC) erro SV ats ee 


(0¢/du) and (d¢/d2) may be obtained from Equation (14). When Equation 
(22) is averaged with respect to time, it becomes 


(23) P=" ff [43 RR cos ma-2 Paty) } 


n=0 m=0 


HE > 720s ma® Pa) Joa-0) p.du.dQ 


ep SEs 


HE Y Tim. sin mo. Paw} he a de - a0. 


n=0 m=0 


Rim. sin m2. Pau) + 


. 


s § 


This expression (23) may be evaluated by using the integrals (1), (3), and (4) 
listed in the Appendix, giving 


(24) P= 22d eh (RR, +7973) 


tnt. 2 (n+m+1)! \ 





+7, 2X (n—m—1)! (2n+1)(2n+3) 


> 
m 


X (RIR a+ Ts Tr): 
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Finally let P, denote the contribution of the term —u.Vq@ (which arises 
from the motion of the sphere under the action of the first-order sound field) 
when integrated over the surface of the sphere and P, its time average, 


. 1 2 
(25) P,, = —2a1’p0 J f [ are 28) | ap a 
0 J-1 ay Ou 


The first term of this integral obviously vanishes upon integration. From 
Equations (4) and (14) the velocity of the center of the sphere may be obtained 
as 


(26) “= —{R} cos wt+T} sin ts fe 
a pi 


When this expression is substituted into Equation (25) and (d¢/du) is found 
from Equation (14), we have for the time average P,, 


(27) P,= —% eff OS (RRTETIT™ ) cos m2. (1—p?) . 5 Pa) 


1 n=0 m=0 
.u. dp. dQ. 


Now J cos ma. a2 = 0 for m > 0, 
0 
=m form =0, 


and making use of integral (4) listed in the Appendix, the expression for P, 
becomes 
(28) f= ire o@. - (Ri RI+-T*T?). 
The radiation force on the sphere is i given by the sum of Equations 
(20), (24), and (28) as 
(29) P= ~2ne] © ——_@E)) __1,3_4(64+2)} . (RORY + TOTS | 
—, (2n+1)(2n+3) ° n4\n+ n+ nt+ 
a As (n+m+1)! 
rol 2 Gam)! (n+) Qn+3) 


n=0 m=1 


xchat (n-+2)-+m( 24-41) (RERD+ TITS A) | 











where 
(RaRavitTn Tat 
(30) a al Laas 98 (Bn+1— Bn) (122) Hpyi(a) .ai"*® 
~ ey. H?,:(a) .@ —gnt6| COS(Pn41— Pn) - n(a) «Hoyt : 


4 
~{at—n(n+2)1") +sin (8r41— Bn) « fa*—n(n-+2)} | 


The Equation (30) is valid for all m. Its derivation involves certain properties 
of the H,,(«) functions which have not been given in this paper. In (29) and 
(30), {a?—n(n+2)} is to be replaced by {a?—3(1—po0/p1)} when m = 1. 

It should be noted that Equation (30) and the first half of Equation (29) 
are identical with the complete expressions for the radiation force which exist 
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for the two cases of plane and spherical waves (King 1934; Embleton 1954). 
The second half of Equation (29) has not vanished from this present analysis 
because of the lower degree of symmetry possessed by cylindrical waves. 


(d) Evaluation of the Coefficients |A;| and By 
't now remains to determine |A;| and 6; from Equation (13) after Equation 
(12) has been made to represent a field of diverging cylindrical waves originat- 
ing from a line source whose nearest point to the origin is at a distance of ro. 
Let the line source have a strength of s per unit length so that the velocity 
potential at P due to the element of line source at S is 
— kz 


(31) ds = sl S— = sil folkx). 
This may be written in terms of the spherical coordinates at O (Lamb 1932) 
as 


(32) $s: = sil (2n+1) . (kp)” . (kr)” .fa(kr) . Wa(kp) . Pa(cos POS). 





The velocity potential arising from the whole line source is obtained by 
integrating Equation (32). To simplify evaluation of the angle POS the point 
P is considered to be coplanar with O and the line source, i.e. 2 = 0. Hence 


(33) $, = 2 s(2n-+1)(ke)". Yn (kp) , ” Cer)". fi(kr) 2 ‘i 0+ sin s) dl. 


Equation (33) is identical with Equation (12) for all values of 6 (or uw) and p 
and this identity may be used to evaluate A;’. For each value of n 
n o —,,2)\4 

(34) 3 AZPE(a) = s(n) [” (br) fuer) Py( M24 =H) "a 

m=0 —a 
where r = (ro?+/?)?. The Legendre polynomial on the right-hand side of 
Equation (34) may be expanded by Taylor’s theorem. Since P,(z) is a poly- 
nomial series in z up to 2", the ¢th differential will be zero for ¢ > n. If Py (u) 
is written in differential form, Equation (34) becomes 


Ss ara-s9""( 2)", 


(35) = s(Qn+1) . (kr)*. futhr) | Pa( #22) 


“I! ae d ae (2) | 
* X r' t! d(uro/r) Pr r dl. 

If the coefficients of powers of y greater than yz” are considered, then it can 
be seen immediately that the terms of the left-hand side of Equation (35) 
for which m is odd are associated with, and only with, the terms of the right- 
hand side for which ¢ is odd. Now consider the right-hand side of the equation 
when ¢ is odd. The expansion of 7 is entirely in terms of even powers of / and 
hence the integral is of the form 
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(36) | ” (odd powers of 1) di 


which is evidently zero. Thus all the coefficients AX are zero when m is odd. 

The Legendre functions of Equation (35) may each be expanded as finite 
polynomial series of their arguments. Then, giving » a number of different 
values, a set of simultaneous equations is obtained (more correctly, one set 
for each value of ) from which the value of each AZ may be found. A typical 
coefficient of the incident sound field is 


(37) Ai = & [" (kr)*. fe(kr) . S| se 10% S41 fat 


This integral may most easily be evaluated by converting to Bessel function 
notation. We have 


HH 
(38) fu(er) = —( 2) Hash) 


where H is a Bessel function of the third kind. From Equations (37) and (38) 
we have, writing y = i, 


(2) 
(0) 4t- — af =) f Til bel 


64k [V(x +R 
X [16 (kro) *y’— 16 (kro)*y*+-y*] dy. 


Watson (1944) gives Sonine’s formula for a similar integral of a Bessel function 
of the first kind and an identical formula may be derived for a Bessel function 
of the third kind. This is 


Hy fav (+2°)) sorry, — 2TO+1) yyw 
Ge =r . 7 
fee ce [V(t +z (VM (l+2)) - a's »-1(az) 
provided that a > 0 and Re(}v»—}) > Re(b) > —1. These conditions are 
satisfied by all the expressions similar to Equation (39) for every coefficient. 
The typical coefficient then becomes 
pi St YI iH@ (kro) ,o tH? (kro) Ere) | 

(41) Ag= 64k 16 - (kro) 48 (kr 2 + 15 (kro) Z 
A complete list of the first 20 non-zero coefficients Ay is given in Table I. 

Since iH®(z) = Y,(z)+7J,(z), every coefficient may easily be split into its 
modulus and argument, as required by Equation (13). It is obvious that if 
any coefficient in Table I is written as 





(40) 


(42) An = —""(2n+1) . G[iH? (bra); m, m1), 
then = [AR] = —SF (Qn +1) (FV m, n]+3 (Tes m, nm), 
i: aes BLY; m, n| sechanii 
(48) aii (PVs m, n]+¥ (Jas m, n)}*? 
lJ 1; m, n] 


and sin Br = 1% V3 m, n|t+s Js m, n}}*° 
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IV. CONCLUSION 
Values of the radiation force P computed from Equation (29) are shown in 
Fig. 2 as a function of the distance of the spherical obstacle from the line 
Fic. 2. 
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source, origin of the diverging sound field. Curves are given for several sizes 





of spherical obstacle, a = ka; both coordinates contain @ as a normalizing 
factor in order that all the curves may be plotted more compactly on the same 
graph. In order that both positive and negative forces, attractions and repul- 


pressure amplitude of the sound field at the same distance ro from the line 


(30). The factor sr/k may be replaced by |p|/kcpo|H'‘? (kro)| where |p| is the 
source as the spherical obstacle. 


The coefficients may then be taken in pairs as necessary to evaluate Equation 
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sions respectively, may be shown on the same logarithmic scale a constant 
has been added to the ordinate so that all the values to be plotted are positive. 

For all finite sizes of sphere the force at an infinite distance from the line 
source is one of repulsion, as expected since the wave fronts then tend to 
become plane. As the distance between the line source and the spherical 
obstacle is decreased so the radiation force of repulsion decreases until at some 
distance, which is a function of the size of the sphere and the frequency of the 
sound field, it becomes zero. At still smaller distances the radiation force 
becomes an ever-increasing attraction. As ro/a approaches unity the force of 
attraction becomes infinitely great; this corresponds to the line source being 
tangential to the spherical obstacle. 

In an earlier paper (Embleton 1954) curves are given for the radiation 
force acting on a spherical obstacle placed in a diverging spherical sound field 
as a function of size of obstacle and its distance from the center of the sound 
field. Qualitatively these are similar to those shown in Fig. 2. However in the 
case of a spherical sound field the distance from the center of the sound field 
to which the force of attraction extends is greater than for a cylindrical field 
providing the size of sphere and frequency of the sound field are the same in 
both cases. Also for a given size of sphere situated at a constant distance from 
the source of the sound field any force of attraction decreases as the frequency 
of the sound field increases and changes to a repulsion at a frequency which 
depends upon both rp and a—this frequency is lower in a cylindrical wave 
field than in a spherical field (for the same rp and a). 
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APPENDIX 

The following integrals are used in the evaluation of the expressions (19), 
(23), and (27). In all cases only the values for N > n were required. 
m/2 for M=m+0 


x for M=m=0 
0 for m ~ M 


(1) f cos m2... cos MQ . dQ 
0 


_ ____ 2(stmt+i)! as le 

3 ~ (n—m)! (2n+1)(2n+3) tor N= a+) 
(2) fu. PRu). PR) dy 
, 0 for N=n or N> 2n+1 
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- (n?—m’+2n) . (nt+m+1)! 
~ (n—m)! (2n+1)(2n+3) 


: 0 0 m 2 
(4) : 2 pay) 2 PR). =u!) nd 
= 0 for N=norN>n+I1 


REFERENCES 


AwartTAnl, J. 1955. J. Acoust. Soc. Amer. 27: 282. 

EMBLETON, T. F.W. 1954. J. Acoust. Soc. Amer. 26: 40. 

Kinc, L. V. 1934. Proc. Roy. Soc. (London), A, 147: 212. 

Lams, H. 1932. Hydrodynamics. 6th ed. Cambridge University Press, London. 

Watson, G. N. 1944. Theory of Bessel functions. 2nd ed. Cambridge University Press, 
London. §13.47. 


for N = n+1 








EVAPORATION NEUTRONS FROM COSMIC RAY NUCLEAR 
INTERACTIONS IN VARIOUS ELEMENTS! 


By K. W. GEIGER 


ABSTRACT 


Production of evaporation neutrons from nuclear interactions by cosmic rays is 
investigated in C, Al, Fe, Cu, Pb, and Bi. Boron trifluoride counters in a paraffin 
moderator in which the element under study is placed serve to detect the evapora- 
tion neutrons. A cosmic ray telescope allows the consideration of ionizing 
neutron-producing particles separately and it is found that of the total neutrons 
14.6% are produced by ionizing particles. The mean neutron multiplicities are 
worked out and they show an increase with the two-thirds power of the atomic 
weight of the element in which the interactions take place. Also, the rate of 
interactions is calculated and this gives the relative interaction cross-sections 
which are found to increase slightly more with atomic weight than expected if the 
cross-section is geometric. An estimate of the intensity of the N-component is made 
and discussed in connection with measurements published by other authors. 


I. INTRODUCTION 


The component of cosmic radiation capable of producing nuclear interactions, 
the N-component, is known to consist mostly of fast neutrons and protons. 
At sea level only a small contribution is due to x-mesons, captured u~-mesons, 
fast 4-mesons, and heavy mesons. For the study of the nuclear interactions 
several methods are available. For example, by the use of nuclear emulsion or 
cloud chamber the interactions can be observed directly as ‘‘stars’” when the 
reaction products are ionizing particles. Another, more indirect, way is the 
use of neutron counters in a paraffin moderator where the interactions are 
studied by detection of evaporation neutrons. In particular, the multiplicity of 
neutrons can be evaluated. The advantage of this method is that it is based 
on a simple counting procedure and that the material in which the neutron 
production takes place (hereafter referred to as the producer) is easily inter- 
changeable. The disadvantage is the relatively low neutron counting efficiency 
in such a recording arrangement which makes the results more difficult to 
interpret. 

Neutron production in lead was first studied by Cocconi Tongiorgi (1949) 
by whom it was established that nuclear interactions were responsible. A 
very thorough investigation was made by Cocconi, Cocconi Tongiorgi, and 
Widgoff (1950) where the production and multiplicity of neutrons from 
energetic nuclear reactions, capable of emitting at least one penetrating 
ionizing particle, was studied in three different thicknesses of lead and in 
aluminum. The authors concluded that owing to the interaction of secondaries 
a cascade process develops because the average multiplicity of neutrons 
originating from a lead producer increases with increasing producer thickness. 
Therefore, it was estimated that the multiplicity for single interactions could 
only be evaluated when the thickness of the producer is not more than 1/10 of 
the interaction mean free path of the N-component. In a study of nuclear 
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disintegrations by u-mesons Cocconi and Cocconi Tongiorgi (1951) obtained 
an average neutron multiplicity in lead of 4.5 due to the entire N-component. 
Ortel (1954), using a liquid scintillator as well as neutron counters to detect 
the nuclear interactions directly, found an average neutron multiplicity of 1.3 
for carbon. Absolute measurements of the neutron production at sea level were 
performed by Tobey and Montgomery (1951). 

Thus far no systematic study had been made on the variation of multi- 
plicities with different elements in which the disintegrations take place. 
This, therefore, was the main purpose of the present investigation. In addition, 
the relative production rates and hence the variation of interaction cross 
sections with atomic number are evaluated. The total N-component is con- 
sidered as well as the ionizing part only. The experiment was performed during 
the period from February to August 1954 in Ottawa, 100 meters above sea 
level. 


II. APPARATUS 


The arrangement of neutron counters and of the producer in the paraffin 
moderator is shown in Fig. 1. The two neutron counters were 6.4 cm. in 
diameter and had a sensitive length of about 50 cm. They were filled with 


)GM COUNTER “ 
TELESCOPE vs 


NEUTRON 
COUNTERS 
PARAFFIN 


_/ MODERATOR 





SCALE o 10 20 3x0 


Fic. 1, Arrangement of producer, neutron counters, and cosmic ray telescope in the parafhin 
moderator. 


boron trifluoride (enriched in B'’) at a pressure of 60 cm. Hg and the maximum 
energy of the recorded neutrons was estimated to be about 10 Mev. Therefore, 
practically, only evaporation neutrons are detected. These are emitted 
isotropically, as was shown in a similar type of experiment by Coccont e a. 
(1950). It is known from nuclear emulsion work that evaporation protons are 
also emitted isotropically (e.g. Harding et al. 1949). 

The producer area was 25 by 38 cm. and when the material did not take up 
the full available height of 8.8 cm. the producer was made up in the torm of 
sheets and spaced evenly to avoid corrections due to a changing geometry 
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A wide angle cosmic ray telescope placed on top of the producer, with the 
area of the lower tray slightly greater than the producer area, served as a 
means for distinguishing between the total and the ionizing neutron-producing 
component. 

The neutron detection efficiency was determined by placing a calibrated 
0.1 mc. Ra+Be neutron source at different locations within the producer so 


that the average efficiency 
1 
c= TS cenededyds 


could be worked out. The energy spectra of the neutrons from a Ra+ Be source 
and from the evaporation processes are similar enough to justify this method 
(Cocconi et al. 1950). For 8.3 cm. of lead in the producer space an integrated 
efficiency of « = 1.7% has been found, for all other producer materials and 
thicknesses an efficiency of 1.5% was determined. The error was estimated 
to be +20%. The neutron source was only known to +13%, the remainder of 
the error being due to experimental uncertainties. 

The neutrons which come from a nuclear interaction in the producer are 
slowed down in the paraffin and have a mean life of about 150 usec. before 
they are captured in the counters or in the paraffin. Hence neutrons which 
have a common origin do not reach the counters simultaneously and groups of 
neutrons arising from one interaction are spread over a time interval compar- 
able with their mean life. The neutrons are detected successively and groups 
can be resolved. Because of the low efficiency of the arrangement generally 
only a very few neutrons from one event are detected, but a statistical 
analysis described later makes it possible to work out the multiplicity. 

The electronic circuit is shown in Fig. 2. It can be divided into two parts; 
one records the neutrons produced by the ionizing and the other by the total 
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Fic. 2. Block diagram of electronic circuit. 
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N-component. The amplified neutron pulse after passing through a pulse 
height discriminator is fed into two gated amplifiers which are normally 
closed. A coincidence pulse from the telescope after a delay of 6 usec. triggers 
a gate pulse, 320 usec. long, about twice the half-life of neutrons in the paraffin 
moderator. The amplifier gate is opened for this time, one or more neutron 
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pulses can pass, and events are selected according to the number of pulses 
within the gate time. The gate of the second amplifier is opened by the first 
neutron pulse of any event and all events of more than one neutron detected 
are analyzed in the second selector. The number of events with one neutron 
detected can be calculated by subtracting the number of multiple events 
from the total number of gates which is recorded separately. 

All information as well as rates from some auxiliary circuits was fed into 
message registers. A 20-pen operation recorder was also connected, and proper 
functioning of the equipment could be verified by analyzing the chart. In 
addition, the selectors were checked periodically by applying a square wave 
of known frequency admitted only during the gate opening. This should record 
an event with the corresponding number of square pulses. Further details on 
the circuit are given in a previous paper (Geiger and Rose 1954) where similar 
equipment was used to study the multiplicity of neutron production from the 
spontaneous fission of uranium. 

The limited gate length reduces the neutron counting efficiency. With a 
mean life of thermal neutrons 7 the gate efficiency is 


t2 
: -if e "dt 
TH 


where ¢, and ¢, represent the opening and closing time of the gate after the 
triggering event has taken place. The neutron mean life has been found 
experimentally to be 154+13 usec. by using three different gate lengths of 
70, 120, and 320 usec. and comparing the corresponding counting rates. The 
gate was triggered by an ionizing particle and lead was used as producer. 

The slowing-down time of fast neutrons to thermal energies is of the order 
of the delay until the gate opens, i.e. 6 usec. (Glasstone and Edlund 1952). 
Hence, slightly different gate efficiencies are to be expected, depending on 
whether the gate is opened by a first neutron pulse or a telescope coincidence 
pulse. In the first case, which applies to the part of the circuit recording total 
neutron production, ¢; and f, are 6 and 326 usec. respectively and a gate 
efficiency of f = 0.84 was calculated. In the second case, that of the ionizing 
producing component, the time before the gate opens is used up to slow the 
neutron down. Therefore, ¢; and fz, are 0 and 320 usec. and the gate efficiency 
is f’ = 0.87. 

Ill. DATA AND CORRECTIONS 


Only neutrons which originate in the producer are of interest and, hence, 
account has to be taken of neutrons from the paraffin wax and the surrounding 
material. The counting rate due to the producer alone is 

R= R,—Ry+ (AR,+ AR,y) 


where R, is the rate with the producer inserted, 
R, is the background rate without producer, 
and the term in brackets takes into account the effect of the presence of the 
producer on the background rate. 
R, is composed of three parts, namely 
Ry = Reat Rut Ra. 
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Rea is the rate when the neutron counters are surrounded by cadmium; 23+3 
counts per hour have been found, and this was, as expected, practically 
constant whether a producer was present or not. R,, is the rate due to neutrons 
produced outside the arrangement and diffusing through the wax to the count- 
ers, and R, is that due to neutrons produced by the N-component in the 
paraffin wax. Absorption or scattering is expected to decrease these contri- 
butions when the producer is present. 

The influence of the producer on R,, is difficult to evaluate. An auxiliary 
experiment has been performed where a 1 mc. Ra+Be neutron source was 
placed on the top center of the pile to investigate scattering from the producer. 
It must, of course, be assumed that the neutron energy spectra from the source 
and from the surroundings are similar after they have been filtered through 
12 cm. of paraffin, the amount on top of the producer space. The observed 
reduction due to producer scattering is shown in Table I. 


TABLE 1 
SCATTERING OF NEUTRONS ON PRODUCER 


Amount scattered 


Producer s 
Cc 8.5cm. 0.37 
Al 8.3cm. 0.29 
Fe 3.8cm. 0.34 
Cu 3.8cm. 0.37 
Sn 3.7cm. 0.19 
Pb 1.9cm. 0.11 
Pb 8.3cm. 0.33 
Paraffin 12.7cm. 0.75 


An idea of what fraction of R,, is due to diffusing neutrons passing through 
the producer space was obtained by determining the reduction in counting 
rate when the space is filled with paraffin. A decrease from 335 to 304 counts 
per hour was found. The minimum value for neutrons passing through the 
producer space, therefore, is 31 counts per hour. As the last line in Table I 
shows, this only represents 75% of the total, and furthermore, the figure is 
still too low because of neutron production in the carbon of the added paraffin 
which makes the above-mentioned decrease greater than measured. From 
these considerations 50+10 counts per hour were estimated. The decrease of 
the background rate due to insertion of the producer is therefore 


AR. = 50s (counts per hour) 


with s taken from Table I. This correction only applies to the rate of events 
with single neutrons, produced by the total N-component. 

The remaining contribution to the background is R,. Most of R, would be 
due to production in the vicinity of the neutron counters where the detection 
efficiency is largest. Therefore, from geometrical considerations as well as by 
taking into account the other contributions to R,, it was estimated that 40% 
of R, or 0.4335 = 134 counts per hour are due to the part of the N-compo- 
nent traversing the producer space. Hence, with the producer present: 


ARy = 134(1—e-?e) (counts per hour) 
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where D is the thickness of the producer and \, the absorption mean free path 
of the N-component. 

Since the ionizing N-component is confined to the producer space by the 
cosmic ray telescope it was estimated that 90% of Ry,; (production by ionizing 
particles) is affected by the presence of the producer: 

AR,,; = 7.9(1—e-?2) (counts per hour). 


Table II shows the observed hourly rates of events consisting of a group of m 
neutrons and also the background corrections. It can be seen that the influence 
of the producer on the background is small for producers of high atomic weight 
but is appreciable for carbon and aluminum. 

R, further includes small corrections for neutron pulses which fall into the 
gate accidentally. Barometric corrections were also included which, however, 
were small since the rates are averages over about seven days distributed 
over the whole time during which the experiment was performed. This brought 
the average pressure for those days near to the mean for the whole time. 

The quoted errors are statistical, with a minimum error of +1% to take 
account of variations of the intensity of the N-component. During the months 
this investigation was performed solar activity was low. Time variations of 
the N-component, as observed by Rose (private communication) at the same 
location, were small enough to justify such a small error. 

IV. DIRECTIONAL INTENSITY OF IONIZING NEUTRON-PRODUCING 
RADIATION 

When the total nuclear component is studied the producing particles can 
enter the producer from all directions. However, for particles which have 
to go through the cosmic ray telescope the angle is limited and the zenith 
angle distribution has to be known to apply the necessary correction. A 
simple method for obtaining the directional intensity is to take the counting 
rate as a function of the height of the top telescope tray. The cosmic ray 
intensity from the zenith angle 6 is assumed to be J = J, cos"6, using J, for the 
vertical intensity. The total flux per unit sphere is then 


= /2 
(1) Jo= onl, [ cos"6 sin 6d0 = 2rI,(n+1)". 
0 
The vertical intensity can be determined: 
d° 
ee a 


where d is the distance between the two telescope trays, A, and A, their areas, 
and & the observed neutron counting rate. In the present case, a lead producer, 
5.7 cm. thick, has been used and the rate without producer subtracted. A 
lower limit for the total flux was determined by increasing the top tray area 
by a factor of two. This was done by moving the top tray to one side so that it 
covered half of the producer and multiplying the coincidence rate by two. 
The use of the lower tray only to trigger the gate is not suitable because of 
the possibility that fast neutrons change into protons within the paraffin 
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between the trays. After corrections for accidentals and barometric effect the 
following intensities were obtained: 


I, = 0.070 hr.—! cm.—? sterad.—!, Jo = 0.069 hr.— cm.—?. 


These values do not represent the actual intensities of the ionizing part of the 
N-component because of multiple neutron production and the limited neutron 
counting efficiency. An exponent of » = 5.4+0.6 follows by inserting these 
values into formula (1). 

A different approach, avoiding the use of the total flux, can be chosen. 
The integral in (1) extends over one half solid angle. When the rate is taken 
with the trays positioned as in Fig. 1, m can be found when J, and the solid 
angle over which to integrate are known: 


6’ 
(2) Jo: ™ orl, { cos"@ sin 6 dé 
0 


where Jo,, represents the neutron rate in that geometry. 6’ could be found from 
an auxiliary experiment on the total ionizing cosmic ray component capable 
of penetrating the paraffin within the telescope. For this component the expo- 
nent is known to be near = 2 (Voisin 1950) and also has been determined 
experimentally using the method described in the previous paragraph to be 
n = 2.1+0.2. To obtain 6’ the lower tray was placed in the center height of 
the producer space to simulate the angle in which the ionizing neutron- 
producing particles are admitted. 6’ = 40° followed using a formula corre- 
sponding to (2) but for the intensity per unit area. Therefore 


J 6’ = 40° 1 
(3) =*t= f cos"@ sin 6 d@ = f cos" d(cos @) 
2nI, 0 cos6’ = 0,766 


ae 1 _n7 an+1 
ae 0.766"""). 


With Jo,, = 0.057 hr.—' cm.~? the same result as before, namely » = 5.3+0.6, 
has been found. 

The result can be compared with a value for protons of » ~ 5, given by 
Mylroi and Wilson (1951), which applies to protons having an energy of 
from 400 Mev. upwards. It is believed that most of the neutrons in the present 
experiment are produced by protons within that energy range. The exponent 
is energy dependent, being high, approximately 7, for very fast protons for 
which angular deviations after interactions are small, and m is based purely 
on the collimation due to absorption along different proton path lengths 
through the atmosphere. 

To make the neutron production rate due to the ionizing component 
equivalent to that from the total, the former has to be multiplied by a 


factor 
Jo/Jo, = 1.21 


to take account of the limited telescope angle. 
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V. ANALYSIS AND RESULTS 


As pointed out under II the number of neutrons m in a recorded event, 
or recorded multiplicity, is not the actual neutron multiplicity from nuclear 
interaction, owing to the low neutron detection efficiency. Two slightly 
different approaches for the analysis have to be made according as the gate is 
opened by the first neutron pulse (total producing component) or by a tele- 
scope coincidence pulse (ionizing producing component). 

(a) Total producing component.—The probability that k of v neutrons 
released from one interaction are detected in the counters follows a binomial 
distribution 


(4) P(k, v) = (:) e‘(1—e)’™, v> k, 


with ¢ being the neutron detection efficiency. Similarly, if a total of k neutrons 
is detected and the gate is opened by the first neutron, the probability that m 
of these & neutrons fall within the gate is 


(5) Pim, #) = (#71 tapes, km 


where f is the gate efficiency. The aim is to work out the average multiplicity 
¢ from the nuclear interactions. To simplify the analysis an assumption for 
the shape of the multiplicity spectrum has to be made. Following Cocconi 
et al. (1950, 1951) an exponential multiplicity spectrum has been used. The 
justification for its use in the present case will be seen later. Thus, 


Ie) =. —1)4 shee S16) = 4. 


The average multiplicity is 
= ae 1 1 
b= Dv(e—-1lhe"” =——y,  p=-+05, a K1. 
I l—e a 
The observed rate of events with m neutrons becomes therefore 


Rn = OD ¥ I) Plk, ») Plm, k) 


vem k=m 
with Q being the rate of the interactions. This summation can be evaluated: 
Rn = O-— C15" 7 = 
m i-—"n-G—<g oT 
To evaluate a and therefore 7 it is convenient to use the logarithmic form of 
this equation and to differentiate it with respect to m 
d log Rm 


= log ef —a log e—log{1— (1—f )e™. 


(6) 


This shows that &, versus m should be a straight line when plotted on semi- 
logarithmic paper, the slope being a function of the mean neutron multiplicity. 








| 
| 
| 
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(6) In the case of a charged particle opening the gate the limited gate length 
decreases the neutron detection efficiency and in (4) k can be replaced by m 
and ¢ by ¢f’. Equation (5) need not be used. Therefore 


Bh, w OT 1) Pas), 


v=m 


which after summation gives 


e*—1)(f’)"e™ 


pe 
Ka = Cheb 





or 


(7) 2108 Rn — jog ef’ —a log e—log{l— (1- fe. 
This also represents a straight line on a semilogarithmic plot. The slope is 
equal to that given in (6) apart from the differences between f and f’. 

Further, it is of interest to know the actual rate of interactions taking place 


HOURLY RATE 


6-3 Cm Pb. 


1OMIZING 





RECORDED MULTIPLICITY m 


Fic. 3. Data on lead in semilogarithmic plot. The first two points of the curves have been 
used to work out the multiplicities. 
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in the producer. The rate can be taken from the formulae above or, more simply, 
is given by 
Q = N/ve 


where N is the total neutron rate. N can be taken from the data by adding up 
the rates multiplied by m and correcting for the gate efficiency. In the case of 
neutrons produced by the ionizing cosmic ray component WN has been multiplied 
by the factor 1.21 to take account of the limited telescope angle as described 
in Section IV. 

Fig. 3 represents the experimental data of Table II for 1.9 and 8.3 cm. 
thickness of lead on a semilogarithmic plot. If the assumed exponential 
multiplicity spectrum is correct straight lines should appear on the graph. 
This is not the case for 8.3 cm. of lead where the curve flattens out, the expected 
effect due to cascades of nuclear interactions as concluded by Cocconi et al. 
(1950). A reasonably straight line is obtained with 1.9 cm. of lead. The 
remaining slight curvature is assumed to be due to non-uniformity of neutron 
detection efficiency over the whole of the producer. The number of cascades 
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must be very small, the lead thickness being one-eighth of the collision mean 
free path of about 160 gm./cm.? Also, in contrast to the experiment of Cocconi 
et al. (1950) which discriminated against low energy events, all neutron- 
producing nuclear interactions are included, the low energy ones, not capable 
of producing cascades, being the most frequent. For the producer materials as 
listed in Table II, in order to eliminate cascades, reasonably small thicknesses 
were always used, thereby sacrificing some statistical accuracy due to lower 
neutron production. To calculate the mean multiplicities the slope through 
the two points R; and R, has been used which is thought to give the most 
realistic value for the multiplicities. 

Fig. 4 shows two graphs in double-logarithmic form representing the mean 
multiplicities for total and for ionizing production as a function of the mass 
number of the producer. For lead, here as well as in the following graphs, 
only the values from 1.9 cm. thickness are included. The dotted lines have, 
tentatively, been drawn according to an A?/*-dependence, and this fits the 
experimental points quite well. 

The errors for > as given in Fig. 4 include only the statistical error from the 
measured rates and not the uncertainty for the neutron detection efficiency, 
the purpose of this work being not so much to obtain absolute figures as their 
comparison for a variety of elements. When account is taken of the additional 
error the mean multiplicity for neutrons produced by the total N-component 
in this lead (¥ = 6.5) compares reasonably well with the value measured by 
Cocconi and Cocconi Tongiorgi (1951) of 4.5 or with Ortel’s (1954) value of 6. 
Ortel also gave a value for carbon of 1.3 which is in good agreement with the 
present value of 1.4. 

As expected, the average multiplicity was found to be higher for neutrons 
produced by the ionizing component than for neutrons produced by the neutral 
part of the N-component (or the total N-component), the figures for 1.9 cm. 
of lead being 9.3 (ionizing) and 6.5 (total). This difference is well known from 
nuclear emulsion work (Brown et al. 1949), the stars with only a few prongs 
being mostly produced by neutral particles. At higher particle energy (many- 
prong stars) the relative number of stars generated by ionizing and non- 
ionizing interacting particles is about equal. Similarly in the present experi- 
ment, the rates for the neutron groups R;, R2, etc. approach a ratio of 1: 2, 
for ionizing to total production, as the number of neutrons in the group be- 
comes large (see Table I1). The explanation for this is as follows: The nuclear 
interactions are mostly caused by fast neutrons and protons which are pro- 
duced locally. The interaction mean free path for these particles in air is close 
to the geometric mean free path of 68 gm./cm.? It can be assumed that the 
number of fast protons and neutrons produced in the air is equal. The protons 
undergo ionization losses which rise quickly when the proton energy falls 
below 500 Mev., the range then being comparable to the mean free path for 
nuclear interactions. Only the energetic protons survive and are capable of 
producing further nuclear interactions. The average energy of the protons, 
therefore, is higher than that of the fast neutrons. 

For the same reason the number of neutrons produced by the total N- 
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component must be more than twice as high as the number of neutrons 
produced by the ionizing N-component only. This is shown in Fig. 5 which 
represents the rate for the total neutron production for both cases as a function 
of the atomic weight of the producer. Taking the mean of the data on neutron 
production for the heavier elements (from Sn up) where the accuracy is not 
too low, 14.6% of the total neutron production is due to the ionizing compo- 
nent. This compares well with 13%, a figure which can be taken from data as 
given by Hogrebe (1952) for lead. The figure quoted in that paper (18.5%) 
is not suitable for comparison because it is based on production at infinite 
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Fic. 5. Total number of neutrons versus atomic weight. 

Fic. 6. Number of interactions per gram atom per hour versus atomic weight. This is 
proportional to the neutron production cross-section. The slope of the dotted line is taken 
from Fox et al. (1950). 


lead thickness, there being appreciably different absorption mean free paths 
for total and ionizing N-component. For very energetic interactions Cocconi 
et al. (1950) have shown that the neutron production is about equally divided 
between neutral and ionizing particles. This is in agreement with the above 
considerations. 

It is of interest to show the relative rates for the two types of interactions Q. 
Assuming the same interaction cross-sections for the neutral and ionizing 
N-components the ratio of their intensities must be equal to the ratio of the 
corresponding interactions. In Fig. 6 the rates of interactions per gram atom 
are drawn against atomic weight of the producer and 10.5% of the interactions 
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have been found to be due to the ionizing N-component. The interaction 
rate per gram atom is proportional to the interaction cross-section of the 
nuclei and hence Fig. 6 gives also the variation of cross-section with mass 
number. A geometric cross-section would require an A?/-dependence. Actually, 
the dashed lines have been drawn following the variation of total neutron 
cross sections as measured by Fox et al. (1950) for artificially produced neutrons 
of 280 Mev. They are somewhat steeper than A?/* which might be expected 
owing to transparency of nuclei of low atomic weight. 


VI. SUMMARY AND DISCUSSION 


The production of evaporation neutrons by cosmic rays in different materials 
has been studied and a distinction has been made between the ionizing and 
the total N-component. It has been shown that fast neutrons are the most 
abundant part of the N-component confirming observations with photographic 
plates. 

The mean multiplicities of the evaporation neutrons and the rates of nuclear 
interactions were evaluated on the basis of an exponential neutron multiplicity 
spectrum. The mean multiplicities have been found to increase as the two- 
thirds power of the mass number of the producer material. For a theoretical 
consideration the conditions in the present experiment are rather complicated 
because of the fact that the incident neutron-producing particles are not 
monoenergetic. Jackson (private communication) computed the mean number 
of evaporation neutrons versus incident nucleon energy for mass numbers > 100 
using excitation energies (McManus, Sharp, and Gellman 1954) obtained with 
the Monte Carlo method. For energies of the incident particles between 90 
and 160 Mev. these calculated values follow closely an A?/*-law and are also 
in satisfactory agreement with the absolute values as given in Fig. 4 for the 
production by the total N-component. The average energy for incident 
neutrons is of this order as it can be taken from a graph given by Rossi (1952) 
for the computed energy spectrum of neutrons in air and considering only 
neutrons above 30 Mev. Not many interactions with neutron production 
could be expected below that energy. 

The number of interactions per gram atom of producer material, which 
is proportional to the interaction cross-section, has been found to increase 
somewhat more steeply with the mass number of the producer material than 
expected if the cross-section is geometric. This indicates a certain transparency 
for lighter nuclei. 

Using the data from the 1.9 cm. lead producer the intensity of the part of 
the N-component capable of producing nuclear interactions with emission of 
evaporation neutrons can be estimated. The omnidirectional intensity is: 


-2t1Q 

Jo= 42 

The factor (n+1): (n+2) converts the intensity per unit sphere to intensity 
per unit area. The value of m is not known but has been assumed to be 5. 


Aime is the interaction mean free path for the N-component in lead. Lacking 
other information Aj, has been assumed to be geometric (160 gm./cm.’). 


d int- 
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Q/W is the number of interactions per gram of lead and is 0.24X10-* gm.-! 
sec.—!, Therefore, the omnidirectional intensity becomes 


Jo = 33X10~‘ sec.-! cm.-? 


with an estimated error of +30%. This represents the intensity at sea level 
since the roof and the paraffin on top of the producer consist of about 15 gm./ 
cm.? of light material, and cancels with the difference in atmospheric depth to 
sea level (13 gm./cm.?). The figure agrees with Jy = 40X10~‘ sec.—! cm.~? as 
given by Cocconi and Cocconi Tongiorgi (1951) who use the same value for 
the interaction mean free path. As pointed out in their paper this value is also 
consistent with estimates for the intensity of the N-component as deduced 
from star production in nuclear emulsions and from measurements of neutron 
production in the free atmosphere. 

A similar calculation can be applied to the ionizing part of the N-component. 
Not all of the interactions are caused by protons. The main additional contri- 
bution is due to nuclear absorption of u~-mesons which is followed by neutron 
emission. The number of stopped y--mesons in the 1.9 cm. lead producer 
can be worked out. According to the differential range spectrum of mesons at 
sea level as given by Rossi (1948) the number of u~-mesons stopped in air is 
I, u- = 2X10-* gm.—' sec.—! sterad.—!. Then, the number of u--mesons stopped 
in the lead producer and passing through the cosmic ray telescope, the effect 
of which is given by equation (3), is 


2a n+1 Ee 
Jou- = n41 01 0-766 ) : Ww 

where p = 1.7 is the equivalent absorber thickness of lead relative to air, W 
the weight of the producer (20.9 kgm.), and m = 3 has been assumed for slow 
mesons (Voisin 1951). Practically all of these stopped u--mesons are captured 
in a lead nucleus followed on the average by the emission of 2.14 neutrons 
(Widgoff 1953). Assuming again an exponential multiplicity spectrum, this 
would give a rate R; = 2.6 and R2 = 0.04 which would change the values 
for 1.9 cm. of lead to Ri = 11.0 and R2 = 1.28 in the row for production by 
ionizing particles in Table II. This brings the average multiplicity up to 
¢ = il and the interaction rate down from Q = 166 to Q = 116 per hour. 

The number of interactions due to fast u-mesons is very small. A rate of 
Q = 7 hr.~! is due to them, using an interaction cross section of 1X10-?° 
cm.?/nucleon and the above multiplicity. 

Neglecting interactions due to 7-mesons the number of which in air at sea 
level is very small and neglecting also interactions due to photons accompanied 
by an electron which triggers the cosmic ray telescope, the rate of interactions 
which are thought to be initiated by protons is Q = 109 hr.-!. With Aj, = 160 
gm./cm.?and » = 5 the vertical proton intensity is 


n+1 = a 2 wei 
I, = — ¢ Nine = 2.3107‘ sec.~' cm.” sterad.~. 
This figure applies to protons with a lower momentum limit of 550 Mev./c 


which is given by absorption in the paraffin on top of the producer and in the 
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roof. The same considerations as in the case of the total N-component apply 
making the above figure the sea level intensity. 

The figure can be compared with the vertical intensity as derived from 
experimental determinations of the differential proton spectrum. Up to 
103 Mev./c Ogilvie (1955) draws a spectrum through his own measured 
points and gives those of previous authors. From 550 Mev./c to 10? Mev./c 
the value J,,, = 0.65 X 10~‘ sec.—' cm.~? sterad.—! can be taken from his graph. 
Mylroi and Wilson (1951) give the intensity for momenta above 10’ Mev./c 
with J,,, = 0.55 X10~4 sec.—! cm.~? sterad.—'. Therefore, at sea level, the verti- 
cal intensity of protons above the energy limit set by the absorber on top of 
the producer is J,,, = 1.20X10~‘ sec.-! cm.~? sterad.—!. 

Although good agreement with the figure derived from the present experi- 
ment can not be claimed, it must be considered that experimental errors are 
quite high. It is felt that because of the various, not too well known, factors 
entering the calculation, the present determination can be out by as much as 
+50%. Also, there is a considerable spread in the intensity determinations 
by different authors as shown by Ogilvie in his graph. 

It is generally assumed that nearly all nuclear interactions which are 
produced by non-iortizing particles are due to neutrons. Therefore, Q(total) 
represents reasonably well the interaction rate when Q(ion.) is subtracted. 
But Q(ion.) is not entirely due to protons and an appreciable proportion is 
due to other ionizing particles of which the N-component is composed. For 
heavier elements, where nuclear absorption of u--mesons is important, this 
is the main additional contribution. One therefore has to bear in mind that 
the results which are given in Section V for the ionizing N-component do 
not apply to the cosmic ray protons only but to a mixture of various ionizing 
particles. 
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THE BEHAVIOR OF FERROMAGNETICS UNDER 
STRONG COMPRESSION! 


By F. D. Stacey 


ABSTRACT 


Specimens in the form of thin disks have been subjected to strong compression 
simultaneously with high magnetic fields by means of a new design of perme- 
ameter, with which their magnetizations have been measured. Magnetization 
curves extending to saturation have been plotted for a number of materials at 
(non-hydrostatic) pressures up to 10,000 atmospheres. The most interesting 
result is the strong increase in saturation magnetization shown by nickel and the 


nickel-copper alloys. 
1. INTRODUCTION 


The application of high magnetic fields to materials under hydrostatic 
pressure presents technical difficulties which have led to several compromises. 
The measurements of Yeh (1925) and Steinberger (1933) were not extended 
above 100 oersteds, while the method used by Ebert and Kussmann (1937) 
can be shown to be invalid because they measured the changes in the total 
flux caused by interchanging specimens in the field gap*of an electromagnet 
with saturated pole-pieces. 

Another approach was suggested by Bridgman’s (1950) method of applying 
high pressures to foils compressed between massive anvils. This does not 
produce hydrostatic pressures, but the present work has shown that there is 
an important hydrostatic component to the compression of soft specimens. 

The surprising increase in the saturation magnetization of nickel, com- 
pressed in this manner, has been reported previously by Jones and Stacey 
(1953). Measurements have now been made on the ferromagnetic nickel- 
copper alloys and a number of other metals, showing that this effect is a 
characteristic both of pure nickel and of nickel-copper alloys. 


2. EXPERIMENTATION 


The permeameter was a low reluctance magnetic circuit, machined from 
black mild steel, and shown assembled in Fig. 1. Specimens were disks of foil 
inserted between the pole-pieces through one of two circular openings in the 
outer case, and pressure was applied to them by placing the whole system 
between the anvils of a hydraulic press. The compression of the specimens 
was estimated from the oil pressure in the press, assuming that the total force 
on the piston was transferred to the specimens. A correction for friction of 
the piston was avoided by allowing a small leak of oil past it, the oil pressure 
being maintained by a backing of nitrogen. 

Magnetizing direct currents between 2 and 60 ma. were passed through 
the primary windings, totalling 3000 turns, in series with the primary coil of a 
standard variable mutual inductance. The mutual inductances M, and M;, 


1Manuscript received November 4, 1955. 
Contribution from the Department of Physics, University of British Columbia, Vancouver, 
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coiL 


INNER SECONDARY 
COIL 


between the primary and each of the secondary coils of the permeameter were 
measured by the automatic reversal of the primary current with each of the 
secondary coils connected in series opposition to the secondary of the induc- 
tometer and to a ballistic galvanometer of fairly long time constant. These 
measurements were used to calculate the flux through the specimens and 
through the whole permeameter. The magnetizing current, 7, was measured 
by balancing the voltage across a resistance box in the primary circuit against 
a standard cell. 

Each secondary coil on the permeameter had 10 turns, giving a low re- 
sistance to the secondary circuit and strong damping to the galvanometer. 
This facilitated the direct measurement of the impure (iron-cored) inductances 
between the permeameter coils by balancing them against the pure (air-cored) 
inductance of the inductometer. At each reversal of the primary current the 
galvanometer first showed a sharp deflection which was reversed and returned 
to a sudden stop, the final deflection indicating the extent to which the bridge 
was off balance. 

Complete reproducibility was obtained by beginning measurements on 
each specimen with the highest pressure, and, at each pressure, subjecting 
the system to a cycling demagnetization process before taking readings. 


3. ANALYSIS OF READINGS 
The magnetic intensity in a specimen is given by 


(1) B = 1M ,/nza, 


where i is the primary current, M, is the mutual inductance between the 
primary coil and the inner secondary of mz turns, and a is the specimen area. 
There is a small correction for the flux outside the specimen but enclosed by 
the inner secondary coil. 

A consideration of the reluctances of the different parts of the magnetic 
circuit gives an expression for the relation between the flux through the 
specimen, which is measured by M,, and the flux through the whole per- 
meameter, measured by M2, the mutual inductance between the primary and 
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the outer secondary. It is convenient to make the specimen reluctance the 
subject of the resulting equation: 


(2) eae (M2/M,)—(M2/M1)o 
Ha 1—a[(M2/M,)— (M2/Mi)o] 

t and yw are the thickness and permeability of the specimen, (M2/M;)o is the 
value of the ratio M2/M, observed with no specimen (at the same value of 
(1M,) or flux through the pole-piece tips), and A is a constant. a is a small 
correcting factor which takes account of the air flux measured with M, and 
of the fact that the pole-piece reluctances vary slightly with specimen re- 
luctance. Since A is a calibration constant of the apparatus a is conveniently 
included in it. 

The constants A, a, and (M2/My)o of equation 2 are characteristics of the 
permeameter, being functions of both pressure and flux. Calibration tests 
were made with standard thicknesses of copper so that graphs of these con- 
stants against (4M) were plotted for each of nine pressures. Tests at constant 
magnetic flux and pressure with seven thicknesses of copper gave exact fitting 
to equation 2. 

Equations 1 and 2 are sufficient to determine the intensity of magnetization, 
I, and field, H, in a specimen, since: 


(3) 


(4) 


Specimen thicknesses were between 0.002 cm. and 0.02 cm. and were meas- 
ured with a micrometer, with an accuracy very much less than that of the 
other measured quantities 1, M,, and M2, which were determined to one part 
in 104 or better. A method was devised of correcting for the errors introduced 
to the results by inaccuracies in the measurements of specimen thicknesses, so 
that the micrometer values were needed only as a first approximation. 

For the non-magnetic specimens no measurement of ¢ was needed since a 
direct comparison of specimen reluctances with and without pressure showed 
that the susceptibilities were unaltered, within the limits of the experiment 
(Jones and Stacey 1953). 

If the micrometer value, ¢’, for the thickness of a magnetic specimen was 
related to the accurate value, t, by: 


t= ?t/(1+e), 
then all of the values of permeability calculated from this were also in error: 


t_ (Ite) 1 _ Ite 


7 = => ’ 7 


a” mn MB fm 
In high fields the specimen was saturated, so that from equations 1 and 4: 


et. 1 
ee, “ae 
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whereas: 


(5) remit) ele 


so that if, instead of saturating in high fields, the specimen appeared to give a 
linear increase in J with (4M), according to equation 5, the correcting factor 
(1+) for t’ was immediately obtained. 

Frequent rupture of the inner secondary coil caused a slight revision of the 
method, so that once calibrations had been completed measurements with this 
coil were no longer necessary. A further consideration of the reluctances of the 
magnetic circuit yielded: 

(6) M, = CM,—D 





where C and D are constants for which calibrations were obtained in the same 
way as for A and a in equation 2. Tests with permalloy, nickel, and copper 
specimens showed that equation 6 was adequate for the determination of M, 
from M2. 

Calibration of the apparatus was found to be constant after the first few 
applications of pressure. As a check on the method, measurements of specimen 
flux and magnetomotive force were made with a fluxmeter and a magnetic 
potentiometer, which gave the same results within the limits of their accuracy. 


4. NICKEL AND NICKEL-COPPER ALLOYS 


The magnetizations of nickel—copper alloys with 0 to 25% copper were 
increased by pressure at all pressures and in all fields. Fig. 2, which shows the 
effect of pressure on the maximum permeability of nickel, gives an indication 
of the nature of the compression imposed upon the specimens. The graphs 
show the effects of three types of compression, linear (Ewing 1900), hydro- 
static (Yeh 1925), and the clearly intermediate type of the present work. 

Further evidence about the nature of the compression was obtained from 
the non-magnetic specimens, which were shown to remain non-magnetic up 
to 10,000 atm. (Jones and Stacey 1953). The values of ¢/u obtained with 
specimens of different hardnesses reflected the elastic compression of the 
specimens, which appeared to be governed by an elastic modulus intermediate 
between Young’s modulus and bulk modulus. Correction for this effect has 
been made to the results obtained with all the magnetic specimens. 

The saturation magnetizations of nickel-copper alloys have been found to 
increase markedly with compression, the increase depending upon the copper 
content, as indicated in Fig. 3, where the proportions of copper are given in 
atomic percentages. 

Measurements were also made on a series of six nickel-copper alloys with 
29.1 to 33.5 atomic per cent copper. These alloys showed very little ferro- 
magnetism at room temperature, which is comparable with their Curie 
temperatures. A small increase with pressure was observed in the saturation 
magnetization of the 29.1% alloy, but in the others the effect was small 
enough to be masked by experimental error. This means that within the limits 
of the experiment, +5°C., pressure was not observed to affect the Curie 
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temperatures of copper-nickel alloys with about 30% copper. Curie tempera- 
tures of nickel and Monel (basically 70Ni-30Cu) under hydrostatic pressure 
have been observed directly by Patrick (1954), who found a very small effect 
with Monel, in agreement with Michels et a/. (1937, 1941), but a much larger 
effect with nickel. 

From measurements of volume magnetostriction, w, at saturation, the 
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variation of saturated intensity of magnetization, 7,, with pressure, p, can 
be calculated, using the thermodynamic relation: 


(dw/dH)» = —(0I/dp)z- 


Calculations from values of the volume magnetostriction of nickel reported 
by Kornetski (1935) and Azumi and Goldman (1954) give —1.6 10-4 e.m.u./ 
atm. and +0.55X10-‘* e.m.u./atm. respectively. The present work gave 
+180X10-* e.m.u./atm. at 1000 atm. and +25X10-‘ e.m.u./atm. at 
10,000 atm. This must have been due to a substantial increase in the occu- 
pation of the 3d band, although the shape of the curves of Fig. 3 may repre- 
sent a gradual transition from purely longitudinal to largely hydrostatic 
compression. The results obtained with the Nilo alloys eliminated the possi- 
bility that this behavior was a characteristic of the apparatus. 


5. IRON AND COBALT 


These metals revealed one of the limitations of the present method. The 
use of magnetic intensities greater than about 18,000 gauss was prevented by 
saturation of the pole-piece tips, so that adequate saturation of iron and 
cobalt specimens was not possible. 

The effect of compression on iron in low fields is strongly dependent upon 
carbon content and thermal history. A specimen of Swedish iron subjected to 
a prolonged hydrogen anneal at 900° C. gave a positive pressure coefficient 
of magnetization in all fields and at all pressures, whereas the same specimen 
only partly annealed had given a negative coefficient in low fields at low 
pressures. Mild steel gave a negative coefficient except for the initial com- 
pression at moderate fields, but a Swedish iron specimen which had been case- 
hardened to the maximum carbon content gave a positive coefficient under all 
conditions. 

It appears that pressure increases the magnetizations of pure iron and high 
carbon steels but that at intermediate compositions, corresponding roughly 
to mild steel, the opposite effect is observed. Work-hardening has the same 
effect as increasing carbon content. This probably explains the disagreements 
between the results of early workers (Nagaoka and Honda 1898; Frisbie 1904). 

Cobalt is very little affected by pressure and measurements up to 5000 
oersteds indicate that its saturation magnetization is independent of pressure 
within the limits of the present method (+2%). In moderate fields magnet- 
ization is increased slightly by compression and in low fields it is decreased, 
the Villari reversal of behavior being a function of pressure and occurring at 130 
oersteds at low pressures. The agreement of this result with that of Chree (1889) 
may be fortuitous; the disagreement with Yeh (1925) is due to the hardened 
state of the specimen, and possibly also to the method of compression. 


6. SOME COMMERCIAL ALLOYS 


Monel and Brightray C (80Ni-20Co) (Mond Nickel Co.) are negligibly 
affected by pressure, but Nilo 36 (36Ni-64Fe) gives a strong decrease in 
magnetization in all fields, the saturation magnetization being given by: 


I, = (1060+20) —(1.14+0.05) X10» e.m.u., 
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p being in atmospheres. J, increases slightly with pressure for Nilo 50 (50Ni- 
50Fe). 

Three austenitic stainless steels based on the composition 18Cr-6Ni—76Fe 
have been tested. They are only slightly magnetic, with a tendency for magnet- 
ization to decrease with pressure. 


7. CONCLUSION 


The method has only been used for measurements up to 10,000 atm., but 
it has been shown to be capable of giving a considerable extension of the range. 
Tests on copper specimens with the pole-piece tips replaced by cobalt steel 
ones have shown that pressures up to 80,000 atm. can be reached without 
rupturing the specimens, but the magnetic hardness of these tips greatly 
reduces the accuracy of measurements. 

Correlation of the results of this experiment with measurements on magneto- 
striction presents some difficulty since the compression which has been applied 
is neither hydrostatic nor linear. Agreement with volume magnetostriction is 
poor (Lee 1955). For comparison with longitudinal magnetostriction we may 
consider the results of Kirchner (1936), who measured the magnetostriction 
of materials under pressure. He found a strong decrease in the longitudinal 
magnetostriction of nickel with increasing pressure, which is in general agree- 
ment with the decrease in d],/0p observed in the present experiment. How- 
ever, direct correlation cannot be attempted, firstly because Kirchner measured 
the changes in length caused by saturating his specimens and not the effect 
of increasing the field on already saturated specimens, and secondly because 
the effects of longitudinal compression on both longitudinal and volume 
magnetostriction were not included. The possibility of repeating measure- 
ments up to 10,000 atm. in a hydrostatic pressure chamber, using a similar 
permeameter method of attaining high fields, is therefore of considerable 
interest. 

Apart from the thermodynamic connection of magnetostriction and the 
effects of stress, one must explain the increase in saturation magnetization of 
nickel—-copper alloys in terms of decreased occupation of the 3d electron 
bands. It is possible to invent band shapes to give this effect when broadened 
by compression, but the nature of the broadening is sufficiently unknown to 
prevent a useful conclusion being drawn. Some light would be thrown on this 
by repeating the nickel-copper experiments at liquid air temperatures, at 
which the 30% copper alloys are strongly magnetic. 
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SELECTION RULES AND TEMPERATURE DEPENDENCE OF THE 
FIRST-ORDER RAMAN EFFECT IN CRYSTALS! 


By O. THEIMER 


ABSTRACT 


Starting from the most general scattering formulae, the current theory of the 
Raman effect in crystals is modified in such a way as to remove the well-known 
discrepancies between theory and experiment concerning the selection rules for 
calcite and similar crystals. A distinction is made between electrons in delocalized 
crystal orbitals and electrons in localized atomic or molecular orbitals and it is 
shown that only the latter produce a Raman scattering in agreement with 
the unmodified theory. The general formula for the scattering by delocalized 
electrons is analyzed and it is found that the magnitude of the components 


[ (Berpe / 395) 094) 544-1 of the first-order polarizability (g; normal coordinate of the 


scattering lattice vibration) depends on the wave vectors Q’ and Q” of incident 
and scattered light. The wave vector dependence of the first-order polarizability 
requires selection rules for the first-order Raman effect which do not correspond 
to the full symmetry of the scattering crystal but only to the symmetry opera- 
tions of the group of Q = Q’—Q” which leave Q unchanged. These modified 
selection rules are shown to be compatible with experiment. The influence of me- 
chanical anharmonicity and of polarizability derivatives of higher order on the 
first-order Raman effect is also discussed. It is found that these non-linear effects 
offer no satisfactory explanation for the great intensity of forbidden lines in the 
Raman spectrum of some crystals. Concerning temperature effects the non-linear 
terms in the intensity formulae are found to be of greater importance and 
are tentatively suggested as being responsible for the anomalous temperature 
dependence of low frequency external lattice vibrations. 


1. INTRODUCTION 


It has been known for several years that there are certain discrepancies 
between the theory of the Raman effect in crystals and some experimental 
results (see p. 314) which, so far, have not been explained.? It is the purpose 
of this paper to modify the current theory such as to remove this deficiency. 
We first try to find a starting point for possible modifications by a critical 
analysis of the current theory and the relevant experimental results. 

In its most elaborate form, the theory of the Raman effect in crystals has 
been given by Born and his collaborators (Born and Bradburn 1947; Smith 1948; 
Theimer 1951, 1952). The theory is a straightforward extension of Placzek’s 
theory for atoms and molecules (1934) and it leads in most essential points to 
the same results as the semiclassical theory used by J. Cabannes and other 
authors who are quoted extensively in a paper by Couture (1947, 1948). 
Born’s theory, referred to as the linear, atomic polarizability theory in this 
paper, is based upon the following assumptions: 


(a) The scattered intensity corresponding to the transition between two 
vibrational states v’, v’’ is proportional to the square of the matrix element 


1Manuscript received July 15, 1955. 
Contribution from the Physics Department of the University of British Columbia, Van- 
a British Columbia, and Summer Research Institute, Queen’s University, Kingston, 
ntario. 
*Some explanations have been tentatively suggested by several authors (Barriol and Fruhling 
1948; Matossi and Mayer 1948; Kastler 1948; Couture-Mathieu 1950a, b, 1952; Hornig 1948) 
which are, however, not quite satisfactory. 
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[M]?.., where M = aK is the electric moment induced by the electric field E 
of the incident light wave and a is the polarizability tensor of the scattering 
crystal. 

(b) M and eare considered to be the sum of the contributions of the individ- 


N 
ual atoms: a = } a(n); N is the number of atoms in the crystal. 
n=l 

(c) The frequencies »’ and v” of the incident and scattered light are far 
removed from all the resonance eigenfrequencies v> of the scattering system, 
such that the condition \pr—v’| > v%. is fulfilled. v%, is the vibrational fre- 
quency corresponding to the transition v’ — v”’. 

(a), (b), and (c) form a set of independent assumptions; they impose some 
further conditions upon the scattering theory which are listed as items (d) 
to (g): 

(d) a isa symmetric tensor. 

(e) @ is a characteristic property of the electronic ground state of the 
scattering crystal and depends on the incident light only via a resonance 
denominator. In particular, it is assumed that a becomes asymptotically 
independent of the wave vectors Q’ and Q” of the incident and scattered light, 
if |Q’| tends towards zero. 

(f) The wave vectors Q’ and Q” enter the scattering formulae only in the 
form of phase factors which take care of the retardation of the induced mo- 
ments at a distant point of observation: 


M = ED a(n) exp[27iQ . r(n)], Q = Q’—Q”, r(n) = position vector of the 


nth atom. 
(g) Only those lattice vibrations are active in the first-order Raman effect 
for which q,; = —Q. q, represents the wave vectors of the normal vibrations 


of a crystal lattice which may be represented as running or stationary plane 
waves. It is generally assumed that, for visible light, the scattering by waves 
with q, = —Q is practically identical to the scattering by the fundamental 
lattice vibrations with q, = 0. 

There are still two additional, independent assumptions contained in Born’s 
theory which are responsible for the label /inear polarizability theory: 

(h) If a is expanded in terms of the normal coordinates* q,, only the first- 
order polarizability (da/4q,)og; contributes substantially to the first-order 
Raman effect (v,’ = v,/ + 1). 

(1) The lattice vibrations are harmonic. 


The selection rules are derived on the basis of the preceding assumptions. 
It follows from (e) and (g) that they correspond to the full symmetry of the 
scattering crystal and that they do not depend on the direction of the incident 
and scattered light; according to the assumptions (/) and (i), they are exclu- 
sively determined by the transformation properties of (da/0q )ogj- 


*The notation g; and q; for normal coordinates and the associated wave vectors has been 
chosen as most closely resembling the different notations introduced by Placzek and Born. 
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The formula for the temperature dependence of the intensity of first-order 
Raman lines is based upon the assumptions (c), (4), and (z). One finds 


1 
~ 1—exp(—hv,/kT) 
jut Jao ees 1 
” exp (hv;/kT)—1 
v, = vets is the vibrational frequency of the scattering lattice vibration. 

The selection rules and the temperature dependence of Raman lines can be 
observed experimentally. Hence, the validity of the assumptions (a) to (2) 
may be checked, at least indirectly, by comparing the experimental results 
with the predictions of the linear, atomic polarizability theory. It is found in 
this way that at least some of the premisses of this theory must be abandoned. 

Perhaps the most important and most interesting argument leading to 
this conclusion is the breakdown of selection rules for calcite and similar ionic 
crystals and for organic crystals like benzene and diphenyl (Michalke 1938; 
Couture 1947, 1948; Aynard 1952; Fruhling 1951). It appears to be well 
established by the work of different authors that some of the tensor compo- 
nents [a@.]>+1 Which should vanish according to the predictions of the polariza- 
bility theory have finite values such that the magnitude of the forbidden 
components is of the same order as the magnitude of the allowed ones. Further- 
more, it is found that the magnitude of the forbidden components depends 
strongly on the direction of propagation Q’ and Q” of the incident and scat- 
tered radiation. Kastler (1948) has suggested that both anomalies have the 
same physical origin. 

Another discrepancy between the theory and experiment is observed for 
the temperature dependence of the scattered intensity. According to the 
linear polarizability theory the intensity of Stokes lines should increase with 
increasing temperature or, for high vibrational frequencies, remain practically 
constant; actually it decreases strongly in many observed cases (Landsberg 
and Mandelstam 1930; Nedungadi 1940, 1941; Ornstein and Went 1935); 
Venkatesvarlu 1941, 1942). This will be referred to as the inverse temperature 
effect in this paper. 

It will be shown in the next sections that these experimental results can be 
understood in terms of the general scattering theory for arbitrary systems of 
nuclei and electrons in which the simplifying assumptions (a) to (7) are 
abandoned. In particular, it will be found that the premisses (b) and (e) are 
mainly responsible for the apparent breakdown of selection rules and the 
assumptions (#) and (2) for the inverse temperature effect. The arguments 
leading to these conclusions are developed in the following sequence: 

Section 2.— Assuming the resonance condition (c) to be valid, the general 
scattering tensor (2.2) is transformed into the polarizability, i.e. a tensor, 
which is a characteristic property of the electronic states of the crystal and 
which depends on the positions only but not on the vibrational states (energies) 
of the nuclei. It is shown that the polarizability may be represented by the 
conventional formula (2.20) only if the electrons of the crystal move in 


Ties (Stokes lines), 


(anti-Stokes lines). 
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localized atomic or molecular orbitals while equation (2.17) is required if 
they move along delocalized crystal orbitals. In this latter case the wave vec- 
tors Q’ and Q” appear in integrals over delocalized electronic wave functions 
and, hence, determine the magnitude of the polarizability components. 

Section 3.—The difference in the behavior of localized and delocalized 
electrons shows up most clearly in the selection rules which are extensively 
discussed in Section 5. The basis for this discussion is the expansion of the 
polarizability in terms of normal coordinates given by the equations (3.3) 
to (3.8) of Section 3. The influence of the second-order and third-order terms 
in this expansion on the first-order Raman effect (v,/’ = vj/ + 1) is analyzed 
in some detail and an upper limit for the magnitude of the non-linear effects is 
estimated (subsections (a) to (d)); the upper limit is found to be big enough 
to offer a possible explanation for the inverse temperature effect which is 
discussed in Section 4. 

Section 5 deals with the selection rules. The general principles of the method 
used for their derivation are outlined in subsections (a) and (b), and applied 
to electrons in localized orbitals in subsection (c). It is shown that the selection 
rules for the fundamental (q; = 0) and Raman active (q,; = —Q) lattice 
vibrations are practically identical in the case of localized electrons and that 
the higher order terms discussed in Section 3 are very unlikely to cause 
a pronounced breakdown of selection rules. This leaves the electrons in de- 
localized orbitals as the most probable origin of forbidden Raman lines. 
The latter hypothesis is discussed in subsection (d) where the selection rules 
for the Raman scattering by delocalized electrons are derived (equations 
(5.8), (5.9), (5.11)). 

Section 6.—The Raman scattering of calcite is discussed in great detail and 
is shown to confirm the hypothesis that electrons in delocalized crystal 
orbitals are responsible for the observed anomalies. 


2. GENERAL THEORY OF THE RAMAN EFFECT IN CRYSTALS 


According to Dirac’s theory of radiation the scattered intensity [4(Q’, Q”’) 
corresponding to a transition between two states k, m of any arbitrary system 
of nuclei and electrons can be represented in the form (Dirac 1927; Fermi 
1932; Kramers 1938; ee 1954) 


2.1) 1(Q',0"de = 855 
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where 

e’,e” = polarization vectors, 

\’, \”” = wave-lengths, —- nih 

of, ww Guat, (° the incident and scattered radiation, 

Q’, Q” = wave vectors, 
r = distance between the scattering system and a distant observer, 

dw = small solid angle in direction Q”, 
I, = intensity of the incident light. 
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As indicated by the notation I%(Q’, Q’’)dw, equation (2.1) refers to the case 
that incident plane waves with wave vector Q’ are scattered and that radiation 
with wave vector Q” is observed within a small solid angle dw. 

The component [C,«]t of the scattering tensor has the form 


Bo coal opin cr ory ENE at | 
(2.2) less = ee Ope — X jie pee de eee ’ 
where 


(2.3) hvt = E,—E, is the energy difference between state r and k, 


(2.4) (P' jt = J exp(—2mi0” .r)(ov)¥ do, dv = dxdydz, 
2.5) (P) = fexp(2niQ’ .r)(ov)s a, 


(2.6) Rit = h fexpenig .¥)(p’/m)s dv, Q = Q’--Q”, 


Zi=Z. Vimy, 252 


(2.7) (ov)§ = DO] | ave eve d +, evi VT] dx, dy, dz, 
a 


2 2 
(2.8) (o°/m)s = l JiOt wv & vDIT dey ay, az, | ; 
i mM; mM; j Liz, yimy, 2i—2 
jAt 
€;, m;, V; = charge, mass, and velocity of the ith particle of the scattering 
system; only the electrons need to be considered in the case of 
visible light. 

According to the equations (2.2) to (2.8), the scattering tensor is, in general, 
complex and non-Hermitian (C,, # C,,) and a complicated function of Q’ 
and Q”. Furthermore it is a property of the whole crystal and can not imme- 
diately be resolved into a sum of atomic contributions. This indicates that the 
assumptions (a) to (e) implied in the polarizability theory can be valid only 
under certain conditions which will be analyzed in the following discussion. 

We try first to replace the general scattering tensor C,, by a tensor which is a 
property of the electronic ground state and depends on the positions of the 
nuclei only but neither on their velocities nor on the vibrational states of the 
crystal. We shall call this tensor the polarizability a of the crystal although it 
will be a function of the wave vectors Q’ and Q”. In order to introduce aye. 
into the scattering formulae we write for the total eigenfunction y, of the 
system of nuclei and electrons the product 


e = electronic 
(2.9) ve = &,.(x, g)Ue..(Q), sapiens \ quantum number, 


where ©,(x,q) represents the electronic and U,,,(q) the vibrational eigen- 
function of the system. (We need not consider rotations in most cases of 
crystals.) x and g stand for all the electronic and nuclear coordinates and the 
notation indicates that ®,(x, g) is assumed to depend on the nuclear normal 
coordinates g which act as parameters in the Hamiltonian for the electrons. 
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With the product notation (2.9) the matrix elements of equation (2.2) can 
be written 


(2.10) (23°) = EB @MSdbe = f Ube @No-T ee (ada, 


where {J>’’(q)}&.. is obtained from the formulae (2.4) and (2.7) after replacing 
y, and ve by ®,-(x,g) and &,..,(x, q). Similar expressions result for (?’); and 


for PM 'n- 
We now consider the denominators of equation (2.2) and write with an 


obvious notation 


1 1 1 wees 
(2.11) of te = yore ev ere 4. yp’ = a =. yp’ TT Gis he ypt ce 


This allows one to split the scattering tensor in the polarizability* 
’ 1 ° etal 
(2.12)  [ape(g)]o" = se Rte 


oe ee (zn ( ) eed pel E ‘( yye | 


” ver mye 


42 eel @lelenl Ele ‘Ie “M) 


roby’ 


, 
e 


and in the tensor 
2 1 ALE iby e’ ro voll ot eo’? 
, ’ i -_ Ne GN ee 
(2 13) [Bre (g v)] Sav v 1 ZA (v' ee ) 


ait ree on FenIp eet ah 
v’) 


which depends on the vibrational frequencies v}.., as ell This dependence 
leads to rather complicated formulae for the temperature dependence of 
Raman intensities if the condition 1(c), i.e. the condition a(q) > B(q, v), is 
not fulfilled. 

The summation over v’” in equation (2. 12) can be removed by means of the 
rules of matrix multiplication ([P(q)]?) = X1-[A(@]?--[B(@]%, P = AB) 
and one obtains 


nr »’ 
(2.14)  [arpe(q)]? =a sal Hi -¥ {2 (g)]e (Je (qe 


Vern —y 


[eg Nel are} | 
vero ty ” 


a(qg) is a function of the nuclear coordinates g and depends, even in this 
simplified form, on the wave vectors Q’ and Q” in a way which is more com- 
plicated than postulated by the assumption 1(f) of the polarizability theory. 
This assumption is based upon the possibility of representing a(g) as a sum of 
*In equation (2.12) and in all the following formulae, it is assumed that the electronic 


eigenfunctions ®,’ and ®,’’ of the initial and final state of the scattering system are identical: 
e’ = e’’, At reasonable temperatures they will refer to the electronic ground state. 
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atomic contributions. Hence, we shall now analyze the conditions under which 
such a representation is possible. 

We shall show that these conditions are related to the question whether the 
electrons of crystals are restricted to move in localized atomic or molecular 
orbitals, or whether they may travel, more or less freely, in crystal orbitals 
spread over the whole crystal. Since the idea of localized atomic and delocalized 
crystal orbitals refers to two extreme limiting cases which are hardly found in 
nature, one can not generally assume that some orbitals are completely 
localized and others are not. We shall therefore represent all orbitals as linear 
superpositions of localized, non-overlapping, atomic orbitals and plane waves, 
the latter having coefficients which measure the degree of delocalization.* 
Not much is known of the magnitude of these coefficients for ionic crystals 
of the type CaCO; but the widely accepted assumption that the ionic bond is 
never perfectly realized, not even in NaCl, indicates that they can not, a 
priori, be neglected in a general theory of the Raman effect in crystals. 

In order to introduce the concept of localized and delocalized orbitals into 
the scattering theory we represent the electronic eigenfunction ®,(x, g) by an 
antisymmetric sum of products of one-electron eigenfunctions: 


(2.15) (x, q) = ~H¥ (—1)"Pwl(x, qwi(e, g)... wile, g)... wilt(x, 9), 


VN! 
with 
N, = total number of electrons in the scattering system, 
in = number of inversions in the permutation P, 
wi(x, g) = orbital / occupied by the electron i; the w;(x, g) are orthogonal 


and normalized and depend on q. 


Since the number of. available orbitals is much larger than the number N, 
of electrons, different functions ®,(x, g) are characterized by a different set of 
occupied orbitals w;. 

We now introduce, in place of y and y,, electronic eigenfunctions ®,-(x, q) 
and #,,,,.(x,q) of the form (2.15) into the matrix element (2.7). This matrix 
element is a sum of terms, each of which depends on the velocity of one 
single electron only. Hence, this matrix element vanishes unless ®,--. (x, g) and 
®,,(x, g) differ by not more than one of the occupied orthogonal orbitals w;. 
More explicitly, one finds that 


(2.16)  (pv)ée = (pv)ier = 4[wih- (x, gevwy (x, g) wr (x, ghev* we (x, g)], 


where wy (x, g) and wy--(x, g) are those occupied orbitals which are different 
in the initial (e’) and intermediate (e’’’) state. 
Equations (2.16) and (2.14) combined give an expression for the polariza- 


*It is always possible to represent the electronic eigenfunctions of a periodic tal exclu- 
sively in the form of plane waves. In this case, the mobility or delocalization of the electrons 
may be measured by the group velocity dE/d(k/d). It is, however, more convenient for our 
purpose to reserve the plane wave representation for delocalized orbitals with group velocity 
different from zero. 

tThe absence of electronic conductivity does not contradict the assumption of delocalized 
orbitals which is a necessary, but not a sufficient, condition for metallic conductivity. 
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bility in which the summation over electrons is replaced by a summation over 
the initially occupied orbitals wy (x, q): 


(2.17) [erye(q)]¢" = wil st J exp2nid.n)(#)’ af 
<< jfexo(= 2miQ” . £) (o0,) dv fexp(2xiQ . r) (ove) i” 


AO 


4 fexp(2niQ' . £) (pve) ido fexp(—2miQ” . r) (p0,)'r'dv x Nivaeh |” 
vor py! ae 

Equations (2.16) and (2.17) have a form which allows one immediately to 
visualize the different behavior of localized orbitals (w,') and delocalized orbit- 
als (w,*). According to the discussion on page 318 we write for each orbital 


(2.18) w(x, g) = w,'(x, g) +cw,*(x, g), 


and substitute (2.18) in (2.16); this gives a sum of terms containing three 
different types of products of eigenfunctions, namely: w'w', w'w, and ww’. 
The last product is a space function spread over the whole crystal while the 
first two products vanish everywhere outside the charge cloud of the localized 
orbital. It is essential to note that the mixed products represent localized 
space functions; hence, (pv); will contain localized terms due to the initially 
occupied orbitals w}, even if all the exited orbitals w,», should be essentially 
delocalized. 

The contributions of the localized terms to the polarizability can be for- 
mally simplified by grouping together all the localized orbitals w,,) belonging 
to the same atom () and by introducing the position vectors r() of the atomic 
nuclei with the help of the relation 


(2.19) r =r(7)+r(n), 

where r(7) is the vector from a nucleus to an electron in the localized orbital 
W y(n) (x, q). Since the linear dimensions of atomic chargecloudsare negligibly small 
compared to the wave-length of visible light the phase factor exp[27iQ’ . r(’)] 
in (2.17) is practically equal to one, over the whole range of integration; while 
the factor exp[27iQ’ .r(m)] may be taken before the integral, r(m) being 
constant by definition. We obtain in this way the following expression for the 
polarizability a'(q) of electrons in localized orbitals: 


(2.20) [ap.(q)]% = Do exp[27iQ . (2) }ape(n)]¥: 


[ape(n) i» = rea =e 4 >» i a /m) > (nd ot, 


-2)> Soe aes saa J (oo. dv | (o0,)' iodo} |" 


Gay «br? 


with 


This may be finally transformed oe 1938; Placzek 1934) into an expres- 
sion resembling | Placzek’s formula for free molecules: 


*The imaginary part of ag¢(n) which is usually negligible has been omitted. 
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(2.22) [erpe(m) Jo 


-? i : TW) ia aa (M, yisPdy | (Me) irqndv . (M=ae) 
hLitm v7 OY ee 


The equations (2.20), fe - (2.22) correspond exactly to aR assumptions 
1(a) to 1(g) of the polarizability theory. However, they are valid only as far 
as the electrons in localized orbitals and their contributions to the polarizability 
are concerned. 

The situation is quite different for the contributions a‘(q) of electrons in 
delocalized crystal orbitals to the total polarizability a. For these contributions 
the general formula (2.17) can not be further simplified. In particular, the 
phase factors can not be written outside the integrals over electronic orbitals 
and the wave vectors Q’ and Q” influence the scattering intensity even if 
|Q’| tends towards zero. This means that the tensor components [a‘, (q)]*- 
depend on the direction of propagation of the incident and scattered light even 
if the wave-length is long compared to the linear dimensions of the elementary 
lattice cell. 

From the preceding discussion of the general scattering formula (2.1) it 
appears that it is the polarizability a‘(g) of electrons in delocalized crystal 
orbitals which is most probably responsible for the breakdown of selection 
rules mentioned in the introduction. We shall elaborate this hypothesis in 
some more detail in Section 5. 

3. EXPANSION OF THE SCATTERING TENSOR IN TERMS OF NORMAL 
COORDINATES 

The problems of light scattering are most readily understood by expanding 
the scattering tensor a,,(q) in terms of the normal coordinates qg, of the 
scattering system. 

We write 


(3.1) Uy = te, (qr) oq (G2) . . » Uo;(Qs) - - - Uogy (Qa) (N = number of nuclei) 


for the vibrational eigenfunction, and 


(3.2) E,= x (v;+4)hvy 


for the vibrational energy of a crystal in its electronic ground state, and we 
expand ay¢(g) in terms of the q;: 


83) ae(0) = (not D (#e). at hd (Pee 


pee of q3' 


wer ( Site a qy 
62s " ak. oaent 


(The subscript zero refers to the equilibrium position of the nuclei.) 
In the case of crystals with localized electrons the terms in the expansion 
(3.3) have the following form (Born and Bradburn 1947; Theimer 1951): 


(3-4) ye(j) = (=) qs = (x DX aren(beny(b)expl2#i(G,+Q) . ran 
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(3.5) Ge(i.i") = (ees a = des ar - ( LDL coeim (ta Yeas lh deny (R') 


i—i%) Ree 
Xexp{2zi[(q,+qy+Q) . r()—(qy—Q). rl’) )a, qs’ 


The index / enumerates the different lattice cells and the index & the different 
nuclei in one cell. e,,(k) exp[2miq, . r(Z)] is the normalized u-component of the 
displacement of the nucleus (&) in the cell (/) with position vector r(/), corre- 
sponding, to the normal vibration g,; with wave vector q,;. The ap.,4(k) and 
Ope uy’ 5’ ') are the first and second tensor derivatives with respect to these 
nuclear displacements. The a,,,,() do not depend on the cell index / and the 


Ope py’ Ge") depend only on the difference (/—/’). Hence the 2 -(7) and 
‘Suli j’) vanish unless q,+Q = 0 or q;+q,+Q = 0. This is a consequence 
of the summation over /. (For more details see Born and Bradburn 1947; Smith 


1948; Theimer 1951, 1952.) 
The scattering intensity is determined by the matrix elements of the 
scattering tensor. For the case of harmonic vibrations one finds: 


(3.6) [ape]2; = att (vy +1)( 224 7) by + Streets b,= V sts ’ 
(3.7) laplas = Veert{ (2) bt Sey (ars Bae Nit...) 


(3.8) [expe ) 3. o/+ = +7’ a I (+4) TT ( (opti? (2 en )i'by Fsace 


We shall now shortly discuss some problems related to the formulae (3.6), 
(3.7), and (3.8) which are of interest for the understanding of selection rules 
and of the temperature dependence of Raman intensities. 

(a) The influence of the third-order terms in the first-order matrix element 
(3.7). 

One can try to find an upper limit for the contributions of the third-order 
terms to the intensity of first-order Raman lines from the fact that the same 
terms determine the intensity of the third-order Raman effect (v,’ = v,/+3 
andy,’ =9,/+1,vy” = vy’+2). 

We write for the intensity of one first-order Raman line 


« da 
(3.9) I(j) = Ivar ~ flepela}? ~a(lty); a= (2s, 
a> ey t1)(5 aon) 4,64 


(y measures the relative contribution of third-order terms to 7(j)) 
and for the intensity of the whole third-order Raman spectrum 

(3) 
(3.10) I~a’y’, 


and obtain 
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(3) () 2 


(3.11) I/I(j) ~ G45) 


One should now distinguish between the case of molecules where the third- 
order spectrum is a line spectrum and the case of crystals where it is a con- 
tinuum. The absence of third-order Raman lines in molecular spectra indicates 
that y < 1 and justifies the assumption 1(h) of the linear polarizability theory. 
In crystals, however, one can not definitely exclude the possibility that the 


(3) 
total intensity J of the continuous third-order spectrum is of the order of 
some per cent of the intensity of one first-order line. Such a small effect, 


spread over a wide frequency region, could hardly be separated from a weak 
3) (1 
continuous background due to other causes. Thus, assuming 1/T(j) < 0.03, we 
find: y?/(1+-)? < 1/30 or y < 0.2. This high value is an upper limit and, 
hence, offers no decisive argument against the assumption 1(h). But it is of 
interest as a possible explanation for the inverse temperature effect, as it allows 
for strong higher order effects with complicated temperature dependence. 
It should, however, be stressed that the third-order terms in equation (3.7) 
can be hardly responsible for a great intensity of forbidden Raman lines, 
even if y has as high a value as 0.2. The reason is that the first-order term a 


1 
of equation (3.9) is zero in this case and the intensity 7 ( j) of the forbidden 
Raman line is of the order y?, i.e. at the most 4% of the allowed first-order 
lines. Only if a ¥ 0 are the third-order contributions enhanced to the value 
27 by forming the square of (1+7). 

(6) For the understanding of temperature effects it is also useful to note that 
the formulae (3.6) to (3.8) are not valid for the matrix elements [8,¢(g, v)]? 
(see equation (2.13)) of this part of the scattering tensor which is negligible 
if the resonance condition 1(c) is fulfilled. Those matrix elements depend on 
the vibrational quantum numbers v, in a rather complicated way which is 
different in each individual case and which strongly modifies the formula for 
the temperature dependence of Raman lines. Hence, breakdown of the con- 
dition 1(c) may produce an inverse temperature effect. 

(c) The equations (3.6) to (3.8) refer only to the case of harmonic vibrations. 
The mechanical anharmonicity of the lattice vibrations introduces the follow- 
ing modifications: 

(i) The formulae for the matrix elements of anharmonic oscillators are 
different from those for the harmonic oscillator. As far as the matrix element 
[aye]e441 is concerned one finds quite generally that the odd derivatives which 
are contained in (3.7) decrease and that the even derivatives, for example 
[(0?ap6/9q 599 y)0979y}e44+1, Which are zero in (3.7) appear and increase with 
increasing anharmonicity. They have probably the same magnitude as the 
odd terms of next higher order. 

(ii) [otpe] 3441 is a function of the vibrational quantum number v,. As far as 
only the term [(Aay6/9q;)09s}s441 is concerned, this function is monotonously 
increasing with v, and tends asymptotically towards a limiting value (Buck- 
master 1952). One has here an analogy to the classical behavior of the vibra- 
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tional amplitude which increases with increasing vibrational energy until 
dissociation occurs. This special property of the function [(dap¢/q,)o9s]1}41 = 
f(v;) is important for a discussion of temperature effects; it warrants that 
only the higher terms in the expansion of a(g) can be responsible for a decrease 
of [ope]2441 with growing v;. 

(iii) As a result of the mechanical anharmonicity the transition frequencies 
vej+1 Of a certain lattice vibration (g,;) depend on the quantum number », 
and tend towards zero with increasing v,;. This produces a strong broadening 
of Raman lines, particularly of those corresponding to highly excited, low 
frequency vibrations of ionic lattices, where the ions are known to move in 
strongly anharmonic potential wells. The broadening, again, reduces the peak 
intensity of Raman lines and may even produce an apparent loss of the integral 
intensity if the wings of broad Raman lines escape observation or are confound- 
ed with a continuous background. 

(iv) A rather important deviation from the case of harmonic normal vibra- 
tions is the hindered rotation of certain atomic groups, e.g. the CO;~-ion in 
calcite, which is a consequence of the existence of several different potential 
minima. We discuss shortly this case as it may have an influence on the 
selection rules for the vibrational Raman effect. 

The problems of hindered rotation can be understood from a discussion of 
the two limiting cases: anharmonic torsional vibration and free rotation. 
Only the latter case introduces some new aspects into the theory of the Raman 
effect. 

We write for the eigenfunction of the rotating group the product 


(3.12) v¥(g,¢) = U(q) R(¢), (@ is angle of rotation about the z-axes) 


where the rotational eigenfunction is given by 
(3.13) R(¢) = Ce*™, (C = normalization factor) 


and we replace the Taylor series (3.3) by the Fourier series 


(3.14) Ope (G, &) = >, ayrer(q, se, (s = integer). 


Some of the a,-.(g, 5s) may be zero because of symmetry. The subscripts 
p’,o’ in (3.14) refer to a coordinate system attached to the rotating group; 
hence, equation (3.14) describes an intrinsic ¢-dependence of the polarizability 
a which is due to the influence of the neighboring atoms and ions surrounding 
the rotating group.* For the understanding of the selection rules it is important 
to note that a,-.:(g, 0) corresponds to the equilibrium position of the nuclei 
and to the full symmetry of the crystal lattice. Most of the other a,-./(g, 5) 
with s # 0 are used to represent the polarizability of the rotating group in 
non-equilibrium positions with lower symmetry. This is quite different in 
the expansion (3.3) where all the terms refer to the equilibrium position. 
We come back to the problem of rotating groups in crystals in Section 5 
dealing with the selection rules. 


*This influence, clearly, excludes a completely free rotation; but the simple assumption of a 
free rotation is sufficient for a qualitative discussion. 
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(d) It is well known that the higher polarizability derivatives and mechan- 
ical anharmonicity of lattice vibrations give rise to a continuous second- 
order Raman scattering which contains overtones and combination tones of 
all lattice vibrations irrespective of their wave-length. This opens up the 
possibility of interpreting a forbidden first-order Raman line with frequency », 
as a second-order continuum with a sharp intensity peak at v;. Excluding 
accidental coincidence the line-like second-order peak must correspond to a 
combination band involving the branch of lattice vibrations with limiting 
frequency v, and elastic vibrations with limiting frequencies », near zero. 
Several authors (Barriol and Fruhling 1948; Hornig 1948) have tentatively 
adopted this interpretation of forbidden Raman lines without closer inspec- 
tion, thereby overlooking the fact that no sharp intensity peak can be found 
with a combination band in the frequency region v,+»,(v, — 0). 

This can be shown in the following way: We assume that v, refers to a 
branch of internal lattice vibrations with frequencies practically independent 
of the wave-vectors and that the intensity of one single second-order transition 
Vj,Ve > Vj+1, ve+1 involving long elastic waves (index e) as one partner 
depends only on the direction but not on the length of the wave vectors q, 
tending towards zero; this latter assumption can be rigorously proved to be 
correct by means of the detailed formulae for the second-order Raman effect 
derived by Born and collaborators (Born and Bradburn 1947; Smith 1948). 
With these assumptions the total intensity J(v,;+»,) of the combination band 
in the frequency region round (v,;+¥»,) is proportional to the number N(»,)dv, 
of elastic vibrations with frequencies between », and »,+d»,. In general the 
frequency distribution in q-space is not known and the theoretical determina- 
tion of intensity peaks is a very difficult problern (Van Hove 1953). However, 
for the simple case of long elastic waves, N(v,) is well known to be 


2 
(3.15) N(v,)dy, = 1 dv. 


3 = mean velocity of sound. N(v,) vanishes for », = 0; hence, the intensity of 
the combination band v,+», (v,— 0) is zero in the frequency region v, No 
intensity peak of the combination band is to be expected there. 


4. THE TEMPERATURE DEPENDENCE OF THE FIRST-ORDER RAMAN EFFECT 


Owing to time dependent thermal perturbations, a crystal in thermal 
equilibrium is in a non-stationary state which may be represented as a linear 
superposition of stationary vibrational states. The squares of the expansion 
coefficients oscillate irregularly around average values which are given by 
the Boltzmann distribution. As far as light scattering is concerned and if a 
small interdependence of thermal and electromagnetic perturbations is neg- 
lected, a crystal in such a state behaves essentially like an ensemble of indi- 
vidual scattering normal vibrations g, with a Boltzmann energy distribution. 
The intensity scattered by this ensemble may be obtained from equation 
(3.7) with the following, rather obvious modifications. 


(a) We consider a Raman line with frequency », scattered by the normal 
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vibration q;. The total intensity /(v,) of this line may be represented as a sum 
I (v3) = Lo; Tei41, where J(4,, is the intensity scattered by the transitions with 
initial vibrational level v,. J;j,1 is proportional to the square of the matrix 
element (3.7)* if (vy +1) is replaced by the thermal average 

(4.1) Gyan = 2 “(vy-+1) exp(—vyhvy/kT) _ 1 

; ; exp(—vyhvy/kT) ~ 1-exp(— hvy /RT) 

(6) The intensity Z(v,;) of a Raman line with frequency v, corresponding 
to all the different transitions v;—> vj, (v; = 0,1,2,...) is obtained from 
the intensity /;4,1 by introducing the thermal average @,+1) instead of one 
special number oka This gives 


ae SRI Manat an |) 
(4.2) 10) ~ eae TB Pe 


1 
+ La 6[1—exp(— 





i i 2 
oan ae aqjaq,/7) 8 °r + -- , 


The exponential outside the curly brackets produces an increase of the 
intensity with increasing temperature which can only be inverted if the 
third-order terms in the curly brackets have negative signs and are not negli- 
gible compared to the first-order term. In order to estimate the influence of 
the third-order terms we introduce an average frequency vy = 3X10" sec.— 
corresponding to a wave number v¥ of 100 cm.—'. This choice implies the 
assumption that the highly excited, low frequency vibrations of the elastic 
branches of the crystal spectrum contribute _ essentially to the third-order 
terms. We now investigate the cases where vs = 500 cm.—! and 1000 cm. 
and y (see equation (3.9)) = 0.1, 0.05, 0.02, al 0.01 at very low temperatures 
(o» = 1). The resulting temperature dependence of J(v,) in arbitrary units is 
given in Fig. 1; it shows that under extreme conditions the third-order terms 
may be responsible for an inverse temperature effect. The effect could be 
further enhanced by reducing the value of v;. 

According to equation (4.2) and the numerical results of Fig. 1 the high 
frequency lattice vibrations are more likely to show an inverse temperature 
effect than the low frequency vibrations for which the exponential outside the 
curly brackets is decisive. This theoretical prediction is not confirmed by 
experiments (Nedungadi 1940, 1941; Ornstein and Went 1935; Venkates- 
varlu 1941, 1942). One may, therefore, hesitate to ascribe the inverse tem- 
perature effect to the contributions of third-order terms which, for unknown 
reasons, should always be opposed in sign to the derivative (da,,/dqo)0; 
but one can not completely exclude this possibility. The situation remains 
essentially unchanged if the mechanical anharmonicity of the lattice vibra- 
tions is taken into account; this follows immediately from the discussion in 


subsection 3(c). 
In addition to the contributions of the higher-order terms, line broadening 


*[etpe 33 +1 determines the intensity of polarized Raman lines produced by polarized incident 
light. In ‘the case of natural light summation over p, ¢ is required. 
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Fic. 1. The temperature dependence of the intensity of first-order Raman lines for different 
values of v? and y. ve = 500 cm.~ (full curves); »* = 1000 cm.~! (dotted curves). 


and resonance may be responsible for the inverse temperature etfect. Con- 
cerning the influence of resonance the following considerations are helpful 
for an estimate of its magnitude. 

(i) The relative magnitude of the tensor component [8,.(q, v)]? which, 
in principle, could produce an inverse temperature effect, can be estimated 
from the approximate relation (compare equations (2.11) and (2.13)) 


[Bre(q. v)}i-4? _ { vw 

~ lae(Q)Ie" J ~ Ue 

where »}, and v%. are vibrational and electronic absorption frequencies of the 
crystal. A contribution of 8(q,¥v) to the intensity of first-order Raman lines 
which is strong enough for producing an inverse temperature effect requires 
to be at least of the order 0.1. The consequences of this condition can be seen 
from the following example: We assume that v}» = 310" sec.-! and »’ = 70 
10" sec.-! (A’ = 4300 A), and obtain from the condition » < 0.1 that vf, 
is about 80X10' sec.-' (Af. = 3800 A). This estimate shows that for the 
scattering of blue light absorption in the near ultraviolet region between 3800 
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and 4000 A implies a breakdown of the resonance condition 1(c) which may 
be partly responsible for the inverse temperature effect. 

(ii) If the absorption frequencies of a crystal are unknown, the validity of 
the resonance condition 1(c) can be checked indirectly in different ways. 
Most useful for this purpose is the relation a,, = a,, expressing the symmetry 
of the scattering tensor and the formula for the intensity ratio of Stokes and 
anti-Stokes Raman lines: 





I(v,)Stokes (v’—v,)* 
__ 4) ee Ee 
£4) Tis, anti eeas” Gran, 

Both relations are a consequence of the resonance condition 1(c) and become 
invalid if the light frequency approaches an absorption frequency of the 
crystal (Placzek 1934). 


5. SELECTION RULES FOR THE FIRST-ORDER RAMAN EFFECT IN CRYSTALS 


(a) General Principles 

Selection rules are rules derived from the symmetry of the scattering 
system according to which certain components of the tensor [a,.]¢ are zero 
and, hence, do not contribute to the intensity of the Raman scattering. 
The vanishing of certain tensor components can be understood and predicted 
from the following argument: The matrix elements [a,.])- are completely 
determined by the Hamiltonian function of the scattering system and must, 
therefore, have the same numerical values in all equivalent coordinate systems; 
i.e. in all those coordinate systems which may be transformed into each other 
by the covering operations which leave the Hamiltonian of the scattering 
system unchanged. In terms of group theory this may be expressed as follows: 
We consider, for convenience, those matrix elements [a,.])- and those of 
their linear combinations which have simple transformation properties, 
€.g. (a@rrtayy), (@rr—Gyy), G22, Ary, Ayz, Wz, in the case of rotations about the 
z-axis.* They must vanish, if they do not transform according to the totally 
symmetric, irreducible representation of the space group of the scattering 
system. More generally this is true if they belong to a reducible representation 
in which the totally symmetric irreducible representation is not contained. 
The question whether or not a representation is totally symmetric in this 
general sense can be answered by means of the standard rules for reducing 
reducible representations and by applying the rule that the character of the 
representation of a direct product is equal to the product of the characters 
of the individual representations. Furthermore, one has to use the rule that the 
derivative (0"a,,/dq")» transforms like the product x,x.q" (x,, x., . . . are the 
coordinates). 


(b) The Representation of Space Groups (Bouckaert, Smoluchowsky, and 
Wigner 1936; Seitz 1936) 
We consider first the Abelian subgroup of the translations. Suitable repre- 


*In this case the linear combinations could be avoided by introducing the circular coordinates 
2, x+4y, x —ty, which will be denoted by the subscripts 0, +1, and —1. It can be shown that 
the circular tensor components a,y (u, » = 0, +1) are multiplied by the factor exp[#(u—») ¢], 
if the coordinate system is rotated through the angle ¢ about the z-axis. It can be further 
shown that (az2+ayy) = 2anand (a::—ayy) = —(ai2+a21). 
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sentations of this group are the matrices corresponding to the displacement 
by the identity period a; (4 = 1, 2, 3); they are diagonal and have the diagonal 
elements M,, = exp(2miq .a;,). Since the eigenfunctions of a mechanical 
system form bases for the representations of the group of operations which 
leave the Hamiltonian of the system unchanged, it follows, for instance, 
that the electronic eigenfunctions of a crystal with nuclei fixed in lattice points 
may be characterized by a wave vector q. The same is true for the vibrational 
eigenfunctions as long as the equilibrium positions of the nuclei form essentially 
a periodic lattice. 

We come now to the non-translatory elements of the space groups. These 
elements may be classified according to their effect on the wave vectors q 
which characterize the eigenfunctions Yq of a system. For each Yq one can 
distinguish between those operations which change q and those which trans- 
form q in itself (or change only its sign).* The latter operations form a group 
which is called the group of the wave vector q or q-group. ~q transforms 
under the operations of the group of q by an irreducible representation of 
this group which is called a small representation of the space group. 

If q = 0, the group of q is isomorphous with the crystallographic point 
group of the crystal and the full symmetry of this point group determines the 
displacement vectors of the normal vibrations and the optical selection rules. 
For general directions of q, however, the group of q contains only the identity, 
and the displacement vectors are arbitrary, i.e. only determined by force 
constants and masses. No selection rules can, in principle, be obtained from 
the non-translatory symmetry elements for these cases. Practically, of course, 
it may be that the selection rules valid for q = 0 are still a good approximation 
for cases where |g| < 1/|a,|. Whether or not this is true is a physical problem 
which will be discussed later. 

(c) Selection Rules for the Scattering by Electrons in Localized Orbitals 

(i) Selection Rules for the Wave Vectors 

Any component (da,,/0qg;)oqg; of the first-order polarizability contributes 

to the matrix element [pe] 9441 a term which transforms by a representation 
with character exp[27i(q,;+Q) .a;]. This follows quite generally from the 
transformation properties of direct products (see the rule at the end of Section 
5(a)), if one considers that a,, has the character exp(271Q . a,) (see equation 
(2.20)) and that g, has the character exp(277q, . a;). The same result may be 
obtained without reference to group theory from the equation (3.4). This 
first-order term is invariant against all possible translations only if 
(5.1) q;+Q = 0 or q; = -Q. 
This selection rule is often interpreted as representing the law of conservation 
of momentum for a collision between a photon and a free phonon. Concerning 
the contributions of the third-order terms in equation (3.7) the same selection 
rule can be obtained in the following way: 

The third-order terms transform by a representation with character 
exp[27i(Q+q,;+qy7+4q,)]. This requires the selection rule 


*The transformation q — —q does not change the geometry of the system; for scattering 
-oblems it is equivalent to a reversal of the optical path of the incident and scattered light. 
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(5.2) Q+q;+Gy+qy = 0. 


The ambiguity in the sign of the last term can be removed by the condition 
that selection rules must be invariant against the transformation Q —- —Q. 
This condition requires the negative sign in the last term and leads again to 
the selection rule (5.1).* 

The selection rule (5.1) is stringent only for harmonic lattice vibrations. 
In the case of mechanical anharmonicity one has to consider the behavior 
of polarizability derivatives of even order, e.g. (0%a,./0q,0¢)oqq7. It is 
readily found that they are subject to the selection rule 
(5.3) q; = —Q+q,, 
according to which lattice vibrations with any possible wave vector q, can 
be active in the first-order Raman effect. The contributions of the anharmonic 
terms (0%a,./0g,0q¥)09q produce a very weak Raman spectrum of the first 
order which is continuous, since the frequency of the lattice vibrations q, 
belonging to one certain branch of the lattice spectrum depends more or less 
strongly on the wave vector q,. Only if this continuous first-order spectrum is 
concentrated within a very small line-like frequency interval may it become 
intense enough for observation. This condition shows that the internal vibra- 
tions of molecules or polyatomic ions in crystals which have practically the 
same frequency for all wave vectors are the only ones which, in case of extreme 
mechanical anharmonicity, may produce a line-like first-order Raman effect 
for which deviations from the selection rule (5.1) and others become observ- 
able. This effect is equivalent to the anharmonic coupling of ‘individual 
molecules’”’ in a crystal with the surrounding crystal field which has been 
discussed by Hornig (1948). 

It should be pointed out that the selection rule (5.3) violates the law of 
conservation of momentum for collisions between photons and individual 
phonons. This difficulty can be easily overcome, if anharmonic lattice vibra- 
tions are pictured as phonons which are not free but coupled to the lattice 
by the anharmonic terms in the potential energy. In this case the law of 
conservation of momentum is only valid for the scattering system as a whole 
but not for collisions between photons and individual ‘‘bound”’ phonons. 


(ii) Selection Rules for the Matrix Elements [aye}s/1 

The selection rules for the wave vectors are valid for all the tensor compo- 
nents [a,«]s341 independent of the subscripts pe. The symmetry elements of 
the group of q, give rise to additional selection rules which are based on the 
different behavior of different components [ape};/,1 under the transformations 
of the group of q,;. The problem of physical interest is to find the group 
of q;. 

We consider first harmonic lattice vibrations for which the q-selection rule 
(5.1) is valid. In this case the question arises whether, for the purpose of 
deriving selection rules, the vibrations with q,; = —Q may be practically 

*The positive sign in the last term gives for positive Q: wt) = —(Q+2q;) and for 


negative Q: q;(—) = (Q—2q;’). The resulting inequality |g;(+) | #|¢;(—)! is forbidden by the 
above condition. 
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identified with the fundamental lattice vibrations with q; = 0. The following 
arguments seem to justify this approximation. 

It is generally assumed that all properties of crystal vibrations approach 
asymptotically the properties of the fundamental lattice vibrations if |q| 
tends towards zero. Particularly, it is assumed that different lattice vibrations 
have very nearly identical optical properties if they belong to the same branch 
of the lattice spectrum and if their wave-lengths are larger than the wave- 
length of visible light. The theory of lattice vibrations justifies this assumption 
as long as only short-range interactions between vibrating lattice points are 
taken into account. However, calculations of Blackman (1952) indicate that 
the long-range electrostatic forces in ionic crystals of lower than cubic sym- 
metry may strongly perturb the asymptotic behavior of lattice vibrations 
such that the restoring forces 0?P/dq,? (P is potential energy) of electrostatic 
origin depend on the direction of q, even if |q,| tends towards zero.* It appears 
that this effect has been discovered experimentally by Mathieu and Couture- 
Mathieu (1950) who mention cases where the frequencies of Raman lines 
depend on the direction of the incident and scattered light, and hence, accord- 
ing to equation (5.1), on q,. It should, however, be mentioned that such 
frequency variations do not occur in calcite and most of those crystals which 
have forbidden lines in their Raman spectrum. 

We come now to the asymptotic behavior of da/dq, whieh i is important for 
the derivation of selection rules. As long as only electrons in localized orbitals 
are concerned the polarizability of an atom or polyatomic molecule depends 
mainly on the position of the immediate neighbors and is insensitive against 
the electric field due to the displacement of distant ions; this follows from the 
fact that the polarizability of molecular systems is practically independent of 
weak external electric fields. We conclude therefore that da/0q, is less sensitive 
against variations of q, (|g,| +0) than 0*P/dg7 and that it should be a constant 
independent of q, if the same is observed for the frequency »(qg;). To make this 
point quite clear we consider the totally symmetric vibration of the CO;~~-ion 
in calcite which offers a good testing ground for the theory of the Raman 
effect in crystals. This vibration is a typical “internal’’ vibration such that 
the vibrational distortion of the individual CO;—~-ions and their neighborhood 
is completely independent of q, if |g,| tends towards zero. The same follows 
necessarily for the polarizability derivative (da/0g;)oq; if only electrons in 
localized orbitals are involved. In this case, the selection rules for the scattering 
by lattice vibrations with q, = —Q are practically equivalent to the selection 
rules for the fundamental lattice vibrations in the sense that Raman lines 
which are strictly forbidden by the latter are extremely weak according to the 
former. 

As far as the first-order term of equation (3.7) is concerned the selection 
rules for the fundamental vibrations have been worked out by several authors 
~ *Blackman’s results refer to wave-lengths long compared to the lattice cell but small 
compared to the macroscopic crystal. They have nothing to do with the fact, well known since 


the work of Lyddane and Herzfeld (1938), that the shape of an ionic crystal determines the 
frequencies of those lattice vibrations the wave-lengths of which are of the order of the crystal 


dimensions. 
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(Brester 1923; Cabannes 1940; Mathieu 1945; Hornig 1948) and need no 
further discussion; they are contained in the selection rule (5.11). Regarding 
the contributions of the third-order terms the selection rules are determined 
by the symmetry character of the third-order derivatives (0%a,,/0q,0q,7), 
which transform like the product xyqg,qy?. If gy is a non-degenerate normal 
coordinate, then g,? transforms always by a totally symmetric representation 
with character 1 and the product xyq,gy? transforms like the product xyq;; 
this means that the selection rules for (da¢/0q;)o and (0°a,./0g¢;0qgy")o are 
identical. However, if g, is a degenerate normal coordinate, g;? belongs toa 
reducible representation which, generally, does not reduce to a totally sym- 
metric representation alone (Rosenthal and Murphy 1936). Hence, symmetry 
character and selection rules are different for (day¢/0q;)o and (0%a,./0q;0q "?)o if 
gy is degenerate; accordingly, first-order Raman lines which are forbidden 
by the /inear polarizability theory may have a non-vanishing intensity pro- 
duced by the third-order terms in equation (3.7). 

We shall now investigate how the selection rules are modified by the 
mechanical anharmonicity of the lattice vibrations. The vibrational eigen- 
functions of a crystal form bases for the representations of the space group 
of the crystal independent of their harmonic or anharmonic nature. Hence, 
the selection rules for the fundamental lattice vibrations are valid for harmonic 
and anharmonic vibrations alike. The only new features introduced by the 
mechanical anharmonicity are the contributions of the terms (0%a,./dg¢,0q¢¥)o 
to the first-order Raman effect. Since the q-selection rules for these terms do 
not exclude any possible wave vector they belong mostly to a group of q 
corresponding to a general direction of q. No selection rules can, therefore, 
be obtained which require the vanishing of these terms for reasons of sym- 
metry. 

The situation is similar for the case of a hindered rotation of certain atomic 
groups in the crystal. For the extreme case of practically free rotation the 
selection rules can be readily obtained by introducing the circular coordi- 
nates described in the footnote on page 327. Using this new system of coordi- 
nates the tensor components a,,(qg,¢) of equation (3.14) (u,» = 0, +1) 
transform under rotations about the z-axis by representations with character 
exp[i(u—v+s)¢] and the matrix element [a,,(q, ¢)]%-m- has the character 
exp[i(u—v-+s+m’—m’’)¢]. This gives the rotational selection rule 


(5.4) m"’—m’' = p—v+s, 


according to which all tensor components a,,(g,s) may contribute to the 
first-order Raman effect. The components a,,(qg, s) do not correspond to the 
full symmetry of the crystal if s # 0 and the same is true for the selection 
rules to which they are subject. 

Summarizing one can say that the selection rules obtained from the linear 
polarizability theory may be violated by the third-order terms of equation 
(3.7) in crystals with degenerate normal vibrations, by second-order terms in 
crystals with strong mechanical anharmonicity, and by rotating groups. 
All these higher order effects strongly increase with increasing temperature. 
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Their magnitude should, however, be independent of the direction of the 
incident light.* 
(d) Selection Rules for the Scattering by Electrons in Delocalized Crystal Orbitals 
(i) Selection Rules for the Wave Vectors 
The selection rules for the scattering by electrons in delocalized crystal 
orbitals are determined by the transformation properties of the general 
scattering tensor [a,(q)]%» of equation (2.17). We investigate first the 
integrals 


(5.5a) {Ie (q)}iew = fexpexid’ 1) (pve) dv 


and 


ve 


(5.55) [R°@lr = J exp(—2niQ” .T)(pv,)} do, 
by expanding the space function (pv) }: 
, ’ 0 . 
(5.6) (pv) iv = [(ov) i+ Jo+ >. [2 owt | Qin 
j qi 0 


According to equation (2.16) [(pv)}-]o transforms like the product of two 
electronic wave functions of the crystal with nuclei fixed in the equilibrium 
positions; hence, [(pv)i-Jo transforms by a representation with character 
exp(27iq;.a@:) . (Qi: = Qv'—qv). Similarly, the first-order terms of the 
expansion (5.6) have the character exp[2x7(q.:+q,) . a;]. This has the effect 
that the product of the integrals (5.5a) and (5.55) which appears in equation 
(2.17) vanishes unless at least one of the following pairs of conditions (a) and 


(d) is simultaneously fulfilled : 
. a condition required for the non-vanishing of 
(5.7) ap . - the product of the zeroth-order terms, i.e. 
~* , selection rule for the Rayleigh scattering; 
(5.8) =n - ’ =. (b) condition required for the non-vanishing of 
arene ; the mixed product of a zeroth-order and a 
first-order term, i.e. selection rule for the 
Q’+a:+q; = 0 (a), , , 
(5.9) QO” +4, = 0 (6), first-order Raman effect; 
condition required for the non-vanishing of 
Q’+q:4+q, = 0 (a), the product of the first-order terms, i.e. 
Q”+q:+q, = 0 (d), selection rule for the second-order Raman 
effect. 


(5.10) 


By eliminating q,, the conditions (5.8) and (5.9) are transformed into the 
selection rules (5.1). 


*The latter statement may be wrong if the thermal equilibrium of the crystal is disturbed 
by the irradiation with an intense beam of light. This irradiation produces strong heating 
effects and possibly a selective excitation of lattice vibrations which may easily depend on the 
direction of the incident light beam. 
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(ii) Selection Rules for the Matrix Elements ape}o4+1 

One question of physical interest is whether or not the selection rules 
depend on the wave vectors Q’ and Q” of the incident and scattered visible 
light. In order to answer this question we consider the conditions (5.8) and 
(5.9) to be selection rules for the electronic wave vectors q;. We find in this 
way that the summation over /’ and /’” in equation (2.17) contains only terms 
for which q; = —Q’ and q; = —Q”. This means that the wave vectors Q’ 
and Q” determine the electronic orbitals which contribute to the polarizability 
[x,e]¢- in such a way that different orbitals are involved for different directions 
of the incident and scattered light; this may cause a strong anisotropy of the 
tensor components [a,.]; even if the wave vectors Q’ and Q” tend towards 
zero. 

We come now to the special selection rules for the [a,.]}. According to the 
selection rules (5.8) and (5.9), the product of the two matrix elements (5.5a) 
and (5.55) is only different from zero if (5.5@) contributes to it the zeroth- 
order matrix element [(pv-)}Jo and (5.55) the first-order matrix element 
[8/dq,(pv,) log) or the other way round. Hence, the product {I*’(q)}i 
x {7*"(g)}t" does not vanish only if [(pv.)iJo and [8/Aq,(ov,) Jo are si- 
multaneously different from zero. This condition leads to the required selec- 
tion rules on the basis of the following argument: Both [(pv,)}-]o and 
[8/dq,;(pv,) 7 Jo refer to the equilibrium position of the nuclei and, hence, trans- 
form by a representation of the space group of the crystal. In particular the 
zeroth-order term transforms like the product w,-v,wy and the first-order term 
transforms like the product w,v,gjwy and both products must be totally 
symmetric. (w, and wy, are electronic orbitals.) Since the velocity compo- 
nents v, (o = 1,2,3) transform like the coordinates x, (o = 1, 2,3), this 
condition gives the selection rule* 


[ape] oy-+1 = 0 unless x(x, qi) = x (xe), 


(5.11) x(x.) = symmetry character of x,, 
x(x.q,) = symmetry character of the product x, q;. 


This selection rule for the first-order Raman scattering by electrons in delocal- 
ized crystal orbitals is formally identical with the selection rule for the scatter- 
ing by localized electrons. There is, however, the essential difference that, 
as far as the properties of the delocalized space functions [(depe/8q;)09]o4+1 
are concerned, the normal vibrations g, can not be approximately identified 
with fundamental lattice vibrations. This may be formally derived from the 
q-selection rules (5.8) and (5.9) which associate different electronic orbitals 
with different lattice vibrations depending on their respective wave vectors. 
More physically speaking, localized atomic or molecular orbitals in crystals 
depend only on the nuclear arrangement in their immediate neighborhood ; 
hence, the polarizability of electrons in such orbitals depends only on the 
local distortions produced by lattice vibrations, which are, in the case of 


*This selection rule is only valid for the first-order derivative and neglects higher terms in the 
expansion (5.6) or (3.3). 
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internal vibrations, practically equal for all q,; (q;—> 0). The polarizability of 
electrons in delocalized orbitals, however, depends on the arrangement of all 
the nuclei in the lattice and this introduces a difference between long waves 
and the fundamentals. It follows that the selection rules for a vibration with 
wave vector q; = —Q are determined by the group of q; and not even 
approximately by the full crystallographic point group. 
6. THE CASE OF CALCITE 

We conclude the theoretical investigation of the Raman effect in crystals 
with an analysis of the Raman scattering of calcite which has been investigated 
most thoroughly and is typical for all cases where deviations from the predic- 
tions of the linear, atomic polarizability theory occur. We give first a summary 
of the relevant experimental results and discuss them later in terms of the 
preceding theory. 


(a) Experimental Results 

(i) The Resonance Condition 1(c) 

The absorption spectrum of calcite (Table I) has been measured by Pfliiger 
and Ornstein (Ornstein and Went 1935a; Pfliiger 1904). It is difficult to predict 
whether or not the weak absorption in the region near 4000 A is compatible 
with the resonance condition 1(c); hence, additional arguments will be used to 
clarify this question. 


TABLE I 
THE ABSORPTION SPECTRUM OF CALCITE 
rAinA 2310 2400 2800 3910 4910 


Absorption 
coefficient 1.00 0.58 0.16 0.092 0.078 





Perhaps the most convincing argument is obtained from an analysis of the 
intensities of Raman lines corresponding to different frequencies v’ of the 
incident light. Ornstein and Went (19350) find that the experimental results 
can be explained by the assumption that only one of the different resonance 
denominators [(v;--)?— (v’)?] in equation (2.17) is decisive; it corresponds to a 
strong absorption in the wave-length region around 1600 A. This absorption 
does not violate the condition 1(c). The same may be concluded from the 
intensity ratio of Stokes and anti-Stokes Raman lines (Ornstein and Went 
19355) and from the experimental fact that the scattering tensor is symmetric 
(Michalke 1938). 


(ii) The Temperature Dependence of the First-order Raman Effect 

Ornstein and Went (19350) find an increase in line width and a decrease of 
frequency and intensity of all the first-order Raman lines if the temperature is 
raised. All these effects are most pronounced with the low frequency lattice 
vibrations and practically disappear in the case of the internal vibration of 
the CO;~~-ion with frequency 1088 cm.~'. Most interesting for our discussion 
is the pronounced inverse temperature effect which is quantitatively given in 
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Table II by the ratio R, = Jea/Jovs(Jons = Observed intensity, J.4; = intensity 
calculated from the formula of the linear polarizability theory; the subscript 
v gives the frequency of the Raman line). 
TABLE II 


THE INVERSE TEMPERATURE EFFECT IN CALCITE 








Reso Rus Rioss 





Eee Riss 

90 — 0.75 7 0.99 
293 1.00 1.00 1.00 1.00 
376 1.30 1.12 1.20 0.99 
426 2.12 1.22 1.60 1.02 


(iii) The Selection Rules 

The experimental results concerning the selection rules for the first-order 
Raman effect in calcite have been most thoroughly discussed by Couture 
(1947, 1948). We repeat some of her results which are necessary for our 
analysis. The experimental situation is described by a coordinate system x 
(direction of incident light), y (direction of scattered light), z which is fixed 
in space and a coordinate system X, Y, Z fixed in the crystal, such that Z is 
the direction of the optical axis. This allows one to distinguish between three 
different relative orientations of crystal and radiation which may be charac- 
terized by the direction Z of the optical axis: 

Z\\x, Z\|y, Z\|z (|| = parallel). 


For each of these three cases one can measure the degree of depolarization 
p = I,/I, and the total intensity of Raman lines J = J,+/,. (J, and J, are 
the components of the scattered light with polarization vector in direction x 
and z.) p and J are proportional to the squares of the matrix elements 
[ape(9)]+341 for which we write, in short, ¢,,. If the €,, are represented in the 
coordinate system fixed in the crystal and if natural incident light is used 
one finds the following relations: 


Z| \x Z\ly Z\\z 

ii eaxterz exy tézx exytezx 

(6.1) éxyt+eéyy : €éyyteyz ; €yzt+ézz : 
I = evyyt+A, evyyt+A, esztA, 


2 2 2 
A= exy teyztezx- 


This scheme allows one to calculate the tensor components ¢,, from the experi- 
mental values of p and J. We discuss only the Raman line corresponding to the 
totally symmetric internal vibration of the CO;~~-ion with frequency 1088 
cm.~! since it is the only one in the calcite spectrum which shows pronounced 
deviations from the predictions of the linear, atomic polarizability theory. 
The experimental results are contained in the scheme 
Z| |x Z\ly Z\|z 
(6.2) p = 0.0, 0.4, 0.2, 
I = 1.4, 1.4, 0.3. 
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This gives 
(6.3) cow = 0.4, tin = das = 0| for the cases 
éyy = 1, Zilx and Z||ly, 
and 
2 _ > ee ai 
(6.4) xy = 008, ers ™ tx = Ol bor the case Z||z. 


éxz = 0.25, J 


These results show that the magnitude of some of the components of the 
scattering tensor depends strongly on the direction of the incident and 
scattered light and that exy is different from zero in contrast to the ordinary 
selection rules. It is, however, important to note that the cases Z||x and Z]||y 
are compatible with one single set of tensor components. 


(6) Interpretation of Experimental Results 


In order to explain the behavior of calcite in a convincing way, we exclude 
first all possible interpretations which are not really satisfactory and show 
later that the only remaining explanation is compatible with all the observable 
facts. 

The resonance condition 1(c) has been already dealt with in subsection 
6(a), (i) and has been shown to be reasonably well fulfilled. 

We come now to the higher order effects of subsection 5(c), (ii). These 
effects could, in principle, be responsible for the appearance of Raman lines 
which are forbidden by selection rules derived from the /inear polarizability 
theory. But the intensity of these forbidden lines should be strongly dependent 
on the temperature and independent of the direction of the incident light.* 
The contrary is true according to the experiments. 

The only remaining effect is the scattering by electrons in delocalized 
crystal orbitals which is most conspicuous if the scattering by electrons in 
localized orbitals is forbidden by the ordinary selection rules. It has been 
shown in the preceding sections that the magnitude of the components of the 
scattering tensor corresponding to electrons in delocalized crystal orbitals 
depends on the directions of the incident and scattered light. This explains 
qualitatively the most striking feature of the experiments. Furthermore, it 
can be shown that the results tabled in the equations (6.3) and (6.4) are 
compatible with the selection rule (5.11). The group of q = —Q is in all 
three cases the point group S which is formed by the two symmetry elements’ 
identity and center of symmetry. The normal coordinate qioss is totally 
symmetric and transforms by a small representation with character 1; the 
cartesian coordinates transform all by the same antisymmetric representation 
with character —1. Hence, no single tensor component e,, is forbidden and the 
finite experimental value of exy is compatible with the selection rules. The 
vanishing of eyz and €zx is not required by a rigorous selection rule, but is an 

"It has been mentioned in the footnote on page 332 that a selective excitation of lattice 
vibrations due to the strong irradiation with an incident light beam could be responsible for the 


anisotropy of the tensor components é. This possibility must be excluded for calcite, since 
the cases Z||y and Z||z are equivalent with respect to the direction of the incident light. 
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accidental property of the scattering system. The reason for the smallness of 
€yz and €zx is most probably the fact that the nuclear displacements corre- 
sponding to the totally symmetric vibration with frequency 1088 cm.~! are 
perpendicular to the z-axis; this leaves the matrix elements [ayz]j4:1 and 
[azx]o41 apparently unchanged and affects mainly the matrix element 
[axy]e41. It appears that simple relations between nuclear displacements and 
the corresponding polarizability changes of the type suggested above are most 
pronounced in the case of delocalized crystal orbitals; in fact, the experimental 
results of Couture-Mathieu (1950), 1952) indicate that they are in general 
true also for the polarizability of localized electrons. 

The equivalence of the cases Z||x and Z||y is also readily explained. These 
two cases can be formally transformed into each other by exchanging Q’ with 
Q”. This corresponds to the transformation Q— —Q which leaves the 
scattering system unchanged. 

Concerning the inverse temperature effect of the low frequency Raman lines 
of calcite, the electrons in localized and delocalized orbitals need not be 
considered separately. We believe that strong line broadening is the essential 
cause of this effect but we do not exclude contributions of the higher-order 
polarizability derivatives as a possible explanation. 
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LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any tssue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


The Resonant Absorption of ~-Rays by N" 


The Ce, y)N* reaction has a broad resonance ([ = 32 kev.) for 554 kev. protons (Seagrave 
1952) which forms N'tina 1—, T = 1 state of 8.06 Mev. excitation (Clegg and Wilkinson 1953). 
This state decays 90% by transitions direct to the ground state (Woodbury, Day, and Tolle- 
strup 1953). The possibility of exciting N™ nuclei to the 8.06 Mev. state by absorption of these 
8.06 Mev. y-rays is suggested by two facts: first, the great width of the level means that shifts 
in the y-ray energy, occurring between production and subsequent resonant absorption, due 
to Doppler and recoil effects are not as serious as in the resonant scattering work of Moon 
(1951); second the recoil losses occurring on emission and absorption and amounting to about 
5 kev. total can be compensated for by the Doppler shift which in effect transfers some of the 
incident proton energy to the y-rays which are emitted in a forward direction. 

The resonant absorption has been observed in the present experiment. A thick C™ target 
containing 65% C" (supplied by Dr. M. L. Smith, A.E.R.E., Harwell, England) was bom- 
barded with 700 kev. protons. The target was sufficiently thick that the protons did not leave 
it until reduced to about 420 kev. or less, so that with little error the y-ray yield corresponded 
to the full integral over the 554 kev. resonance. 

The inverse resonant absorption of the 8.06 Mev. y-rays in N“ was detected using a glass 
proportional counter, with an active volume 27 cm. long by 4 cm. in diameter, filled to a 
pressure of one atmosphere with nitrogen gas and operated with a gas amplification of 10. A 
metallic layer of sodium evaporated into one end of the counter maintained the purity of the 
gas. The energy scale was calibrated by means of the N"(n, p)C™ reaction produ in the 
counter gas by thermal neutrons and by injecting a-particles into the counter via a thin mica 
window. 





¥—RAYS 
TO SCINTILLATION 
COUNTER 


Fic. 1. The proportional counter and target used to observe the resonant N“(y, p)C™ reaction. 


The counter was placed adjacent to the C™ target as shown in Fig. 1 and the pulses from it 
were recorded on a 100 channel kicksorter (Wilkinson 1950). Simultaneously the y-ray flux 
was measured with a Nal scintillation counter in conjunction with a 60 channel kicksorter 
(Hutchinson and Scarrott 1951). The efficiency of this counter for 6 Mev. y-rays had been 
determined by simultaneous @ and y counting on F!*(p, a y)O"* (Van Allen and Smith 1941). 
The efficiency for 8 Mev. y-rays was obtained L extrapolation using the theoretical absorption 
cross sections. Typical proton and y-ray spectra are shown in Figs. 2 and 3. That the 550 kev. 
peak was due to resonant absorption of y-rays by the reaction N"(y, p)C™ was confirmed by 
the absence of a peak in the spectrum when the counter was bombarded with 6, 12, or 17 Mev. 
y-rays. Calculations are in progress to determine the N*(y, p)C™ cross section from the data; 
results and experimental details will be published later. 

I would like to thank Dr. D. H. Wilkinson who suggested this experiment to me and Prof. 
O. R. Frisch who extended to me the hospitality of the Cavendish laboratory where this work 
was done. I am also indebted to Mr. D. B. Stewart who operated the Cavendish high tension 
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Fic. 2. Proportional counter spectra produced by 1.12 Mev. test a-particles (right-hand curve) and by the 
N'4(y, p)C® reaction. 
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Fic. 3. Typical pulse height spectrum produced by the C"(p, 7)N™ 4-rays. 


set. The work was carried out while the author held a Royal Society of Canada Rutherford 
Memorial Fellowship and a National Research Council of Canada Fellowship, for which he is 
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BUCKMASTER: TRANSITIONS 341 


4M = +3 Transitions in Dilute Gadolinium Ethyl Sulphate 


The seven AM = +1 and six AM = +2 paramagnetic resonance transitions in dilute 
gadolinium ethyl sulphate have been studied by Bleaney, Scovil, and Trenam (1954) and 
Buckmaster (1956) respectively. The sensitivity of the new spectrometer developed by Buck- 
master and Scovil (1956) is sufficient to also allow the AM = +3 transitions to be observed 
at 90° K. with signal-to-noise ratios of about 10: 1. Unfortunately, this sensitivity is insufficient 
to observe these transitions within 25° of the parallel direction. The value of the spectroscopic 
splitting factor (g) for the perpendicular direction is given in this letter along with a com- 
parison of the calculated and observed positions for the same direction. 

The AM = +2 and +3 transitions both arise from the same admixture of wave functions 
describing the energy levels due to the largest off-diagonal term P3(S,) in the spin-Hamiltonian 
(Elliott and Stevens 1953). It is probable that other off-diagonal terms such as Pi(S;) are also 
important; however, they do not improve the theoretical agreement with the observed values 
given below. The transition probabilities for the AM = +3 transitions are smaller than those 
for the AM = +2 transitions by the same factor that the AM = +2 are smaller than the 
AM = +1. The observed intensities of the AM = +1, +2, and +3 transitions are in the ratio 
10‘: 10?: 1 confirming this. Four of the five AM = +3 transitions are shown in Fig. 1. 





Fic. 1. Four of five AM = +3 transitions near crossover point. 


The spectroscopic splitting factor for the perpendicular direction is given in Table I. 


TABLE I 
v = 24.520 kMc./sec.. T = 90° K. 








Transition Mean field Correction Corrected 

4M = +3 (gauss) factor (P2) field gi 
+7/2<> +1/2 2932 +4 2936 1.9886 
+5/2¢5 +1/2 2897 +41 2938 1.9875 
+3/2<> -3/2 2893 +54 2947 1.9814 


Mean value g,: uncorrected 2.0085 +0.001, corrected 1.9858 +0.001 
Mean value g,: corrected excluding +3/2 <> —3/2 transition 1.9881 +0.001 


The values of gi corrected for (P3}(S,)) obtained with the same crystal from the AM = +1 
and +2 transitions (Buckmaster 1956) at 90° K. are 
AM = +1, g, = 1.9901+0.001; 
AM = +2, g, = 1.9901+0.001. 


This shows that the effect of other off-diagonal terms such as P3(S,) is becoming important. 

It is not possible to calculate the splitting parameters from the AM = +3 transitions. The 
observed positions for the perpendicular direction are compared in Table II with the positions 
calculated from the values of the splitting parameters obtained from the AM = +1 and +2 
transitions in this direction (Buckmaster 1956). The agreement is not within the experimental 
error. 








TABLE II 
v= 5 520 | kMc. ‘sec., T = 90° K. 

Transition Observed position Correction Corrected Calculated Calculated 

4M = +3 (error +1 gauss) factor (PR) position AM = +1 data AM = +2 data 
+7/2 <> +1/2 3364 + 1.7 3366 3374 3369 
+5/2¢> -1/2 3101 +38.8 3140 3141 3138 
+3/2¢> -3/2 2893 +53.8 2947 2934 
+1/2 49 —5/2 2692 +44.2 2736 2731 
-—1,2¢5 -7/2 2500 + 6.8 2507 2500 
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Recently Cooke (1955) has reported that a resonance from an excited state in dilute dys - 
prosium ethyl sulphate has been observed with g,, = 5.80+0.02 and g, = 8.40+0.2 which 
disappeared below 13° K. because of suggested depopulation and above 20° K. because of 
spin-lattice broadening. A weak resonance (S/N = 10:1) has been observed at 4.2° K. with 
gi, = 5.85+0.05 using the high sensitivity spectrometer (Buckmaster and Scovil 1956) in 
agreement with the above value and confirming the depopulation hypothesis. 
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